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This reference manual contains both technical and programming guides for the expert system, 
Thermal Analyst's Help Desk. The Thermal Analyst's Help Desk provides the user with a capability 
for performing the preliminary analyses required to determine the thermal characteristics of a 
spacecraft. It gives the analyst a first look at the parameters required and provides an assessment 
of the feasibility of a preliminary thermal system design. Analyses can be performed on a single flat 
plate or a rectangular box enclosure. The enclosure analysis can be performed at a single orbit point 
or at a user defined number of orbit points. The orbit itself may be any elliptical conic section around 
any planet. 

This reference manual is intended to give the knowledge engineer insight into the theory and 
development of the Thermal Analyst’s Help Desk. It is not the purpose of this document to describe 
capabilities and functions that are contained in the EXSYS® manual. Therefore, there are locations 
in this document that are lacking in information. In these locations there are references to the 
EXSYS® manual and other readings about expert systems. It is not intended to replace those 
references in this manual. In certain instances, EXSYS® options are described more fully to give the 
knowledge engineer further insight into specific aspects of EXSYS® that are used in Help Desk. 
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INTRODUCTION 


The following document describes the background and implementation of the theory used to create 
the Thermal Analyst's Help Desk. 

The Thermal Analyst’s Help Desk (also referred to as Help Desk) was created as a prototype to 
determine if expert systems can be used in the design and analysis of spacecraft. Its function is 
to provide first order approximations to thermal design parameters subjected to a radiation 
environment. 

Its features include: 

• sizing a radiator with a known surface temperature 

• finding the temperature of a surface with a known area 

* finding the sink temperature of a surface with only environmental effects 

• determining the temperature and heat loads of an enclosure that is affected by environmental 
radiation 

* determining the temperature and heat loads on an enclosure as it proceeds around an orbit 

* assistance in determining the required emissivity and absorptivity for a surface. 

The Thermal Analyst's Help Desk runs within an expert system shell called EXSYS® from EXSYS, 
Inc. in Albuquerque, New Mexico. 

This document is divided into two parts: 

1) Technical guide: describes the equations and theory that pertain to Help Desk 

2) Programmer's guide: describes features about Help Desk that are not explained in the 

EXSYS® user’s manual 


Technical Guide 
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PARTI: TECHNICAL GUIDE 


1 EXPERT SYSTEMS 

An expert system uses artificial intelligence (AI) techniques that allows the computer to perform 
as an expert in a specific and well defined field. 

An Expert system has three basic parts, 

1) a collection of knowledge on a given topic (knowledge base) 

2) a set of rules pertaining to the knowledge (rule base) 

3) a method to apply the rules (inference engine). 

The knowledge and the rules can pertain to any topic that requires an expert to determine a 
result. The Thermal Analyst’s Help Desk is an expert system that is used to help determine the 
thermal characteristics of a spacecraft. 

In this text the expert is the source of the expertise which is included in the knowledge and rule 
bases. The knowledge engineer is the developer of the expert system, in this case the Thermal 
Analyst's Help Desk. The user refers to the individual running an analysis. 

A knowledge base contains a series of symbolic parameters, called qualifiers, and possible valid 
responses, called values. The knowledge base also contains the definitions of variables which 
include variable name, type, numeric bounds and initialization values. The Help Desk knowledge 
base contains information about spacecraft and thermal analysis. Qualifiers and their values are 
created by the knowledge engineer using information known about the system. For example, a 
sample knowledge base is shown in Box 1. It shows the qualifiers and their values for the Surface 
Properties Analysis in Help Desk. Qualifier 1, "Value to calculate is" has six values listed. 
Qualifiers 9, 10 and 13 each have two values. These qualifiers do not assume an application or 
an order, they are simply a collection of knowledge. 


/* Qualifier 1 

/♦ Qualifier 10 

Q> Volue to calculate is 

Q> the known surface information is 

V> rodiator area (given temperature) 

V> surface use 

V> surface temperature 

V> surface material 

V> enclosure temperatures and heat loads 
V> orbital average temperatures and heat loads 

/* Quolifier 13 

V> surface properties 
V> none of the above 

Q> Surface properties are 


V> calculated 

/* Qualifier 9 

Q> the performance time is 

V> unknown 

V> 80L (beginning of life) 
V> EOL (end of life) 

- 


Box 1 Sample Knowledge Base 


A rule base is a set of IF-THEN rules assigning values to qualifiers and numbers to variables. 
The IF-THEN rules in any expert system are comparable to IF-THEN statements in most higher 
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level computer languages. An example comparison is shown in Box 2 using the qualifiers and 
values from Box 1. A sample rule base, also from the Surface Properties analysis in Help Desk 
is shown in Box 3. The premise clauses (IF statements) for the four rules are similar. They 
represent all possible combinations of the two qualifiers. If one premise clause is false, the rule 
is marked False and is never executed. When all the premise clauses in a rule are true, the rule 
is marked True and is a candidate for execution by the expert system. Since any rule can be true 
at anytime, the rule base defines a state-it does not imply an order. The rules are ordered to 
increase efficiency, but not to define a sequence for evaluation. 



EXPERT SYSTEM 


FORTRAN 

RULE NUMBER: 39 (FRAME EPS BOL) 



IF: 

Value to calculate is surface properties 

IF 

((01 .EQ. VI) .AND. 

and the performance time is BOL (beginning of life) 

(09 .EQ. VI) .AND. 

and the known surface information is surfoce use 

(010.E0.V1)) 

THEN: 


THEN 


[ 

EPSILON] IS GIVEN THE VALUE.... 


EPSILON] = ... 

and 

ALPHA] IS GIVEN THE VALUE .... 


ALPHA] = ... 

and 

MATERIAL] IS GIVEN THE VALUE 


MATERIAL^ ... 

and Surface properties are calculated 


013 = VI 



ENDIF 



Box 2 FORTRAN Equivalent Statements 


The inference engine determines which rules are used, and the order in which they are applied. 
The inference engine defines the order by "firing" the rules. A rule is "fired" when the THEN 
clause is performed. The selection of applicable rules can be accomplished by multiple means. 
This expert system uses the method of chaining, which is explained in detail in the next section. 
In the example, this process reduced the examination of all possible rules to the examination of 
just four rules shown . When the inference engine examines a premise clause it asks the following: 

Is the value defined? 

yes: Is the premise true or false 

no: can the value be derived from an additional rule? 

yes: derive the value 

no: prompt the user for the answer. 

Figure 1 shows a flow diagram of the procedure described. 

This method is applied to the example shown in Box 3. The inference engine would start with 
Rule 39 because, when there is no implied order in the chaining process, the inference engine 
evaluates the rules sequentially. For qualifier numbers and values, see Box 1. 
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RULE NUMBER: 39 (FRAME EPS BOL) 

IF- 

RULE NUMBER: 41 
ir. 

ir . 

Value to calculate is surface properties 

Ir . 

Value to calculate is surface properties 

and the performance lime is BOL (beginning of life) 

ond the performance time is EOL (end of life) 

and the known surface information is surface use 

and the known surface information is surface 



materiol 

THEN: 




[ 

EPSILON] IS GIVEN THE VALUE.... 

THEN: 


and 

ALPHA] IS GIVEN THE VALUE .... 

[ 

EPSILON] IS GIVEN THE VALUE .... 

and 

MATERIAL] IS GIVEN THE VALUE 

and 

[ALPHA] IS GIVEN THE VALUE .... 

and Surfoce properties are calculated 

and 

APPLICATION] IS GIVEN THE VALUE 



and Surface properties are calculated 

RULE NUMBER: 40 (FRAME EPS EOL) 



IF: 


RULE NUMBER: 42 

Value to calculate is surface properties 

IF: 


and the performance time is EOL (end of life) 

Value to calculate is surfoce properties 

ond the known surface information is surface use 

and the performance time is BOL (beginning of life) 



and the known surface information is surface 

THEN: 


materiol 

[ 

EPSILON] IS GIVEN THE VALUE .... 



and 

ALPHA] IS GIVEN THE VALUE .... 

THEN: 


ond 

MATERIAL] IS GIVEN THE VALUE .... 

[ 

EPSILON] IS GIVEN THE VALUE.... 

and Surface properties are calculated 

and 

ALPHA] IS GIVEN THE VALUE.... 



ond 

APPLICATION] IS GIVEN THE VALUE.... 



ond Surface properties are calculated 


Box 3 Sample Rule Base 


EXAMPLE: 

Qualifier 1: "Value to calculate is" 

Is the value defined? 

no: con its value be derived? 

no: prompt the user with a screen for the information. 

When running Help Desk, the user is prompted ot this point with the Main Menu, the menu 
with the analysis choices. If the user selects "surfoce properties", the inference engine 
ossigns "surfoce properties" to qualifier I (or V5 to Q 1 ) ond proceeds to the next premise 
clause. 

NOTE: In Help Desk, once o qualifier hos been ossigned o volue, it cannot be chonged. 

Qualifier 9: "the performance time is" 

Is the value defined? 

no: can its volue be derived? 

no: prompt the user with a screen for the information. 

The user is prompted for o choice of performance times: Beginning of Life or End of Life. 
If the user selects Beginning of Life, the inference engine ossigns VI to 09 ond proceeds to 
the next premise douse. (If the user selects End of Life, the inference engine ossigns V2 
to 09, marks the rule os Folse ond proceeds to the next rule.) 
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Qualifier or Variable 


Is the value defined? 




Proceea to next 
premise clause 
or mark rule True 


Mark rule False 
proceed to next 
rule 


Figure 1 Inference Engine Decision Diagram 


Ouolifier 10: "the known surface information is" 

Is the volue defined? 

no: con its volue be derived? 

no: prompt the user with o screen for the information. 

At this point, the user is prompted for the known surface information. If the user selects 
surface use, the inference engine assigns VI to 010 and marks the rule as True. (If the user 
selects surfoce material, the inference engine ossigns V2 to 01 0, marks the rule os false and 
proceeds to the next rule.) 

When a premise is completely satisfied the inference engine "fires" the conclusion. In this case it 
defines the value for EPSILON (emissivity), ALPHA (absorptivity), determines the unknown surface 
information, ond defines 013 to have the value "calculated". 

This sample analysis can be seen by running the "surface properties” option from the main menu 
of Help Desk. 

1.1 Chaining 

The Help Desk inference engine uses a deductive procedure called chaining. Chaining is generally 
one of two forms or a combination of the two forms: 

1) forward 

2) backward 
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Forward chaining starts with a premise to prove a conclusion. The user is asked questions until 
all the applicable information has been derived from the expert system. The inference engine 
flows through all the rules sequentially firing all true rules. In this manner it determines a set 
of all possible conclusions. 

Backward chaining starts with an assumed conclusion and shows that the premise is true or false, 
proving or disproving the assumption. In this form the user specifies a conclusion. The inference 
engine then tries to prove the premise of the rule containing the conclusion. It does this by finding 
other rules in the rule base with conclusions that will verify the premise under inspection. The 
inference engine then tries to prove the premise of the new rule. The rules are thus chained 
together until the initial premise is proven or disproved. This method is similar to the AI method 
of recursion. Box 4 contains a set of rules taken firom Help Desk that are chained together in this 
manner. 

The rules displayed are used in the environmental temperature calculations in Help Desk. The 
rules shown are only those that will be fired by the expert system when running the analysis. 
There are more rules available that will calculate the qualifier value for "Thermal analysis is" and 
a value for the variable [EPSILON]. To simplify the example, only the rules that will be fired are 
shown. 


A RULE NUMBER: 1 

IF- 

A RULE NUMBER: 10 

IF- 

Ir . 

Thermal analysis is jcompleted| 

Ir . 

Volue to calculate is [surface temperature! 


and: 

EB] > 0 

THEN: 

and: 

[SIGMA] > 0.0 

X> REPORT(RADIATE.OUT) 


ond: > Thermal anolysis is complete - 

THEN: 


Confidence=l 


[TEMP] IS GIVEN THE VALUE ... 

A RULE NUMBER: 5 

IF* 

A RULE NUMBER: 12 
ir- 

Ir . 

Value to calculate is jsurfoce temperature! 

ir . 

Surface temperature to be calculated is 

ond: Surface temperature to be calculated is 

[environmental temperature! 

{environmental temperature! 

and: 

EPSIL0N]>0 

ond: [TEMP] > 0 

THEN: 

THEN: 

[EB] IS GIVEN THE VALUE ... 

[REPORT] IS GIVEN THE VALUE [REP0RT3] 

and: Thermal analysis is {completed} 

A RULE NUMBER: 22 

IF- 


ir . 

Value to calculate is [radiator area (given 


temperature)! 


or: Value to calculate is [surface temperature! 

- 

THEN: 

[SIGMA] IS GIVEN THE VALUE... 


Box 4 Example of Backward Chaining 
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In this example the user wants to prove that the "Thermal analysis is complete" (conclusion of rule 
1). The inference engine tries to prove this conclusion using backward chaining as defined above. 
It starts with Rule 1, which contains the conclusion in question. 


EXAMPLE: A decision tree for this example is shown in Figure 2. 
Surfcce 
temperature 



Figure 2 Decision Tree 


Rule 1: 

Qualifier: "Thermal anolysis is" 

Is the volue defined? 

no: can the volue be derived from on additional rule? 

yes: derive using Rule 5. 

Rule 5: 

Qualifier: "Value to calculate is" 

Is the volue defined? 

no: can the value be derived from on additional rule? 

no: prompt the user. 

Help Desk prompts the user with the Main Menu, which contoins o list of the analyses. If the user 
selects "surfoce temperoture", the inference engine assigns the value "surface temperature" to 
Q1 and proceeds to the next premise clause. 
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Qualifier: "Surfoce temperate to be calculated is” 

Is the value defined? 

no: can the value be derived from an additional rule? 

no: prompt the user. 

Help Desk prompts the user with o screen containing two choices: 

1) Surfoce temperature (given area) 

2) Environmental temperature 

If the user selects "Environmental temperature" the value is assigned to the qualifier and the 
inference engine proceeds with the next premise clause. (If the user selects "Surface temperature 
(given area)" the value is assigned to the qualifier, the rule is marked False, and the inference 
engine proceeds to the next rule that con satisfy the previous premise.) 

Varioble: [TEMP] (the temperature of the surface) 

Is the volue defined? 

no: can the value be derived from on additional rule? 

yes: derive using Rule 10. 

Rule 10: 

Qualifier: "value to calculate" 

Is the volue defined? 

yes: Is the premise clause true or false? 

true: continue with next premise clause 

Varioble: [EB] (blackbody radiation) 

Is the volue defined? 

no: can it be derived from on additional rule? 

yes: derive using Rule 12 

Rule 12: 

Qualifier "Surface temperature to be calculated is" 

Is the value defined? 

yes: Is the premise clause true or folse? 

true: continue with next premise clause 

Variable: [EPSILON] (surface emissivity) 

Is the volue defined? 

no: con it be derived from on additional rule? 

no: prompt the user 

— -This user js-prompted with a screen -that contains a location to input a value -for- 

[EPSILON], The screen also contains positions for the value of other anolysis constants. 
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If the volue for [EPSILON] is greater thon zero, the rule is marked os True. The THEN clause 
is fired, defining the volue for [EB], 

Note: [EPSILON] < 0 hos no physical meaning. .Help Desk will continue to prompt for 

this variable until it is properly satisfied. This is due to the variable definition, not 
the inference engine. 

(A tye 10 continue# 

Once the inference engine is finished with a rule, it returns to the colling rule. Since Rule 10 
colled Rule 12. control now returns to Rule 10. 

Varioble: [SIGMA] (Stefan Boltzmann constant) 

Is the volue defined? 

no: con it be derived from an additional rule? 

yes: derive using Rule 22 

Rule 22: 

Qualifier: "value to calculate is" 

Is the volue defined? 

yes: Is the premise douse true or folse? 

true: mark the rule os True and fire the THEN douse. 

Rule 22 contains on "or" douse. The first clause is folse, but the second "or" clause 
is true, moking the entire premise true. The THEN douse ossigns the value to [SIGMA], 

Currently EXSYS* does not handle scientific notation. Therefore, a varioble cannot be 
initialized to such o small volue, ond consequently the value must be assigned through 
the use of o rule. 

( Rule JO continued] 

All the premise douses for Rule 10 are defined ond True. The inference engine morks the rule 
os True ond fires the THEN douse, which defines the volue for [TEMP], 

fRu/e 5con/inued) 

All the premise douses for Rule 5 ore defined. If the variable [TEMP] > 0, then the entire premise is 
True. The inference engine morks the rule as True ond fires the THEN douse. The THEN clause defines 
the variable [REPORT] (which displays a screen containing the report information) and ossigns the 
volue "completed" to the Qualifier "Thermal onolysis is". 

flfu/e I continued) 

The premise for Rule 1 is derived and True. The inference engine marks Rule 1 as true and fires the THEN 
douse. The THEN douse executes the report generator found in the file RADIATE.OUT, ond defines "Thermal 
Analysis" os "complete", proving the conclusion. 

This brief description of Help Desk and the inference engine should be sufficient for a cursory 
examination for the knowledge and rule bases of Help Desk. For more information on expert 
systems and chaining, refer to the text in References 3,4 and 5. 
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2 MATHEMATICAL APPROACH 


The following section covers the derivations of the equations used in Help Desk. 
2.1 Environmental Heat Balance on a Surface: 



The net heat flux through the surface, Qj, equals the sum of the environmental effects on the 
surface minus the heat radiated from the surface, Qrad- The environmental effects include the 
radiation on the surface due to the sun, Q SUN ; the planet, Q PLANET ; and the albedo of the planet, 
Qalbedo- 


QrQsUN*QpLAMT*QALBEDO~QllAD 


Equation 1 


The surface of a radiator must be able to remove all of the environmental heat loads plus any heat 
load generated by internal equipment, Qinternal- An example of an internal heat load is the energy 
that must be removed to cool an instrument to a defined temperature. The addition of this term 
makes the above equation: 


QcQsm+Q f u^ + QAiJUM + QimE*M-Q*AD Equation 2 

For a single surface, Qj is zero when the surface is in equilibrium. 

2.2 Surface Heat Loads 

Q sun is the heat load on the surface due to the radiation of the sun as defined in Equation 3. It 
is assumed that the solar rays are essentially parallel in the vicinity of a planet, and consequently 
solar flux is constant throughout an orbit (except for planetary shadowing effects). However, 
variations in flux that are due to yearly variations in the distance from the sun are accounted for 
with the user defined solar flux value. Help Desk also calculates changes associated with the 
surface properties. 

QsmTAfa a* Equation 3 

oc, solar absorptivity Fraction of the total solar energy (in the 0.2 to 2.8 pm range) that is 

incident on a surface and is absorbed by that surface. 

S Solar flux, W/m 2 The radiative flux due to the sun. Varies as the inverse square of the 

distance from the sun, which changes throughout the a year. 
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solar view factor Fraction of the surface field of view that sees the sun due to variations 

in the angle of incidence of the surface. See Section 5.2, Figure 7 for 
more information. 

Aj surface area, m 2 Area of the surface under investigation. 

Q planet * s h eat l° a d on the surface due to the thermal radiation from the planet as defined in 

Equation 4. The surface of a planet absorbs energy from the sun and emits it as thermal radiation 
in accordance with the Stefan-Boltzmann law. It is assumed that thermal radiation is emitted 
from a planet at a constant rate. However, it does vary with the surface infrared absorptivity, the 
area of the surface, and the view factor from the surface to the planet. 

Qpunet= a FeP«ir Equation 4 

E Earth flux, W/m 2 Radiative flux of the planet. A portion of the solar radiation incident 

on the planet is reemitted as thermal radiation. 

F Ei planet view factor: Fraction of the surface that sees the radiation from the planet. For 

example, in Low Earth Orbit, fully facing the planet, F E - .88, due to 

the surface seeing space off to the sides of the field of view. See 

Section 5.3 Figure 11 for more information. 

a IR Infrared absorptivity Fraction of the total energy in the infrared wavelength band (5 to 50 

pm) incident on a surface, that is absorbed by that surface. Assuming 
a gray body, this is equal to the emissivity of the surface. 

Aj surface area, m 2 Area of the surface under investigation. 

Qaibedo * s the heat load incident on the surface due to the solar energy that is reflected off the 

surface of a planet as defined in Equation 5. The reflection is due to the planet's surface and the 

atmospheric conditions. It is assumed that the albedo is constant over the planet's surface. It is 
a function of the surface infrared absorptivity, the area of the surface, the view factor from the 
surface to the planet, and the view factor to the sun. 

Qalbedo sA ^-^ ALB a s Equation 5 

ALB albedo coefficient Fraction of the total solar radiation that is reflected off the planet's 

surface. Depends on the atmospheric conditions of the planet. For 
Earth it ranges from .1 to .8 with and average of .3. 

E Earth flux, W/m 2 Radiative flux of the planet. A portion of the solar radiation incident 

on the planet is re-emitted as thermal radiation. 

S Solar flux, W/m 2 The radiative flux due to the sun. Varies as the inverse square of the 

distance from the sun, which changes throughout the year. 

(1-Pi) view factor Approximation for the reflected view factor. 

Fbj planet view factor Fraction of the surface that sees the radiation from the planet. For 

example, in Low Earth Orbit and when fully facing the planet, F E - .88. 
This is due to the fact that the surface sees space near the edges of its 
field of view. See Section 5.3, Figure 11 for more information. 
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a IR Infrared absorptivity Fraction of the total solar energy in the infrared wavelength band (5 

to 50 pm) incident on a surface that is absorbed by that surface. 
Assuming a gray body, this is equal to the emissivity of the surface. 

Aj surface area, m 2 Area of the surface under investigation. 

Qinuimaiis the internal heat load that must be removed from the system. Examples of the internal 
heat load are: 

• The electrical power of the internal systems that must be removed from the system 

• The heat load on an instrument due to an enclosure. This can be calculated using the 
Enclosure analysis discussed later in this document. 

Q rad is the heat radiated from the surface as defined in Equation 6. It depends on the blackbody 
radiation of the surface, the surface properties, and the surface area. 

<?«» ■ Afflrt . Affi, Equation 6 

Emitted power in W/m 2 for a body that absorbs all incident radiation. 
Temperature of the surface in K. 

= 6.78E-8 


The ability of a surface to radiate energy compared with a blackbody. 
For example, the emissivity of a blackbody is 1. Assuming a gray body 
this is equal to the infrared absorptivity of the surface. 

Area of the surface under investigation. 

2.3 Help Desk Surface Analyses 

In this subsection a single surface is considered, and the i subscript is omitted 
2.3.1 Surface Area 

The surface is in equilibrium. Therefore, the net heat load through the surface is equal to 
zero, Q = 0. 

Solve Equation 2 for the surface area, A: 

0 = ASpctg + AF E Ea 1R +AF E (l-p)S ALB ctg +Q internal * AeE b 


E bi blackbody radiation 

Tj temperature 

a Stefan-Boltzmann 
constant 

e emissivity 
Aj surface area, m 2 


. (~QqrTERNAl) 

(5p cc s +F gEttm+F S ALB 


Equation 7 


Equation 7 is the basis for the "Surface Area (Given Temperature)” option of the Thermal 
Analyst's Help Desk. 

2.3.2 Surface Temperature 

The surface is in equilibrium. Therefore, the net heat load through the surface is equal to 
zero, Q = 0. 
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Solve Equation 2 for the surface temperature, T (or the blackbody radiation, E b since E b = oT*): 


0 — ASpotg + AF E Ea IR +AF g(l-p)S ALB ctg +Q internal ■ AeE b 
E b = aT 4 


V 


Spa s +F^a a +F/1-»)S AIA 


Equation 8 


Equation 8 is the basis for the ’’Surface Temperature (given Area)" option of the Thermal 
Analyst’s Help Desk. 

2.3.3 Environmental Temperature 

When there are only environmental effects on the surface due to the sun, planet and albedo, 
the surface is in equilibrium. The total heat load through the surface is equal to zero, Q = 0. 
Solve Equation 2 for the surface temperature, T (or the blackbody radiation, E b since E b = oT 4 ): 

0 = ASpotg + AFEEa^+AFgd-pJS ALB otg - AeE b 
E b - aT 4 

ALB <*$) Equation 9 

Equation 9 is the basis for the "Surface Temperature (given Area)” option of the Thermal 
Analyst's Help Desk. 

2.4 Enclosure Heat Loads 

2.4.1 Heat Balance for Each Surface of an Enclosure 



Figure 4 Enclosure Geometry and Surface Numbers 


Expand Equation 2 for each surface of a six sided rectangular box shown in Figure 4. 


1: South facing 
2: Right 
3: Top 
4: Left 
5: Bottom 
6: North facing 


Qi = A 1 Sp 1 Og 1 + AjF^Eocjkj + AjF^fl-pjlS ALB ctgj -A^E^ 
Q 2 = ^2Sp 2 ctg 2 + A^^gEct^ + A,F 2E (l-p 2 )S ALB o tg2 -A^E^ 
Qg = AgSpgOtgg + AgFggEOtjjjg + AgFggd'PslS ALB Ctgg 'AgCgE^g 
Q 4 = A 4 Sp 4 o S4 + A 4 F 4E Ea IR4 + A 4 F 4E (l-p 4 )S ALB a* -A^E^ 
Qs = AgSpgOtgs + AgFggEajjjg + AgFggd-p^S ALB ot^ -A 5 e 5 E b5 
Qe = ^eSpgOtgg + A^Ea^ + A^d-PelS ALB Ot^ -A^E^ 
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In matrix form: 


0, 


Qsuni *Qeakthi +Qalbi 


~XraDI 

0 

0 

0 

0 

0 


hi 

0 2 


QsUN2 + QeaKTH2 + QaLB2 


0 

_ Y 
A FAD2 

0 

0 

0 

0 


hi 

03 


Qsum ^earths *Qalb3 


0 

0 

~^RAD3 

0 

0 

0 


F b3 

04 


QsUN4+QeaB1H4*QaLB4 


0 

0 

0 

~XraD4 

0 

0 



05 


Qsuns +QeAKTH5 *Qalbs 


0 

0 

0 

0 

~*RADi 

0 


hs 

06 


QsUN6 +QeAKTH6 +QaLB6 


0 

0 

0 

0 

0 



E * 


Equation 10 


or: 

Q = X+XRAD ~E b 


Equation 11 


Where : 

X contains the environmental effects. 

XRAD contains the area times the emissivity 

2.4.2 Enclosure Heat Balance 

In Section 2.4.1 the heat balance is performed on the external surfaces of the enclosure. In 
this section the heat balance is performed on the internal surfaces of an enclosure. It is 
assumed that the temperature is constant across an enclosure wall. 

An enclosure heat balance is the net heat load on an internal surface of an enclosure due to 
the other surfaces of the same enclosure, and to itself if it is concave. A heat balance on single 
surface i yields Equation 13. For the enclosure there are six such equations having values for 
i equal to one through six. The heat balance includes energy arriving on surface Aj from all 
surfaces of the enclosure, including A ^ if it is a concave surface. 

NOTE: The following equations include the capability for a concave surface. However, 

Help Desk currently contains only a rectangular box enclosure. Therefore, there 
are no concave surfaces. 


N (X 1 -p 

E t'-'.vT' 

J ml e J *J 


^"E (h~ F i-j) al j 
A j J m i 


Equation 12 


in a different form: 


N 


Qr°T?A r Y, 

J‘ i 


Equation 13 
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which simplifies to: 




Qr A iZ{ E „rW‘-J 


Equation 14 


Expand for each of the six surfaces: 

Qi = + (Eb^E^JF^ + (Ebi-E^F^ + (E bI -E b4 )F 14 + (E bl -E b5 )F 15 + (E bl -Ebg)F 16 ] 

Q 2 = A2[(E b2 -E bl )F 21 + (E b 2*Eb2)F 22 + (E b2 -Eb3)F 23 +■ (E^-E^JF^ + (E b 2"E b5 )F 2 5 + (E b 2‘Ebg)F 26 ] 

Q 3 = AgKEbg-Eb^Fa! + (Eb3“Eb 2 )F3 2 + (E^-E^JF^ + (E b3 -E b4 )F 34 + (E b 3-E b5 )F 35 + (E^-EbgJFsg] 

Q 4 = A 4 [(E b4 -E bl )F 41 + (E b4 -E b2 )F 42 + (E^-E^JF^ + (E^-E^F^ + (E b4 -E b5 )F 45 + (E b4 -E b6 )F 4€ ] 

Q s = AstCEbs-Eb^F^ + (Eb5-E b2 )F 52 + (E^E^F^ + (Ebs-E^JFg, + (E b5 -E b5 )F 55 + (E b5 -E b6 )F 56 ] 

Qe = AgtXEbg-Eb^Fg! + (E b 6*E b2 )Fg 2 + (Ebg-E^Fgg + (E b6 -E b4 )Fg 4 + (E b6 -E b5 )Fg 5 + (Ebg-Eb^Fgg] 


where 

E bi is the blackbody radiation intensity for the surface 


F 


is the view factor from surface j to surface i 


~Eb5E 25 


J b6 r 16 

-E^F, 


‘EbgFgg 


■Ev^Fo 


J b6 x 26 


Remove the zero terms (F u coefficients) and collect coefficients for the E bi terms: 

Qi = A 1 [E bl (F l2 +F 13 +F 14 + F 15+ F 16 ) -E b2 F 12 -E^F^ -E b4 F 14 -Eb5F 15 -E>^F 

Q 2 = A2[-E bl F 21 +E b2 (F 21 +F 23 +F 24 + F^+Fjg) -E b3 F 23 -E^F^ 

Q3 = A 3 [-E bl F 31 *E b2 F 32 +E b3 (F 31 +F 32 +F 34 + F 3S +F 36 ) -I 
Q4 = A 4 [-E bl F 41 -E b2 F 42 -EbgF^ +E b4 (F 41 +F 42 +F 43 + F^+F^) -E^F* -E^FJ 

Qs = ^st'^bi^si ~E b2 F 52 'E b3 F 53 -E b4 F 54 +E b5 (F 51 +F 52 +F 53 + F 54 +F 56 ) -EbgF 5 g] 


^b6 r 36^ 

-E^F, 


Qg - Ag[-E bl F 61 -E b2 Fg2 


'Eb3F 5 3 

"EbaFga 


-Eb4F 54 

*E b4 F 64 


'E^F^ 


+ E b g(Fg 1 +Fg 2 +F 63 + Fg 4 +F 65 )] 


In matrix form: 


f/4, 0 0 0 0 0 

0^ 0 0 0 0 

0 0 A 3 0 0 0 

0 0 0 A 4 0 0 

0 0 0 0 a 5 0 

0 0 0 0 0 A, 


_j -F n -F u ~F 14 -F l5 -F 16 
L M 


-F 21 J F 2 _j -F 2 3 -F^ 

n 

■F 31 -F 32 £ F^J -F M 

7=1 


1 41 


1 51 


1 25 


35 


-f 42 -f 43 £f 4 _, -f* 

7=1 


-F- 


26 


-F 


36 


46 


“^52 “F 33 F» £ ^S-j ~ F 56 


”F<si -F^ 




' bl 


*b2 


'b3 


"b4 




*66 


Equation 15 


The sum of all view factors on one side of the surface equals 1. 


E f w - 1 


Equation 16 
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Substitute into Equation 15 




a i 

0 

0 

0 

0 

0 

<? 2 


0 

A 2 

0 

0 

0 

0 

<?3 


0 

0 

a 3 

0 

0 

0 



0 

0 

0 

A 4 

0 

0 

Q s 


0 

0 

0 

0 


0 

q 6 . 


0 

0 

0 

0 

0 

A e 



Equation 17 


Simplify 

0 * A F E b 


Equation 18 


Solve Equation 11 and Equation 18 for heat load, Q, and blackbody radiation, E b . The most 
direct solution is to set Equation 11 and Equation 18 equal and solve for E b (Equation 19). 

X = [-XRAD+ A f\ T b Equation 19 


This is in the form Ax = b. Solve for E b using Ax = b, and then solve for Q using Equation 18. 


This derivation is valid, but limiting in the analysis. Since Q is not explicitly stated in the 
equation, a single side cannot be constrained to have zero net heat flux (Q t = 0) across it. This 
feature is necessary to decouple one side from the rest of the enclosure. 

Help Desk allows one or more sides of an enclosure to be decoupled from the environmental 
effects. A decoupled side is one with no net heat flux across it. An example is the side of an 
enclosure (or spacecraft bus) that has the thermal radiator attached to it. There is no net heat 
flux from the enclosure through the side with the radiator. The side is effectively decoupled 
from the enclosure and is subject to its own environmental fluxes and effects. Physically this 
is accomplished with insulation between the radiator and enclosure surface. 

A different formulation results when Equation 11 is solved for E b , placed into Equation 18 and 
solved for Q. 

XRAD-' \Q-X\ Equation 20 

[a F XRAD - 1 - [1]) Q = A f XRAD~i X Equation 21 

This is also in the form Ax = b. This form has Q explicitly stated allowing one term to be 
constrained to zero. 

A area matrix (known) 

— F enclosure view factor matrix (known) - - 

XRAD matrix containing the heat radiated from each surface 
X total environmental heat load on each surface 
Q net heat load through a surface 
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To constrain one or more surfaces to have 0 net heat flux across them, set the XRAD term 
equal to 9.9E+9. When the inverse is taken the coefficient of Q will be zero. This also makes 
the right hand side of the equation for Q equal to zero. 

Equation 21 is the basis for both the Enclosure Temperature and Heat Loads analysis and for 
the Orbital Temperatures and Heat Loads analysis. 

This concludes the fundamental mathematical approach to the analyses in the Thermal 
Analyst’s help Desk. The following sections discuss the procedures and equations used to 
derive the terms found in the main equations described above. 
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SURFACE PROPERTIES 


Solar absorptivity and infrared emissivity are two physical parameters that describe the 
surface properties. Solar absorptivity describes the amount of ultraviolet radiation that is 
absorbed into the surface. Infrared emissivitv is the amount of infrared radiation that is 
emitted from a surface. Assuming a gray body, the infrared absorptivity is equal to the 
infrared emissivity. 

These parameters depend on the surface material. For example, white paint has an 
absorptivity of 0.2 and an emissivity of 0.9, while black paint has an absorptivity of 0.9 and 
an emissivity of 0.9. The surface properties also depend on the time in the life of the material 
that is used in the analysis. In the space environment surface degradation occurs and is a 
major factor in the determination of the surface capabilities. Performance parameters are 
generally specified for both beginning of life (BOL) and end of life (EOL). When choosing a 
surface a value is selected from manufacturers specifications, and entered in the appropriate 
box in Help Desk. 

Help Desk allows the user another option. A table containing the surface properties of several 
surface materials is contained in an AI tool called a "Frame". By definition a frame contains 
an "object” plus "slots" for related information. In this manner it is very similar to a small 
database, and it allows reference to a "parent" or "child” frame. A parent or child frame 
contains information about the object: parent information is more general and child 
information is more specific. This referencing is called inheritance. 

When the frame is referenced it compares a user defined text value describing the surface 
application, or the surface material, to either the Surface field or to the Typical.application 
field. Once a match has been found the value for solar absorptance and emittance are 
returned to Help Desk. 

NOTE: The frame call returns the first occurrence of a correct match. If a different value is 
preferred, the user may enter it directly during an analysis. 


Technical Guide 


19 



The frame for the surface properties is displayed below. 


Surface 


'ypicaLappiication 


Solor.obsorptonce 

Emittance 

BOL 

EOL 

BOL 

EOL 

0.90 

0.90 

0.90 

0.90 

0.20 

0.60 

0.90 

0.90 

0.08 

0.21 

0.80 

0.80 

0.08 

0.21 

0.80 

0.80 

0.84 

0.84 

0.85 

0.85 

0.84 

0.84 

0.85 

0.85 

0.35 

0.50 

0.60 

0.60 

0.60 

0.60 

0.60 

0.60 

0.12 

0.18 

0.06 

0.06 

0.12 

0.18 

0.06 

0.06 

0.12 

0.18 

0.06 

0.06 

0.20 

0.60 

0.80 

0.80 

0.70 

0.70 

0.82 

0.82 

0.25 

0.25 

0.05 

0.05 

0.05 

0.05 

0.05 

0.05 




black paint 


white point 


optical solor 
reflector 


opticol solar 
reflector 


graphite/epoxy 


grophite/ epoxy 


aluminized kopton 


tiodized titonium 


oluminum 


aluminum tope 


deposited aluminum 


anodized oluminum 


solar cells 


gold 


multilayered 

insulation 


BOL. Beginning of L 


surfoce data taken fror 


"terior structure 


c.iienna reflector 


"orth panel radiator 


south panei radiator 


solar ponel 


cntenno structure 


•nermal insulation 


coogee motor thermal shield 0.60 


cropellant insulation 


cropellont insulotion 


cropellant insulotion 


nterior structure 


solar panels 


surfoce insulation 


. ; e; EOL, End of Life, 7 years 


r. reference 1 Page 275 


TABLE 1: Surface Properties Frame 
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3.1 Updating Help Desk surface properties option 

The following procedure may be used to add more information to the frame: 

•use an ASCII editor to open the file "SURFACE.TXT" in the directory containing Help Desk 
files. 

•Enter another line of data in any location, making sure to exactly repeat the format of the 
other lines. 

•save and exit 


Rrofessionci 

Choose the surface material: 
[•; block pqint 

I II white .point 71 


jv-|-;graphite/epo)(y 


s.ipio mini zed kopton 
. loiummurn ; •: . 


ilgphloye red-: in suiohds;: 


Figure 5 Surface Material and Application screens 



EXSYS . Prof essionol 

Choose the surface application: 

1 



[ 

[ 

anl <2^00 r-eftec* or- 



. north pane! radiator 

[ 

p s sout ft : panel rod sot or pp- ; : 




[ 

. sctor panel- 

I 

[ 

|, p antenna structure:;.'p^^| 

[ 

.omaf '• m^s 

[ 


[ 

p Wpf opetlo nt : tnsuio ban : ; |§§||| 

f 





To update the user interface for the Surface Properties option, two screens must also be 
modified. The first screen contains the list of surface materials, and the second screen 
contains the list of surface applications (Figure 5). One more box must be added to each 
- screen for each new material and surface application pair in the frame. For more information 
on updating screen files, see Section 2. 
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ORBITAL INFORMATION 


The user is not constrained to perform a thermal analysis about the Earth. Currently the 
equations are valid for any elliptical or parabolic conic section for a "two body problem". A two 
body problem has two assumption: 

1) Each body is spherically symmetric. The mass can be considered as a single point 
mass at the geometric center of the body. 

2) There are no external or internal forces acting on the system other than those defined 
by the two masses. 

The environment of the system is dependent on orbit specific parameters. 

• orbital position 

• spacecraft orientation 

* planet heat flux 

* solar flux 

The user specifies the values for the spacecraft orbit in terms of the six classical orbital 
elements: 

inclination: angle between the momentum vector and the global Z axis. 

nodal crossing Oongitude of ascending node) the angle between the vernal 

equinox direction and the line of nodes, in the equatorial plane. 

argument of the angle between the line of nodes and the direction of periapsis, 

periapsis: in the orbit plane. 

true anomaly at the angle between the direction of periapsis and the position of the 

epoch: spacecraft at a specific time. 

semi-major axis: defines the size of the conic section. 

eccentricity: defines the shape of the conic section (0 for a circle; >0 and <1 for 

an ellipse; 1 for a parabola; >1 for a hyperbola). 

Four parameters associated with the orbit are also used in this analysis: 

Planet Flux Heat load per unit area incident on the surface due to the planet's 

thermal radiation. 

Solar Flux: Heat load per unit area incident on the surface due to the sun's 

radiation. 

Albedo coefficient Percentage of the total solar radiation incident on the planet, that 

is reflected back to space. 

Planet Radius: Radius of the planet under investigation. 


For a given a planet the user may select values for the above variables and enter them into 
the proper boxes. Alternatively, Help Desk allows the user access to its own table called a 
"frame". 
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The frame for the planet information is displayed below. The values in the frame are mean 
values for planet flux, solar flux and albedo. 


PLANET 

PLANET-FLUX 

S0LAR_FLUX 

ALBEDO 

RADIUS 

MERCURY 

4201 

8920 

.058 

2487 

VENUS 

154.2 

2570 

.76 

6187 

EARTH 

236.6 

1352 

.3 

6378 

MOON 

603.25 

1352 

.047 

1738 

MARS 

123 

577.3 

.140 

3380 

JUPITER 

6.1 

49.6 

.51 

71370 

SATURN 

1.8 

14.7 

.5 

60400 

URANUS 

.31 

3.65 

.66 

23530 

NEPTUNE 

.14 

1.48 

.62 

22320 

PLUTO 

.0 

0 

.0 

2284 


TABLE 2: Different Planet Frame 


The following procedure may be followed to add more information to the frame: 

• use an ASCII editor to open the file "PLANET.TXT", which is in the directory 
containing the Help Desk files. 

• Enter another line containing the information for an additional astronomical body (e.g. 
the sun) in any location and duplicate the displayed format. 

• save and exit 

To update the "different planet" option, the screen containing the list of planets must also be 
modified (Figure 5). One more box must be added to each screen to account for the new 
astronomical body. For more information on updating screen files see Section 2. 


24 


Technical Guide 






5 


VIEW FACTORS 


Three sets of view factors are calculated: 

1) enclosure view factors dependant on enclosure geometry 

2) solar view factors dependant on the spacecraft orbit and orientation 

3) planet view factors dependant on the spacecraft orbit and orientation 

5.1 Enclosure View Factors 

Enclosure view factors are the fraction of the total field of view of one side of an enclosure that 
sees each of the other sides of the same enclosure. If the surface under examination is 
concave, it includes the fraction in which it sees itself. 

NOTE: Help Desk currently contains only a rectangular box enclosure, therefore there are no 
concave surfaces. 

The enclosure view factors for Help Desk, Fj t , are calculated using an equation from reference 
6 Appendix C page 825. Equation 22 is a corrected version of the published equation. The last 
two terms in Equation 23 are not raised to any power. 



Figure 7 Two finite rectangles of same length, having one common edge, and at an 
angle of 90* to each other 


kW [ w h JiFZw 2 4 


Equation 22 


where: H = h/1 W = w/1 


(l+W^l+H 2 ) 

' W^UW^H 2 )' 

h 2 {i*w 2 *h 2 ) 

1+W 2 +H 2 

(1 +w\w 2 +h\ 

(i +H z iw i +H 2 ). 


Equation 23 


Technical Guide 


27 


RAG£ SiANK 





Due to the geometry of an enclosure only a few view factors for each surface of an enclosure 
must be solved this way, the remaining can be solved using the following two equations: 


AiFjj= AjFj, 

IFy=l 

Aj is the area of surface i 

Fjj is the view factor from i to j 


5.2 Solar View Factors 

The solar view factor is the fraction of the total radiation emitted by the sun that is directly 
incident on the surface, (e.g. 1: fully facing sun, 0: parallel to the sun’s rays, 0: facing away 
from the sun). 

Sun 

> 

> 

> 

> 

> 

View Factor = 1 View Factor = dos 

Figure 8 Solar View Factor Geometry 



The solar view factors (pj) are found by taking the dot product of the sun vector in local 
coordinates and the unit normal vector of each side. The dot product is: 


S*u-|5||u|cos6 


Equation 24 



u 

Figure 9 Dot Product Geometry 


The solar view factor, p, is defined as cos9, where 9 is the angle between the sun vector and 
the unit normal. 
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Therefore: 





In local coordinates the unit normal vectors are constant. 


Equation 25 



Ui = 

Oi 

+ Oj - lk 

«2 = 

Oi 

- lj + Ok 

U 3 = 

-li 

+ Oj + Ok 

U 4 = 

Oi 

+ lj + Ok 

«S = 

li 

+ Oj + Ok 

U 6 = 

Oi 

+ Oj + lk 


where i, j and k are the direction vectors along the X, Y and Z directions respectively. 
The sun vector in local coordinates is defined as: 

S = X. i+ Y j + Z 8 k 

where i j,k are the unit normals in the coordinate directions, X,Y,Z respectively. 

The magnitude of S is: 

MAGS = SQRTOC , 2 + Y , 2 + Z, 2 ) 
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The solar view factors in local coordinates are: 

Pj = -Z/MAGS 
p 2 = -Y/MAGS 
p 3 = -X/MAGS 
p 4 = Y/MAGS 
p 5 = X/MAGS 
Pg = Z/MAGS 

They are constant in local coordinates since the coordinate axis stays fixed with respect to the 
enclosure. 

5.3 Planet View Factors 

Planet view factors are the fraction of the total thermal radiation emitted by the planet that 
is directly incident on the surface. (e.g. for Low Earth Orbit: -.88 for fully facing planet, .287 
for perpendicular to the planet's surface, 0 for facing away from the planet). It is not one, 
because the field of view out each side, sees space. See Figure 11. 


Planet Surface Planet Surface 



View Factor = .88 View Factor < .88 


Figure 11 Planet View Factor Geometry 


Planet view factors have normal and tangential components. These terms are calculated using 
an equation from reference 12 Appendix C page 829. 

The geometry for the normal component is shown in Figure 12. Equation 26 shows the normal 
component for the planet view factors. 



Figure 12 Plane element dAl to sphere of radius r; normal to center of element passes 
through center of sphere 
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Equation 26 


The geometry for the tangential component is shown in Figure 13. 
tangential component for the planet view factors. 


Equation 28 shows the 



Figure 13 Plane element dAl to sphere of radius r; tangent to element passes through 
center of sphere 



r 


Equation 27 


F„ 

*( n ! j 


Equation 28 
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SPACECRAFT ORIENTATIONS 

The orientation of the spacecraft is defined relative to the sun and the planet about which it 
is orbiting. Two orientations are accommodated in Help Desk (Figure 14): 

• solar inertial orientation 

• planet orientation 

In solar inertial orientation the same side is facing the sun at all times. In planet orientation 
the same side is facing the planet at all times. The orientation is determined by the mission 
of the spacecraft. The effects of the orientation are seen in the solar and planet view factors. 



Solar Inertial Orientation Planet Orientation 

Figure 14 Orientation Geometry 


6.1 Planet orientation: 

Planet orientation has solar view factors that change as the spacecraft travels around the 
planet. The global sun vector rotates with respect to the local spacecraft coordinates and to 
the dot product of the sun vector and the unit normals performed at the orbit position. 

The planet view factors are always normal or tangential to the planet's surface. To account 
for the orbits eccentricity, the planet view factors are calculated at every user defined position 
in the orbit. They are constant only if the orbit is circular. 

6.2 Inertial orientation: 

Solar inertial orientation has constant solar view factors. The global sun vector is rotated to 
local spacecraft coordinates at the line of nodes (a) = 0, and u = 0 in the transformation from 
local to global see Section 7). In this orbit position the dot product of the sun vector and the 
unit normals is performed. This value is kept constant at every point in the orbit. The planet 
view factors change as the spacecraft travels around its orbit. At each location, the normal 
and tangential components are calculated. The planet view factors are then calculated by 
rotating the surfaces through the spacecraft position as follows (see Figure 10 for side number 
referencing): 

FE(1) = FET 

FE(2) = -FEN*SIN© + FET*ABS(COSco) 
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FE(3) = FET*ABS(SIN to) - FEN*COS© 

FE(4) = FEN*SIN© + FET*ABS(COSco) 

FE(5) = FET*ABS(SIN©) + FEN*COS© 

FE(6) = FET 

FET tangential view factor component 
FEN normal view factor component 

© spacecraft position in the orbit. The sum of the argument at periapsis and the true 
anomaly. 

FE(i) view factor for the i lh side 


Technical Guide 



7 


COORDINATE TRANSFORMATIONS 


7.1 Global to Local: 

This coordinate transformation is never needed. When the sun vector is transformed into local 
coordinates, the solution is found solving the equation Ax = b. See Section 5.2 on Solar view 
factors given Global sun vector and Section 6.2 on inertial orientation, for more information. 

7.2 Local to Global: 

When calculating shadowing effects and spacecraft position, the appropriate vectors must be 
rotated from local to global coordinates. The transformation used is from Perifocal coordinates 
to geocentric coordinates. The transformation is taken from reference 2 page 80-83. 


r l 


r P 


Rll R12 R13 

r J 

r t 

* R 

r « 

r w 

where : R = 

Rll R22 R23 
R31 R32 R33 


Where: 
i j and k 
p,q and w 
i 

Q 

(0 


the unit vectors in global coordinates 

the unit vectors in local coordinates. 

inclination 

line of nodes 

line of periapsis 


Equation 29 


Rll = cosQ cosco - sinft since cos i 
R12 = - cosQ sin© - sinQ cos© cos i 
R13 = sinQ sin i 

R21 s sinQ cos© + cosQ sin© cos i 

R22 = -sinQ sin© + cosQ cos© cos i 

R23 = -cosQ sin i 

R31 = sin© sin i 

R32 = cos© sin i 

R33 = cos i 

7.3 Coordinate Transformations Modifications 

To be consistent with the coordinate conventions in IDEAS TM Thermal Model Generator 
(TMG tm ), there is a slight difference between the local to perifocal and global to geocentric 
coordinates. The local and perifocal coordinate systems differ by a 180* rotation about W (see 
Figure 15 and Figure 16). 


Pp 


-h 


G, 


~8, 

P Q 

= 

~ l q 

also 

Gj 

s 

~8j 

Pw 


V 


G* 


8 k 


Equation 30 


Where: 

P p perifocal coordinates P direction 

L p local coordinates P direction. 
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The global coordinates relate to the geocentric coordinate system in the same way. 
Gj the geocentric coordinates in the I direction 
g t the global coordinates in the I direction 



Local Coordinates Perifocal Coordinates 

Figure 15 Local Coordinate Geometry 



VERNAL EQUINOX VERNAL EQUINOX 

DIRECTION DIRECTION 

Global Coordinates Geocentric-equatorial Coordinates 

Figure 16 Global Coordinate Geometry 


Substituting into the transformation relationship. 



Equation 31 
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substituting the negatives through the transformation matrix R giving: 


r i 




Rll R12 -R13 

r J 

^mod 


where R. = 

woa 

R21 R22 -R23 

r k 


v 


-R31 -R32 R33 


Equation 32 


Rmod : 


Rll = 

cos Q. coso) - sinQ sino) cos i 

R12 = 

-cos Q sino> - sinft cosco cos i 

-R13 = 

-sinft sin i 

R21 = 

sinft coso) + cosQ sino) cos i 

R22 = 

-sinQ sino) + cosQ coso) cos i 

-R23 = 

cosQ sin i 

-R31 = 

-sino) sin i 

-R32 = 

-coso) sin i 

R33 = 

cos i 


Equation 32 states the transformation from the perifocal system coordinate system to the local 
coordinate system. One more rotation is needed to transform the system into local spacecraft 
coordinates. The rotation is about the true anomaly at epoch. 


LS is the transformation from local spacecraft to perifocal coordinates, Is is the local spacecraft 
coordinates. 


r P 

r « 

= LS 

K 

where LS = 

cosv -sinv 0 
sinv cosv 0 

Equation 33 

r w 


K 


0 0 1 
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SHADOWING 

Shadow calculations are performed using the Orbital Temperature and Heat Loads option in 
Help Desk. The procedure calculates the spacecraft position, and then determines if the 
spacecraft falls in the shadow of the planet. 

8.1 Position Vector 

The position vector, r, points to the location of the spacecraft in the orbit. It is used to 
determine the spacecraft position relative to the shadow. 

Calculate the position vector in perifocal coordinates, r, and transform into Global coordinates, 
R. r is defined in the perifocal coordinate system as: 

R=r cosv P+r sinvG Equation 34 

P and Q are the directional unit normals in the perifocal coordinate system. 

where the scalar r (the magnitude of the position vector in perifocal coordinates) can be 
calculated with the polar equation of a conic section (Equation 34): 

p is the semi-latus rectum 

e is the eccentricity 

v is the true anomaly at epoch 

a is the semi msgor axis 


Figure 17 Position Vector Geometry 

r= Equation 35 

1+e cosv 

The semi-latus rectum is found using the geometry of the conic section. 

p = a (1-e 2 ) 

8.2 Shadow Calculations 

The geometry for shadowing is shown in Figure 18. For a spacecraft to be in the shadow the 

spacecraft position must satisfy two conditions. 

1) It must be behind the vertical shadow plane defined by a plane that passes through 
the center of the planet and is perpendicular to the suns rays. 

2) It must be inside a cylindrical shadow defined by a cylinder having the same diameter 
as the planet and parallel to the sun's rays. The cylindrical shadow extends away 
from the sun starting at the vertical sun plane where D = 0. 
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D 1 = 0 


CYLINDRICAL 
SHADOW PLANE 
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> 
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* PLANE 

Figure 18 Shadowing Geometry 


In Figure 18, R1 is the position of a spacecraft SI, defined by components D1 and D2. SI is 
on the same side of the vertical plane as the sun placing it outside the shadow. The cylindrical 
shadow plane effect does not need to be calculated since the cylindrical shadow does not exist 
on the sun-side of the vertical shadow plane. R2 is the position of a spacecraft S2, as defined 
with similar components D1 and D2. D1 is the distance from the position to the vertical 
shadow plane and D2 is the distance from the position to the cylindrical plane. 

Define in global coordinates: 

R position vector 

S solar vector 

O origin vector (0,0,0 in our case) 

RE radius of the planet 

V = R-O 


D1 = VS 


D2=<J |K 2 |-|ZW 2 | 

D2<RE A Dl < 0 SHADOWED 
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RADIATE 

An expert system running withing the EXSYS® shell has two main files "expert. system.RUL" and 
"expert_system.TXT". In the case of the Thermal Analyst Help Desk they are "RADI ATE. RUL" 
and "RADIATE.TXT". Both files together define the knowledge base and rule base (see the 
EXSYS® Manual for more details). Help Desk also contains the following supplementary files: 

• RADIATE.SCR contains all the screen descriptions 

• RADIATE.OUT contains the report generator information 

• RADIATE.CFG contains the configuration data for the system such as the location of the 

PASS and RETURN files. 


• RADIATE.COM sets up the command environment 


The rule base contains 47 Rules. All the rules are listed and explained in Appendix A. 

The knowledge base contains 21 qualifiers and 143 variables. All qualifiers are listed in Appendix 
B. 


The screen file, RADIATE.SCR, contains 36 customized screens. All screens and file descriptions 
are displayed in Appendix C. 

The screens use default data files to set the default input values. There are 13 data files. Each 
file is listed and described in Appendix D. 

The report generator file, RADIATE.OUT, has 7 different reports listed. The file uses GOTO 
statements to point to the file description for each report. The report generator file and a sample 
file output for each analysis are listed in Appendix E. 

There are 47 different FORTRAN routines used by Help Desk. They are contained in 20 ASCII 
files, and produce 6 executable files. A list of the files, their routines with the calling arguments, 
and the link commands necessary to produce the executable files are contained in Appendix F. 

The FORTRAN source code is listed in Appendix G. 
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2 SCREENS 

Most of the screens that appear during an analysis run in Help Desk are custom screens, created 
with the EXDESIGN software and an ASCII editor. Only two screens are default EXSYS® screens, 
the description screen, and the HOW? screens. The rest were specifically designed for Help Desk, 
and may be modified for future applications. 

The formats for the screens are in an ASCII file called "RADIATE.SCR". This file is located in the 
same directory as the expert system and data files. The file contains all the screens, and is 
organized by qualifier number or variable name. All screens are displayed in Appendix C. 

The default values displayed on the screens are stored in ASCII data files that are defined for each 
screen in the expert system. Any screen value may be changed by the user. Once a value is 
changed on the screen, it remains at the modified value for the extent of the run. At the start of 
the next run the values will again be the alues specified in the default data files. To modify the 
default values for a set of analysis runs, the information in the data files must be changed. 

There are two types of data files. The first, "filename.DAT" is used during an analysis run and 
contains the data used by the expert system. The second, ”filename.DT!" contains the original 
default values. At the end of an analysis sesion, the values may be reset to the original default 
values by copying the default file "filename.DT!" to the working file "filename.DAT". 

The format for the input data files is one variable per line. Each line corresponds to a variable 
defined in the screen file, in order of appearance. When EXSYS® displays a screen, it places the 
information from the data file into the available input boxes. All of the applicable data files are 
shown in Appendix D. 

A screen is attached to a qualifier or a variable. When the expert system prompts the user for a 
qualifier or variable, it searches the screen file for a screen attached to the qualifier or variable. 
EXSYS® displays the screen defined immediately following the variable name. There are 
approximately 30 different screens defined in the screen file for this expert system. 

The following is a sample screen with the screen description listed and each command explained. 
While the user creates the screen shown, EXDESIGN creates the file described. 


EXSYS Professional 
Planet Constants 


Albedo coefficient: 
Solar Flux: W/m n 2 
Planet Flux: W/m A 2 
Planet radius: km 
II different planet 



OK 


Figure 19 Sample Data screen 
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the screen is ottoched to Qualifier 17 (~Q 1 7) 


~Q17 

~SET_WIN 88. 0. 359, 225 the window is located at 88,0 and has a size of 359 wide 

and 225 toll. 


"LINE-COLOR DKGRAY 
"FOREGROUND DKGRAY 
"BACKGROUND CYAN 
~CLS 

~CURSET 102, 28 
Planet Constants 
"CURSET 10. 60 
Albedo coefficient: 

"CURSET 10. 100 

Solor Flux: W/m-2 

"CURSET 10. 140 

Planet Flux: W/m-'2 

"EDIT "[ALB]". 225, 45. 103. 35 -M:1 

"EDIT "[S]". 225. 85, 103. 35 -M:1 

-EDIT "[E]". 225, 125. 103. 35 -M:1 

"EDIT "[RE]". 225. 165, 103. 35 -M:l 


The line color is dark groy 
foreground color is dark gray 
the background is cyan. 

In order to make the entire screen background one color 
instead of just the background of the text boxes, the "CIS 
commond must follow the BACKGROUND command. 

The first cursor position is at 102. 28 

the text "Planet Constants" 

the cursor is placed at 10,60 

the text "Albedo coefficient:" 

the cursor is placed at 10,100 

the text "Sotor Flux: W/m~2" 

the cursor is placed at 10,140 

the text "Planet Flux: W/m~2" 

the edit box is ottoched to the variable [ALB] (albedo). Its 
is located at 225, 45. Its width is 103 and height 35. The 
-M:1 option states that the box is only to allow 1 line of 
entry. 

the edit box is attached to the variable [S] (solar flux). 
Its is located at 225, 48. Its width is 103 and height 35. 
The -M:1 option states that the box is only to allow 1 line 
of entry. 

the edit box is ottoched to the variable [E] (planet flux). 
Its is locoted at 225, 125. Its width is 103 and height 35. 
The -M:1 option states that the box is only to allow 1 line 
of entry. 

the edit box is attached to the variable [RE] (plonet 
radius). Its is located at 225, 165. Its width is 103 and 
height 35. The -M:1 option stotes that the box is only to 
allow 1 line of entry. 


44 


Programmer's Guide 



^BUTTON "different olanet", "017 1", 23. 205. 15S. 35 -I 


the button box is labelled "different planet". If the button 
is pressed Ouolifier 17 (017) is assigned the value 1 (a 
different planet is requested). The box is located at the 
position 23,205, it is 158 wide and 35 high. The -I option 
says to return the value immediately after the button is 
pressed. 

~8UTT0N "OK", "017 2", 229, 205, 75, 35 -I the button box has the word OK on it. It returns the 

volue 2 to 017. It is located at 229, 205, has a width of 
75 and a height of 35 and returns immediately if pressed. 

~E01T_FILL ALBSCR.DAT the ~E01T_FILL command specifies the input data file to be 

accessed to fill in the edit boxes. 

~END defines the end of the screen 


2.1 Default Data Files 

The Planet Constants screen shown in Figure 19 is an example of a data screen using default data 
values. This data screen is used to enter the variables for albedo coefficient, solar flux, planet flux 
and planet radius. The following is the associated default data file: ALBSCR.DAT 

.3 

1356 

236 

6378 

The value .3 is placed in the edit box for [ALB] (albedo), 1356 is for [S] (solar flux), 236 is for [E] 
(planet flux), and 6378 for [RE], the planet radius. 
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3 RUNNING EXTERNAL PROGRAMS 


Throughout Help Desk, external programs are used to perform complex calculations. They are 
used to minimize the complexity of the expert system and to increase the capability of Help Desk. 
When executing Help Desk a user can determine when the expert system is running an external 
FORTRAN code, because the screen goes blank momentarily. This blanking is due to the 
Windows-DOS interface. 


The executable Fortran files are called using external DOS batch files. The batch files used are: 


FETANOR.BAT 

PLANET.BAT 

ORBIT.BAT 

QENCL.BAT 

SOLAR.BAT 

SUNVECT.BAT 


Calculates the normal and tangential components for the planet view factors 

Calculates the planet view factors for each side of an enclosure 

Performs the orbital analysis 

Performs the enclosure analysis 

Calculates the solar view factors 

Rotates the global sun vector into local coordinates and calculates the solar 
view factors 


A typical DOS batch file, PLANET.BAT, is shown below. 

TYPE NUL > EXSYS.EP.END 
EARTHV.EXE 

RENAME EXSYS.EP.END EXSYS.EP.FIN 

The first line creates a null file called EXSYS.EP.END. 

The second line runs the executable EARTHV.EXE. 

The final line renames the file to EXSYS.EP.FIN. 

This batch file procedure is required to run non-windows applications. EXSYS® realizes that the 
analysis is finished when it encounters the file EXSYS.EP.FIN. As soon as this file exists in the 
working directory. Help Desk continues it’s analysis. If the file exists prior to completion of the 
executable file. Help Desk will take control of the analysis. This may place Help Desk into an 
infinite loop. The file, EXSYS.EP.FIN, is created in the method shown above to ensure that the 
file exists only when the external analysis is complete. If the file is created during the executable, 
it may be created too early. 

The variables needed to execute the code are passed to a file called RADIATE.DAT. The file is 
defined in the configuration file for the expert system. The path to the file must be defined as 
the path to the directory containing the batch files and executable files. 

An example of a rule call statement to the batch file PLANET.BAT is shown below. 

IF: 

The spacecraft orientation is {inertial! 
and: The spacecraft orbit is {LEO (Low Earth Orbit)} 

THEN: 

X> RUN(PLANET.BAT [INC], [ECC],[WOMEGA] ,[SMA],[V0],[RE] M /Z) 
and: Planet view factors are (calculated) 

This rule states that if the spacecraft has an inertial orientation, and it is in a low earth orbit, 
then the system will run the file PLANET.BAT. Six variables are passed to the file 
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"RADIATE.DAT". Two parameters are also specified, /M and /Z. /M states that multiple variables 
will be returned in the file "RETURN.DAT” and /Z states that a Ctrl-Z will be added to the end 
of the file "RADIATE.DAT". Notice that /B option, which the EXSYS® manual suggests using, is 
not used in this case. 

NOTE: The /B option is not used for this type of batch file. The /B option indicates that a batch 
file, containing only DOS commands , is to be run. In this example the batch file includes 
a FORTRAN subroutine, if the /B option were used in this case, data would not be passed 
to the RADIATE.DAT file. 
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4 ADDING A NEW FEATURE 


The following procedure was used during the development of Help Desk and is presented here to 
illustrate the way in which new features can be added. 

The only planet data included in the early versions were the default values for Earth. For other 
planets the user was required to enter the necessary values for the planetary parameters. The 
expert system has been modified to provide the values for the other planets and the Moon while 
retaining the option of entering user selected values. 

User provided values were entered into the two screens shown in Figure 20. One was used in 
surface analyses and the other in enclosure analyses. No logic was present to direct the expert 
system to ask for either a planet selection or the solar and planet flux values. The values were 
entered only as additional information on screens, which were activated by logic that was 
searching for emissivity or albedo, respectively. 


X s I S . P f of e S 3 KOT Of 
Analysis Constants 
Surface emissivity: 
Surface absorptivity: 
Solar Flux: W/m~2 
Planet Flux: W/m~2 
Bolar view factor: 

Planet view factor: 
Albedo coefficient: 



i£|sts|||pce : ! 



: .EXSX$fProf essi o no 11 


1356 


Analysis Constants 
Solar Flux: W/m~2 

Planet Flux: W/m-2 
Albedo coefficient: 


236 


.3 


Planet radius: km 


6378 


plonet 


liilil 


Figure 20 Original Constants screen 


To add the other planets it was necessary to split up the information on the screens and add a 
button called "different planet". This button is similar to the "assistance” button in appearance 
and function. 

The new analysis contants screens called "Surface Constants" and "Planet Constants" are shown 
in Figure 22. The two screens are now used in both the surface analyses and the enclosure 
analyses. A new frame containing the values for all the planets was also created. The following 
sequence describes the process used to modify the expert system. 

The first step in adding a new feature is to determine the new qualifiers and variables required 
for the modification. In this case it requires only one new qualifier, one new rule and a new frame. 
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The new qualifier is: 

Qualifier A different planet is 

V requested 

V not requested 

and the new rule 

Rule: If a different planet is requested 

Then: [E] IS GIVEN THE VALUE .... 
and [S] IS GIVEN THE VALUE... 
and [ALB] IS GIVEN THE VALUE... 

In this example it is also necessary to define a new frame called "PLANET.TXT", shown below. 
The frame contains each of the planets, the average values of solar flux and planet flux, and the 
radius of the planet. Different values may be entered directly at the Planet Constants screen. The 
frame may also be editted to modify the values for all future analysis runs. 


PLANET 

PLANET.FLUX 

SOLAR.FLUX 

ALBEDO 

RADIUS 

MERCURY 

4201 

8920 

.058 

2487 

VENUS 

154.2 

2570 

.76 

6187 

EARTH 

236.6 

1352 

.3 

6378 

MOON 

603.25 

1352 

.047 

1738 

MARS 

123 

577.3 

.140 

3380 

JUPITER 

6.1 

49.6 

.51 

71370 

SATURN 

1.8 

14.7 

.5 

60400 

URANUS 

.31 

3.65 

.66 

23530 

NEPTUNE 

.14 

1.48 

.62 

22320 

PLUTO 

.0 

0 

.0 

2284 


The variables used to retrieve the information from the frame must be defined. The format of a 
frame statement is: 


FRAME(”filename’’,test_expression,slot_reference) 


In this example the frame statements are: 

FRAMECPLANET.TXT”, $PLANET$==[PLANET],#SOLAR_FLUX#) 
FRAMECPLANET.TXT", $PLANET$==[PLANET],#PLANET_FLUX#) 
FRAME(”PLANET.TXT",$PLANET$==[PLANET],#ALBEDO#) 
FRAMECPLANET.TXT” ,$PLANET$==[PLANET],#RADIUS#) 
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Using the frame statements, the entire new rule is: 

Rule: If a different planet is requested 

Then: [S] IS GIVEN THE VALUE 

FRAMECPLANET.TXT” ,$PLANET$==[PLANET],#SOLAR_FLUX#) 
and [E] IS GIVEN THE VALUE 

FRAME("PLANET.TXT",$PLANET$==[PLANET],#PLANET_FLUX#) 
and [ALB] IS GIVEN THE VALUE... 

FRAME("PLANET.TXT”,$PLANET$==[PLANET],#ALBEDO#) 

The variable [PLANET] contains the planet name in a string format. To allow the user to choose 
a planet name, a screen is attached to [PLANET]. When the inference engine processes the frame 
statements above, it tries to derive the variable [PLANET], If it cannot derive it, it prompts the 
user with the screen that is attached to the variable name. The planet screen is shown in 
Figure 21. 


EXSYS Professional.::::-: 

Choose the Planet: 


-Cf iMEPCOPY : : 






i • MOON . ' 

i 

MARS . ; 


'If. : SAT0PN||:;||| 

;J. URANUS 


PLSJ 


Figure 21 Planet Selection Screen 


The new screen definition is shown below. The format is created by the following procedure: 

• Define the window size, line color, background and foreground colors. 

• Place the curser at 41, 36 and enter the text "Choose the planet:" 

• Place a series of nine buttons each with 

NAME 

VARIABLE and its returned value if the button is pressed 
POSITION OF THE BOX: 

X and Y (component of the upper left hand comer of the box), WIDTH, HEIGHT 
-I return immediately after the button is pressed. 

• END 
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-[PLANET] 

-SET.WIN 88. 0. 255. 460 
-LINE-COLOR DKGRAY 
-EOREGROUND OKGRAY 
-BACKGROUND CYAN 
-CLS 

-CURSET 41. 36 


Choose the plonet: 

-BUTTON "MERCURY". "[PLANE 
-BUTTON 'VENUS". "[PLANET] 
-BUTTON "EARTH". "[PLANET] ' 
-8UTT0N "MOON ", "[PLANET] 
-BUTTON "MARS ". "[PLANET] 
-BUTTON "JUPITER". "[PLANET 
-BUTTON "SATURN". "[PLANET' 
-BUTTON "URANUS". "[PLANET 
-BUTTON "NEPTUNE". "[PLANE 
-BUTTON "PLUTO". "[PLANET] ' 
-END 


T] 'MERCURY'". 25, 50, 192. 33 -I 
'VENUS'". 25. 90. 192. 33 -I 
EARTH’", 25, 130. 192. 33 -I 
'MOON'", 25, 170. 192, 33 -I 
MARS’". 25, 210, 192, 33 -I 
'JUPITER'", 25, 250. 192, 33 -I 
'SATURN'". 25. 290. 192, 33 -I 
’URANUS’". 25. 330, 192, 33 -I 
] 'NEPTUNE'". 25. 370. 192, 33 -I 
PLUTO’". 25, 410, 192, 33 -I 


The original screens were then modified to allow the use of the planet button. The original 
analysis screens were split into two seperate screens with no repitition of user inputs. The first 
contains all of the surface parameters: emissivity, absorptivity, solar and planet view factors. The 
second contains all of the parameters related to the planet: albedo coefficient, solar and planet 
fluxes, and planet radius. Each of the two screen contain a button that allows the user to access 
information from a frame, as well as a button accepting the user input values. 


EXSYS Professional 


Surface Constants 
Solar view factor 
Planet view factor 

Surface emissivity: 

Surface absorptivity: 


248 


.05 


.05 


assistance 


EXSYS Professional 
Planet Constants 


Albedo coefficient: 
Solar Flux: W/m A 2 
Planet Flux: W/m A 2 
Planet radius: km 
|| different planet 



Figure 22 New Analysis Constants screens 


Both of the screens are displayed by the inference engine prompting the user for values to 
qualifiers. The "assistnace” and "different planet" buttons are attached to qualifiers. These button 
values drive the logic flow to the screens. 
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SURFACE CONSTANTS SCREEN 

Ouolifier 16 assistance is: 

V needed 

V not needed 


~Q16 

~SET_WIN 88, 0, 380, 280 
-LINE-COLOR BLACK 
-FOREGROUND BLACK 
"BACKGROUND CYAN 
-CLS 

-CURSET 110.35 
Surfoce Constants 
"CURSET 30, 75 
Solar view factor: 
"CURSET 30, 115 
Planet view factor: 


"CURSET 30, 155 
Emissivity: 
-CURSET 30, 195 


Absorptivity: 

"EDIT "[MU]". 225. 50. 103, 35 -M:1 
"EDIT "[FE]", 225, 90, 103, 35 -M:l 
"EDIT "[EPSILON]". 225. 130, 103, 35 -M:1 
~EDIT "[ALPHA]". 225. 170, 103, 35 -M:l 
"BUTTON "OK". "016 2". 261, 228, 75. 35 -I 
"BUTTON "assistance". "q16 1", 38. 228, 158, 
"EDIT JILL ALPSCR.DAT 
"END 
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PLANET CONSTANTS SCREEN 

Qualifier 17 A different plonet is 

V requested 

V not requested 

"Q17 

"SET-WIN 88, 0, 360, 271 
"LINE-COLOR BLACK 
"FILL-COLOR BLACK 
"FOREGROUND BLACK 
"BACKGROUND CYAN 
~CLS 

"CURSET 102, 30 
Plonet Constants 
"CURSET 10, 70 
Albedo coefficient: 

"CURSET 10. 110 
Solor Flux: W/m-2 
"CURSET 10, 150 
Plonet Flux: W/m-2 
"CURSET 10, 190 
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Planet radius: km 

-EDIT "[ALB]", 225, 45. 103. 35 -M:1 

-EDIT ”[S M . 225. 85, 103, 35 -M:1 

-EDIT M [E]". 225, 125, 103,35 -M:1 

-EDIT "[RE]". 225, 165, 105. 35 -M:1 

-BUH0N "OK", "017 2". 241. 217. 75, 35 -I 

-BUTTON "different planet", “q17 1", 23, 216, 158, 35 -I 

-EDIT-FILL ALBSCR.DAT 

-END 

Once the above modifications are complete, change the input data files to delete the variable 
moved. 

Finally, run the system through every option to make sure that there are no new conflicts. 
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FUTURE DEVELOPMENT 


The Thermal Analyst's Help desk can be expanded to incorporate a number of additional 
features. Some of the features under consideration are listed here. 

• include double sided radiator 

• add the planet vector calculations 

• add the beta angle calculations 

• simplify the parametric trade-off capabilities 

• add the planet shield calculations to the Orbital Thermal Analysis 

• add the ability to modify the solar view factors for the Enclosure Analysis and allow solar 
shields to the Orbital Thermal Analysis 

• add shadow calculations to the Enclosure analysis 

• incorporate the Thermal Electric Cooler expert system into Help Desk 

• expand the Surface Properties Analysis 
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Purpose: 

Finish the thermal analysis and write the report file. The rule keeps 
thermal analysis separate from other analyses that may be added, (e.g. 
Thermal Electric Coolers). 

Premise: 

The qualifier "Thermal analysis is" (Q20) has the value "completed". 
The value for Q20 is derived using Rules 2-8. 

Conclusion: 

Execute the report generator file "RADIATE.OUT", see the report 
generator section of this document for further explanation. 

Define the Choice "Thermal analysis is complete" to have a confidence 
value of 1. 


Rules 


/* RULE NUMBER: 1 
RULE: THERMAL COMPLETE 

IF: 

Thermal analysis is (completed) 

THEN: 

X> REPORT(RADIATE.OUT) 
and: > Thermal analysis is complete - Confidence=l 
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Rules 


A4 


Purpose: 

Define the area of a surface. When the inference engine is proving a 
premise, it finds the first rule, in sequential order, that can define the 
variable or qualifier under investigation. In the example, the inference 
engine is trying to define the qualifier "Thermal analysis is", as show 
in the premise of Rule 1. 

Premise: 

The qualifier "Value to calculate is" (Ql) has the value "radiator area 
(given temperature)". The user is prompted with the screen containing 
the main menu. When the user selects an analysis, the value is 
returned and assigned to Ql. 

The variable [AREA], the area of the surface, must be within .1% of 
the equation shown. The value for the variable [AREA] is derived 
using Rule 9. The terms on the right hand side of the equation are 
explicitly stated to require each term in the equation to be evaluated 
regardless of the initialized values. See Rule 9 for more information 
on the derivations. The characters ~= stand for "approximately equals", 
which is defined to be ±X% of the value. 

The variables [QRAD] the heat load radiated from the surface, [QSUN] 
the heat load on the surface due to the sun, [QALBEDO] the heat load 
on the surface due to the planet’s albedo, [QPLANET] the heat load on 
the surface due to the planet, must sum to a value greater than zero. 
Each of the variables are initialized to zero, thus the sum equals zero. 
If a value cannot be derived, the expert system will use the initialized 
value prior to prompting the user for the information. Since the values 
for each variable is displayed on the final report page, the expert 
system needs a value for them. If they are not initialized to zero, the 
values would be derived, regardless of the user selection on the 
radiation effects screen. Thus, the expert system is required to derive 
all the terms from existing rules before accepting the initialized value. 

Conclusion: 

Display the report page [REPORT1]. See Appendix C for the figure of 
[REPORT1]. 

The qualifier "Thermal analysis is" is assigned the value "completed". 


f* RULE NUMBER: 2 
RULE: AREA FOUND 

IF: 

Value to calculate is {radiator area (given temperature)! 
and: [AREA] ~= 

-([QP0WER]+[QENCL])/([XSUN]+[XPLANET]+[XALBED01-[XRAD]) 
and: [QRAD] + [QPLANET] + [QSUN] + [QALBEDO] > 0 

THEN: 

[REPORT] IS GIVEN THE VALUE [REPORT1] 
and: Thermal analysis is [completed) 



Purpose: 

Defines an area has been calculated less than zero. The rule allows 
the completion of the analysis if Rule 2 defines a negative area. 

Premise: 

The qualifier "Value to calculate is" (Ql) has the value "radiator area 
(given temperature)”. The user is prompted with the screen containing 
the main menu. When the user selects an analysis, the value is 
returned and assigned to Ql. 

The variable [AREA), the area of the surface, has been calculated to 
have a value less than zero. 

Conclusion: 

The user is informed of a bad area using [REPORT2]. See Appendix 
C for the figure of [REPORT2]. 

The qualifier "Thermal analysis is" is assigned the value "completed". 


Rules 


/* RULE NUMBER: 3 
RULE: BAD AREA 

IF: 

Value to calculate is (radiator area (given temperature)! 
and: [AREA] < 0 

THEN: 

[REPORT] IS GIVEN THE VALUE [REP0RT21 
and: Thermal analysis is (completed) 
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Purpose: 

Define the temperature of a surface. 

Premise: 

The qualifier "Value to calculate is" (Ql) has the value "surface 
temperature", which is derived from the main menu with the analysis 
chosen by the user. 

The qualifier "Surface temperature to be calculated is” (Q12) has the 
value "radiator temperature given area". The user is prompted with 
a screen containing two choices, radiator temperature and 
environmental temperature. Once the user selects a value, it is 
assigned to Q12. 

The variable [TEMP], the temperature of the surface, must be 
calculated and greater than zero. The value for [TEMP] is derived 
using Rule 10. 

The variables [QRAD] the heat load radiated from the surface, [QSUN] 
the heat load on the surface due to the sun, [QALBEDO] the heat load 
on the surface due to the planet’s albedo, [QPLANET] the heat load on 
the surface due to the planet, must sum to a value greater than zero. 
Each of the variables are initialized to zero, thus the sum equals zero. 
If a value cannot be derived, the expert system will use the initialized 
value prior to prompting the user for the information. Since the values 
for each variable is displayed on the final report page, the expert 
system needs a value for them. If they are not initialized to zero, the 
values would be derived, regardless of the user selection on the 
radiation effects screen. Thus, the expert system is required to derive 
all the terms from existing rules before accepting the initialized value. 

Conclusion: 

Display the report page containing the final information [REPORT1]. 
See Appendix C for the figure of [REPORT1]. 

The qualifier "Thermal analysis is" is assigned the value "completed". 


/* RULE NUMBER: 4 
RULE: SUR TEMP FOU 

IF: 

Value to calculate is [surface temperature] 
and: Surface temperature to be calculated is (radiator temperature given 
area) 

and: [TEMP] > 0 

and: [QRAD] + [QSUN] + [QPLANET] + [QALBEDO] + [QENCL] > 0 
THEN: 

I REPORT | IS GIVEN THE VALUE IREPORT1I 
and: Thermal analysis is (completed) 



Purpose: 

Define the environmental temperature of a surface. 

Premise: 

The qualifier "Value to calculate is" (Ql) has the value "surface 
temperature", which is derived from the main menu with the analysis 
chosen by the user. 

The qualifier "Surface temperature to be calculated is” (Q12) has the 
value "environmental temperature". The user is prompted with a 
screen containing two choices, radiator temperature and environmental 
temperature. Once the user selects a value, it is assigned to Q12. 

The variable [TEMP], the temperature of the surface, must be 
calculated and greater than zero. The value for [TEMP] is derived 
using Rule 10. 

Conclusion: 

Display the report page containing the final information [REPORT3]. 
See Appendix C for the figure of [REPORT3]. 

The qualifier "Thermal analysis is" is assigned the value "completed". 
Purpose: 

Define the enclosure surface temperatures and heat loads. 

Premise: 

The qualifier 'Value to calculate is" (Ql) has the value "enclosure 
temperatures and heat loads", which is derived from the main menu 
with the analysis chosen by the user. 

The qualifier "The temperature values are" (Q7) has the value 
"calculated". Help Desk derives Q7 using Rule 24. 

Conclusion: 

Display the report page containing the final information [REPORT4]. 
See Appendix C for the figure of [REPORT4]. 

The qualifier "Thermal analysis is" is assigned the value "completed”. 


Rules 


/* RULE NUMBER: 5 
RULE: ENV TEMP FOU 

IF: 

Value to calculate is [surface temperature) 
and: Surface temperature to be calculated is [environmental temperature) 
and: [TEMP] > 0 

THEN: 

[REPORT] IS GIVEN THE VALUE [REPORT3] 
and: Thermal analysis is (completed) 


/* RULE NUMBER: 6 
RULE: ENCL FOUND 

IF: 

Value to calculate is (enclosure temperatures and heat loads) 
and: The temperature values are (calculated) 

THEN: 

[REPORT] IS GIVEN THE VALUE [REPORT4] 
and: Thermal analysis is (completed) 
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Purpose: 

Define the enclosure surface temperature and heat loads over an entire 
orbit. 

Premise: 

The qualifier "Value to calculate is" (Ql) has the value "orbital average 
temperatures and heat loads", which is derived from the main menu 
with the analysis chosen by the user. 

The qualifier "Orbital analysis is" (Qll)has the value "finished". Help 
Desk derives the value for Qll using Rules 25-26. 

Conclusion: 

Display the report page containing the final information [REPORTS]. 
See Appendix C for the figure of [REPORT5]. 

The qualifier "Thermal analysis is" is assigned the value "completed". 


Purpose: 

Calculate the surface emissivity and absorptivity. 

Premise: 

The qualifier 'Value to calculate is" (Ql) has the value "surface 
properties", which is derived from the main menu with the analysis 
chosen by the user. 

The qualifier "Surface properties are" (Q13) has the value "calculated". 
Help Desk derives the value for Q13 using Rules 38-41 or Rules 43-46. 

Conclusion: 

Display the report page containing the final information [REPORT6]. 
See Appendix C for the figure of [REPORT6]. 

The qualifier "Thermal analysis is" is assigned the value "completed". 


A8 

/* RULE NUMBER: 7 
RULE: ORBIT FOUND 

IF: 

Value to calculate is (orbital average temperatures and heat loads) 
and: Orbital analysis is (finished) 

THEN: 

[REPORT] IS GIVEN THE VALUE [REPORTS] 
and: Thermal analysis is (completed) 


/* RULE NUMBER: 8 
RULE: PROP FOUND 

IF: 

Value to calculate is (surface properties) 
and: Surface properties are (finished) 

THEN: 

[REPORT] IS GIVEN THE VALUE [REPORT6] 
and: Thermal analysis is (completed) 



Purpose: 

Calculate the area of the surface (see Rule 2). 

A « P power* Q enclosure) 

^SUN*^ PLANET* ^ALJSEDO~^ RAD 

Premise: 

The qualifier ,r Value to calculate is" (Ql) has the value "radiator area 
(given temperature)", which is derived from the main menu with the 
analysis chosen by the user. 

The sum of the variables [XSUN] the heat flux on the surface due to 
the sun’s radiation, [XPLANET] the heat flux on the surface due to the 
planet's radiation, [XALBEDO] the heat flux on the surface due to the 
albedo of the planet and -[XRAD] the heat flux radiated from the 
surface is greater than zero. This requires all of the heat flux terms 
to be evaluated and that the sum is not zero, preventing division by 
zero in the THEN clause. [XSUN] is derived using Rule 19. 
[XPLANET] is derived by Rule 20. [XALBEDO] is derived by Rule 21. 
[XRAD] is derived using Rule 23. 

Conclusion: 

The variable [AREA] is defined. [QPOWER], the internal power to be 
dissipated, is input by the user. [QENCL] the heat load from the 
enclosure is derived using Rule 17. The remaining terms are derived 
in the premise. 


/* RULE NUMBER: 9 
RULE: AREA 

IF: 

Value to calculate is (radiator area (given temperature)) 
and: [XSUN] + [XPLANET] + [XALBEDO] - [XRAD] <> 0 

THEN: 

[AREA] IS GIVEN THE VALUE 

-([QPOWER]+fQENCLl)/((XSUN1+|XPLANETl+[XALBEDO]-[XRADl) 


Rules 
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Rules 


A10 


Purpose: 

Calculate the temperature of the surface (see Rule 4 and Rule 5) 

E b = oT* =» r=^j 4 

Premise: 

The qualifier "Value to calculate is" (Ql) has the value "surface 
temperature", which is derived from the main menu with the analysis 
chosen by the user. 

The variable [EB], the blackbody radiation on the surface, must be 
greater than zero. It is derived using Rule 11. 

The variable [SIGMA], the Stephen-Boltzmann constant must be 
greater than zero. It is derived using Rule 22. EXSYS® does not use 
scientific notation and cannot initialize a variable as small as 
[SIGMA]. However, it can be assigned a small value in a rule. 

Conclusion: 

The variable [TEMP] is assigned the value as shown. 


/* RULE NUMBER: 10 
RULE: TEMP 

IF: 

Value to calculate is (surface temperature) 
and: [EB] > 0 
and: [SIGMA] > 0.0 

THEN: 

[TEMP] IS GIVEN THE VALUE SQRT(SQRT([EB]/[SIGMA])) 



Purpose: 

Calculate the blackbody radiation on the surface (see Rule 10). 

17 _ if v j_y _ POWER* Q ENCLOSURE} ] 

C t> “ e \ a sun a planet a albedo" ^ I 

Premise: 

The qualifier "Surface temperature to be calculate is" (Q12) has the 
value "radiator temperature given area". The user is prompted with 
a screen containing two choices: radiator temperature and 
environmental temperature. The value selected is assigned to Q12 (see 
Rule 4 and Rule 5). 

The variable [AREA] is greater than zero. The inference engine cannot 
derive a value for the area, so it prompts the user for a value with a 
screen that attached to the variable [AREA], 

The variable [EPSILON], the emissivity of the surface, is greater than 
zero, which is derived using Rules 43-46. Even though Rules 38-41 can 
also derive values for [EPSILON], the rules were marked False as soon 
as Q1 was assigned a value from the main menu, so the inference 
engine will never try to "fire" them. The user will be allowed to enter 
the known value or ask Help Desk for "assistance". 

Conclusion: 

The variable [EB] is defined. [XSUN] the heat flux on the surface due 
to the sun is calculated using Rule 19. [XPLANET] the heat flux on 
the surface from the radiation of the planet is derived by Rule 20. 
[XALBEDO] the heat flux on the surface from the albedo of the planet 
is derived by Rule 21. [QPOWER] the internal power to be dissipated 
is input by the user. [QENCL] the heat load from the enclosure is 
derived using Rule 17. The remaining terms are derived in the 
premise. 


/* RULE NUMBER: 11 
RULE: BLACKBODY 

IF: 

Surface temperature to be calculated is Iradiator temperature given 
area) 

and: [AREA] > 0 
and: [EPSILON] > 0 

THEN: 

[EB] IS GIVEN THE VALUE 

(fXSUN]+[XPLANETl+fXALBED01+(([QP0WER]+fQENCLl) 

/[ AREA )))/! EPSILON I 
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Purpose: 

Calculate the blackbody radiation on the surface (see Rule 10). 


/* RULE NUMBER: 12 
RULE: ENVIRON TEMP 


^b = g (^SUN + ^PIANET + ^ ALBEDO) 


Premise: 

The qualifier "Surface temperature to be calculate is" (Q12) has the 
value "environmental temperature". The user is prompted with a 
screen containing two choices, radiator temperature and environmental 
temperature. The value selected is assigned to Q12 (see Rule 4 and 
Rule 5). 

The variable [EPSILON], the emissivity of the surface, is greater than 
zero, which is derived using Rules 43-46. Even though Rules 38-41 can 
also derive values for [EPSILON], the rules were marked False as soon 
as Q1 was assigned a value from the main menu, so the inference 
engine will never try to "fire" them. The user will be allowed to enter 
the known value or ask Help Desk for "assistance". 

Conclusion: 

The variable (EB1 is defined. TXSUNl the heat flux on the surface due 
to the sun is calculated using Rule 19. [XPLANET] the heat flux on 
the surface from the radiation of the planet is derived by Rule 20. 
[XALBEDO] the heat flux on the surface from the albedo of the planet 
is derived by Rule 21. [EPSILON] is derived in the premise. 


IF: 

Surface temperature to be calculated is (environmental temperature] 
and: [EPSILON]>0 

THEN: 

[EB1 IS GIVEN THE VALUE 
([XSUN]+[XPLANET]+[XALBEDO])/[EPSILON] 


/* RULE NUMBER: 13 
RULE: SUN HEAT LOAD 


Purpose: 

Calculate the heat load on the surface due to the sun’s radiation. 

Qsun ~ ^SUN^ 


Premise: 

The qualifier "radiation effects to include" (Q2) has the value "sun”. 
The user is prompted with a screen containing all the possible 
radiation effects on the surface. The user selects one or multiple 
values, which are assigned to Q12. 

Conclusion: 

The variable [QSUN] is assigned a value as shown. [AREA] is derived 
using Rule 9. [XSUN] is derived using Rule 19. 


IF: 

Radiation effects to include (sun) 

THEN: 

[QSUN] IS GIVEN THE VALUE [XSUN]*[AREA] 



Purpose: 

Calculate the heat load on the surface due to the radiation of the 
planet. 

Q PLANET = ^ PLANET^ 

Premise: 

The qualifier "radiation effects to include” (Q2) has the value "planet". 
The user is prompted with a screen containing all the possible 
radiation effects on the surface. The user selects one or multiple 
values, which are assigned to Q12. 

Conclusion: 

The variable fQPLANETl is assigned a value as shown. (AREA) is 
derived using Rule 9. [XPLANET] is derived using Rule 20. 

Purpose: 

Calculate the heat load on the surface due to the albedo of the planet. 

Q ALBEDO = ^ ALBEDO^ 

Premise: 

The qualifier "radiation effects to include" (Q2) has the value "albedo". 
The user is prompted with a screen containing all the possible 
radiation effects on the surface. The user selects one or multiple 
values, which are assigned to Q12. 

Conclusion: 

The variable [QALBEDO] is assigned a value as shown. [AREA] is 
derived using Rule 9. [XALBEDO] is derived using Rule 21. 


Rules 


t* RULE NUMBER: 14 
RULE: PLANET HEAT LOAD 


IF: 

Radiation effects to include (planet) 

THEN: 

[QPLANET] IS GIVEN THE VALUE (XPLANET]*[AREA) 


/* RULE NUMBER: 15 
RULE: ALBEDO HEAT LOAD 

IF: 

Radiation effects to include (albedo) 

THEN: 

[QALBEDOl IS GIVEN THE VALUE (XALBED01*(AREA1 
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Rules 


Purpose: 

Calculate the heat load radiated from the surface. 

Qrad ~ ^racA 

Premise: 

The value of [XRAD] must be greater than zero. This forces the 
evaluation of the variable TXRAD1. This state the radiative capability 
of the surface. The value is calculated regardless of which radiation 
effects have been included. 

Conclusion: 

The variable [QRAD] is assigned a value as shown. [AREA], is derived 
using Rule 9. [XRAD] is derived using Rule 23. 

Purpose: 

Calculate the heat load to be radiated from the surface, due to an 
enclosure. 

QenCL = ^lna z U\a^B m 

Premise: 

The qualifier "radiation effects to include" (Q2) has the value 
"enclosure radiation". The user is prompted with a screen containing 
all the possible radiation effects on the surface. The user selects one 
or multiple values, which are assigned to the qualifier. 

The qualifier "Enclosure radiation effects are to be” (Q21) has the 
value "calculated" The user is prompted with a screen with a place to 
enter a numeric value and two buttons. The screen satisfies the rule 
regardless of the input. Either the user will input a value for QENCL 
or will return a value of "calculated" which is assigned to Q21. 

Conclusion: 

The variable [QENCL] is assigned a value as shown. [AINST], area of 
instrument surface, and [EPSINST] , the emmissivity of the surface 
are input by the user with the screen [AINST]. [EBENC] is derived 
with assistance from Rule 18. 


A14 

I* RULE NUMBER: 16 
RULE: RAD HEAT LOAD 

IF: 

[XRAD] >0 
THEN: 

[QRAD1 IS GIVEN THE VALUE [XRAD]*[AREA1 


/* RULE NUMBER: 17 
RULE: ENCL HEAT LOAD 

IF: 

Radiation effects to include {enclosure radiationl 
and: Enclosure radiation effects are to be {calculated! 

THEN: 

[QENCL] IS GIVEN THE VALUE [AINST]*[EPSINST]*[EBENC] 



Purpose: 

Calculate the temperature and blackbody radiation of an enclosure. 
»p4 (li+it+ii+Ti+it+ii) 

I ENCL 7 


AND: 


oft, 


ENCL 


Premise: 

The qualifier "The temperature values are” (Q7) has the value 
"calculated". The qualifier is derived using Rule 24. 

Conclusion: 

[TENCL], the average calculated temperature of the sides of an 
enclosure, is assigned a value as shown. The variables [Tl], [T2], 
[T3], [T4], IT5], and [T6] are the temperatures of the respective 
surfaces. The values are calculated using Rule 24. 

[EBENCJ, the average blackbody radiation for the enclosure, is 
assigned a value as shown. [SIGMA] is defined in Rule 22. 


/* RULE NUMBER: 18 
RULE: ENCLOSURE TEMP 


IF: 

The temperature values are (calculated) 


THEN: 

[TENCL] IS GIVEN THE VALUE 
([Tl] A 4+[T2] A 4+[T3] A 4+lT4] A 4+[T5] A 4+[T61 A 4)/6 
and: [EBENC] IS GIVEN THE VALUE [SIGMA]*[TENCL] 
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Rules 


A16 


Purpose: 

Calculate the heat flux on the surface due to the effects of the sun. 

x sun = 


Premise: 

The qualifier "Radiation effects to include” (Q2) must include the value 
"sun". The user is prompted with the screen attached to Q2. The user 
is allowed to select one or all of the possible radiation effects listed. 
Q2 is the only multiple value qualifier in the expert system. 

Conclusion: 

The variable [XSUN1 is assigned the value as shown. The inference 
engine derives the value for [ALPHA], the surface absorptivity, using 
Rules 43-46. Even though Rules 38-41 can also derive values for 
[ALPHA], the rules were marked False as soon as Q1 was assigned a 
value from the main menu, so the inference engine will never try to 
"fire" them. Help Desk displays the screen attached to Q16 "assistance 
is". It displays variables including [ALPHA] and [MU]. The user can 
input a value or ask Help Desk for "assistance". Thus [MU] the solar 
view factor is defined as additional information from the screen. 

The inference engine derives the variable [S] the planet solar flux with 
Rule 47. 


/* RULE NUMBER: 19 
RULE: SUN FLUX 

IF: 

Radiation effects to include (sun) 

THEN: 

[XSUN1 IS GIVEN THE VALUE [ALPHA1*[S1*[MU1 



Purpose: 

Calculate the heat flux on the surface due to the effects of the planet's 
radiation. 

X PLANET = 

Premise: 

The qualifier "Radiation effects to include" (Q2) must include the value 
"planet". The user is prompted with the screen attached to Q2. The 
user is allowed to select one or all of the possible radiation effects 
listed. Q2 is the only multiple value qualifier in the expert system. 

Conclusion: 

The variable [XPLANET] is assigned the value as shown. The 
inference engine derives the value for [EPSILON] the surface 
emissivity using Rule 43-46. Even though Rules 38-41 can also 
derive values for [ALPHA], the rules were marked False as soon as Q1 
was assigned a value from the main menu, so the inference engine will 
never try to "fire" them. Help Desk displays the screen attached to 
Q16 "assistance is". It displays variables including [EPSILON] and 
[FE]. The user can input a value or ask Help Desk for "assistance”. 
Thus [FE] is defined as additional information from the screen. 

The inference engine derives the variable [E] the planet flux with Rule 
47. 


Rules 


/* RULE NUMBER: 20 
RULE: PLANET FLUX 

IF: 

Radiation effects to include [planet] 

THEN: 

[XPLANET] IS GIVEN THE VALUE [EPSILON]*[FEl*rE] 
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Rules 


A18 


Purpose: 

Calculate the heat flux on the surface due to the planet's albedo. 
\ albedo = ooFjSfl-p) 


Premise: 

The qualifier "Radiation effects to include” (Q2) must include the value 
"albedo". The user is prompted with the screen attached to Q2. The 
user is allowed to select one or all of the possible radiation effects 
listed. Q2 is the only multiple value qualifier in the expert system. 

Conclusion: 

The variable [XALBEDO] is assigned a value as shown. The inference 
engine derives the value for [ALPHA], the surface absorptivity, using 
Rules 43-46. Even though Rules 38-41 can also derive values for 
[ALPHA], the rules were marked False as soon as Q1 was assigned a 
value from the main menu, so the inference engine will never try to 
"fire” them. Help Desk displays the screen attached to Q16 
"assistance is". It displays variables including [ALPHA], [MU] and 
[FE], The user can input values or ask Help Desk for "assistance". 
Thus [MU] and [FE] are as additional information from the screen. 

The inference engine derives the variable [S], the planet solar flux, and 
[ALB], the planet albedo, with Rule 47. 


/* RULE NUMBER: 21 
RULE: ALBEDO FLUX 

IF: 

Radiation effects to include [albedo] 
THEN: 

[XALBEDO] IS GIVEN THE VALUE 
[ALPHA]*[ALB]*[FE]*[S]*(1-[MU]) 


Purpose: 

Define the value for the Stephen Boltzmann constant. EXSYS® does 
not have the ability to use scientific notation. Therefore, it has a 
problem with values such as 5.67E-8. The problem includes the 
inability to initialize the variable [SIGMA], The value must therefore 
be derived by a rule, as it is done here. 


/* RULE NUMBER: 22 
RULE: SIGMA 

IF: 

Value to calculate is (radiator area (given temperature)] 
or: Value to calculate is [surface temperaturel 


THEN: 

[SIGMA] IS GIVEN THE VALUE .0000000567 



Purpose: 

Calculate the radiative flux leaving the surface. 

X MD = eoF 4 

Premise: 

The qualifier "Value to calculate is" (Ql) has the value "radiator area 
(given temperature)” or "surface temperature". The user is prompted 
with the screen containing the main menu. When the user selects an 
analysis, the value is returned and assigned to Ql. 

Conclusion: 

The variable [XRAD] is assigned a value as shown. The inference 
engine derives the value for [EPSILON] the surface emissivity using 
Rule 43-46. [SIGMA], the Stephen-Boltzmann constant, is derived 
using Rule 22 and [TEMP], the temperature of the surface, is input by 
the user using the screen attached to the variable [TEMP]. 


/* RULE NUMBER: 23 
RULE: RAD FLUX 

IF: 

Value to calculate is (radiator area (given temperature)) 
or: Value to calculate is (surface temperature) 

THEN: 

[XRAD] IS GIVEN THE VALUE [EPSILON]*[SIGMA]*[TEMP] A 4 


Rules 
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Rules 


A20 


Purpose: 

Calculate the blackbody radiation generated by the enclosure. 
Premise: 

The qualifier "Value to calculate is” (Ql) has the value "enclosure 
temperatures and heat loads". The value is assigned with the main 
menu, 
or: 

The qualifier "Radiation effects to include" (Q2) has the value 
"enclosure radiation". The value is assigned with a screen containing 
all the possible radiation effects. 

The qualifier "Solar view factors are” (Q6) has the value "calculated", 
which is derived using Rules 28-31. The variable definitions of [MU1], 
IMU21, [MU31, [MU4], [MU51 and [MU61 are defined when deriving the 
qualifier value. 

The qualifier "Planet view factors calculations are" (Q15) has the value 
"finished". This is derived using Rules 32-33. The variable definitions 
of [FE1], [FE2], [FE3], [FE4], [FE5] and [FE6] are defined when 
deriving the qualifier value. 

Conclusion: 

QENCL.BAT is an external DOS batch file, which passes the 
information shown, and returns the six temperatures ([Tl] - [T6]) and 
heat loads ([Ql] - [Q6]>. 

The variables are taken from the expert system's known information. 
If any variable is not previously defined, the expert system will try to 
derive the value using rules in the rule base. If the values cannot be 
defined, the expert system will prompt the user, the expert system 
first looks for a screen attached to each of the variables. If a screen 
exists it is shown to the user, if not, the expert system displays its own 
screen using the defined prompt for the variable. 

In this way all the variables are defined. Instead of defining a screen 
for each variable, variables are grouped together. For example, [INS1] 
to [INS5] are all defined on the same screen. The screen will always 
be called by the first variable [INSl]. The knowledge engineer must 
make sure that every time any of the other variables are called that 
[INSl] is always listed first. 


/* RULE NUMBER: 24 
RULE: ENCLOSURE 

IF: 

Value to calculate is [enclosure temperatures and heat loads) 
or: Enclosure radiation effects are to be (calculated) 
and: Solar view factors are (calculated) 
and: Planet view factors calculations are (finished) 

THEN: 

X> RUN(QENCL.BAT [ALB],[S],[E],[X1,[Y],[Z], 
[FE1],[FE2),[FE3],[FE4],[FE5],[FE6], 
[ALPIR1],[ALPIR2],[ALPIR3],[ALPIR4],[ALPIR5],[ALPIR61, 
[ALPS1],[ALPS2],[ALPS3],[ALPS4],[ALPS5],[ALPS6], 
[MU1],(MU2],[MU3],[MU4],[MU5],[MU6], 
[NUMINS1,[INS1],[INS2],[INS3],[INS4],(INS5] fM /Z) 
and: The temperature values are (calculated) 



Purpose: 

Calculate the orbital temperature and heat loads on an enclosure. 
Rule 25 is for a planet oriented spacecraft, and Rule 26 is for a solar 
inertial oriented spacecraft. Otherwise they are identical. 

The analysis is similar to the enclosure analysis shown in Rule 24. 
Exceptions include the two qualifiers that must be defined in Rule 24 
and not here. The solar and planet view factors are calculated at every 
point in the orbit in the Fortran Code, so Help Desk is not required to 
calculate them first. 

Premise: 

The qualifier "Value to calculate is" (Ql) has the value "orbital average 
temperatures and heat loads”. The value is assigned with the main 
menu. 

The qualifier "spacecraft orientation is" (Q4) has the value "planet" for 
Rule 25 and "inertial" for Rule 26. The user is shown a screen 
containing two buttons. Once the user selects one, it is assigned to Q4. 

Conclusion: 

ORBIT.BAT is an external DOS batch file, which passes the 
information shown, and returns the orbital average information shown 
on the screen [REPORT5], 

The variables are taken from the expert system's known information. 
If any variable is not previously defined, the expert system will try to 
derive the value using rules in the rule base. If the values cannot be 
defined, the expert system will prompt the user, the expert system 
first looks for a screen attached to each of the variables. If a screen 
exists it is shown to the user. If not, the expert system displays its 
own screen using the defined prompt for the variable. 

In this way all the variables are defined. Instead of defining a screen 
for each variable, variables are grouped together. For example [INS1] 
to [INS5] are all defined on the same screen. The screen will always 
be called by the first variable [INSl]. The designer must make sure 
that every time any of the other variables are called that [INSl] is 
always listed first. 


/* RULE NUMBER: 25 
RULE: ORBITAL PLANET 

IF: 

Value to calculate is (orbital average temperatures and heat loads! 
and: The spacecraft orientation is (planet) 

THEN: 

X> RUN(ORBIT.bat [ALBU8UE], [GS1],IGS2],[GS3], 
riNCl,rOMEGA],fWOMEGAl,[ECCl,(SMAl,0, 

[X],[YJ,[ZJ,[REJ, 

[ALPIR1J,[ALPIR2],[ALPIR3],[ALPIR4],[ALPIR5],[ALPIR6], 

[ALPS1],[ALPS2],[ALPS3],[ALPS4],[ALPS5],[ALPS6], 

[NUMINS], [INSl], [INS2], [INS3], [INS4],[INS5], [NUMPTS] /M 
/Z) 

and: Orbital analysis is (finished) 


/* RULE NUMBER: 26 
RULE: ORBITAL INERTIAL 

IF: 

Value to calculate is (orbital average temperatures and heat loads) 
and: The spacecraft orientation is (inertial) 

THEN: 

X> RUN(ORBIT.bat [ALB], [S],[E], [GS1],[GS2],[GS3], 

[INC], [OMEGA], rWOMEGA],[ECC],[SMA],l, 

[X],[Y],[Z] ,[RE], 

[ALPIR1],[ALPIR2],[ALPIR3],[ALPIR4],[ALPIR5],[ALPIR6], 
[ALPS1], [ALPS2],[ALPS3],[ALPS4],[ALPS5],[ALPS6], 

[NUMINS], [INS1UINS2], [INS3], [INS4],[INS5], [NUMPTS] /M /Z) 
and: Orbital analysis is (finished) 
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Rules 


Purpose: 

Restrict the mode of entering the solar vector to global coordinates 
when performing orbital analysis. 

Premise: 

The qualifier "Value to calculate is" (Ql) has the value "orbital average 
temperatures and heat loads”. 

Conclusion: 

The qualifier "Orbit information will be” (Q3) is assigned the value 
"user specified global sun and Planet vector with orbit information 
(orbital elements)". See Rules 28-31. 


Purpose: 

Calculate the solar view factors for an enclosure. 

Premise: 

The qualifier "Orbit information will be" (Q3) has the value "user 
specified local sun and Planet vector". The sun vector is entered in 
local coordinates. 

Conclusion: 

Run the external DOS file "SOLAR.BAT", which passes the local sun 
vector [XS11, [XS2], and [XS31. The six solar view factors are returned, 
IMU1I, IMU2I, IMU3I, IMU41, |MU5| and |MU6|. 

The value "calculated" is assigned to the qualifier "Solar view factors 
are". 


A22 

/* RULE NUMBER: 27 
RULE: ORBIT-GLOBAL 

IF: 

Value to calculate is (orbital average temperatures and heat loads) 
THEN: 

Orbit information will be (user specified global sun and Planet vector 
with orbit information (orbital elements)) 


/* RULE NUMBER: 28 
RULE: SOLAR VECTOR 

IF: 

Orbit information will be (user specified local sun and planet vector) 
THEN: 

X> RUN(SOLAR.BAT [XS],[YS],[ZS] /M /Z) 
and: Solar view factors are (calculated) 



Purpose: 

Rules 29 and 30 have similar purposes: to calculate the local sun 
vector and the solar view factors. Rule 29 is for a spacecraft in planet 
orientation, and Rule 30 is for a spacecraft in solar inertial orientation. 

Premise: 

The qualifier "Orbit information will be" (Q3) has the value "user 
specified global sun and Planet vector with orbit information (orbital 
elements)”. The sun vector is entered in global coordinates. 

The qualifier "The spacecraft orientation is” (Q4) has the value 
"planet" or "inertial". The user is prompted with a screen that 
contains two buttons, one for planet orientation, and one for solar 
inertial orientation. 

Conclusion: 

Run the external DOS file "SUNVECT.BAT", which passes the global 
sun vector [GS1], [GS2], and [GS3J. The local sun vector, [XS1], [XS2J, 
and [XS3], and the six solar view factors are returned [MUll, [MU2], 
[MU3], [MU4], [MU5] and [MU6]. 

The qualifier "Solar view factors are" is assigned the value 
"calculated". 


/* RULE NUMBER: 29 
RULE: ORBIT ANGELS 

IF: 

Orbit information will be {user specified global sun and Planet vector 
with orbit information (orbital elements)) 
and: The spacecraft orientation is (planet) 

THEN: 

X> RUN(SUNVECT.BAT (GSU,[GS21,(GS31, 

IINCJ, IOMEGA), IWOMEGA) , IvO] /M /Z) 
and: Solar view factors are (calculated) 


/* RULE NUMBER: 30 
IF: 

Orbit information will be (user specified global sun and Planet vector 
with orbit information (orbital elements)) 
and: The spacecraft orientation is (inertial) 

THEN: 

X> RUN(SUNVECT.BAT (GS1UGS2UGS31, 

(INC), IOMEGA], 0,0 /M /Z) 
and: Solar view factors are (calculated) 
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Purpose: 

Define the defaults for the enclosure analysis. 


/* RULE NUMBER: 31 
RULE: ORBIT DEFAULTS 


Premise: 

The qualifier "Orbit information will be” (Q3) has the value "defaults". 


IF: 

Orbit information will be (defaults) 


Conclusion: 

Use the system defaults which are: 

semi-m^jor axis of 6785 km 

eccentricity of 0 deg 

true anomaly of 0 deg 

inclination of 0 deg 

line of nodes of 0 deg 

line of periapsis of 0 deg 

local sun vector of (-1,0,0) 

the solar view factors are (0,0,1, 0,0,0) 

The qualifier "Solar view factors are” is assigned the value "calculated". 


THEN: 

[SMA] IS GIVEN THE VALUE 6785.0 
and: [ECC] IS GIVEN THE VALUE 0 

and: [V0] IS GIVEN THE VALUE 0.0 

and: [INC1 IS GIVEN THE VALUE 0 

and: [OMEGA] IS GIVEN THE VALUE 0 

and: [WOMEGA] IS GIVEN THE VALUE 0 

and: [XS] IS GIVEN THE VALUE -1.0 

and: [YS] IS GIVEN THE VALUE 0.0 

and: [ZS] IS GIVEN THE VALUE 0.0 

and: [MU1] IS GIVEN THE VALUE 0 

and: [MU2] IS GIVEN THE VALUE 0 

and: [MU3] IS GIVEN THE VALUE 1 

and: [MU4] IS GIVEN THE VALUE 0 

and: [MU5] IS GIVEN THE VALUE 0 

and: [MU6] IS GIVEN THE VALUE 0 

and: Solar view factors are (calculated) 



Purpose: 

Finish the planet view factors calculations after the user decides to not 
modify the calculated planet view factors. 

Premise: 

The qualifier "planet view factors are" (Q8) has the value "calculated”, 
which is derived from Rules 34-36. 

The qualifier "planet view factors are to be" (Q14) has the value 'left 
alone". The user is allowed to modify the planet view factors 
calculated by Help Desk. A screen attached to Q14 is displayed asking 
the question if the user wants to modify the planet view factors. If the 
user selects "no", the value "left alone" is returned. 

Conclusion: 

The qualifier "Planet view factors calculations are" is assigned the 
value "finished". 


t* RULE NUMBER: 32 
RULE: LEAVE EVF 

IF: 

Planet view factors are (calculated) 
and: Planet view factors are to be (left alone) 

THEN: 

Planet view factors calculations are (finished) 


Purpose: 

Allows the user to modify the calculated planet view factors. 

Premise: 

The qualifier "planet view factors are" (Q8) has the value "calculated", 
which is derived from Rules 34-36. 

The qualifier "planet view factors are to be" (Q14) has the value 
"changed". The user is allowed to modify the planet view factors 
calculated by Help Desk. A screen attached to Q14 is displayed asking 
the question if the user wants to modify the planet view factors. If the 
user selects "yes”, the value "changed” is returned. 

Conclusion: 

The planet view factors are displayed, and the modified values are 
assigned to the variables [FE1], [FE2], [FE3], IFE4],. [FE5], and [FE6] 

The qualifier "Planet view factors calculations are” is assigned the 
value "finished". 


/* RULE NUMBER: 33 
RULE: EVF CHANGED 

IF: 

Planet view factors are (calculated) 
and: Planet view factors are to be (changed) 


THEN: 

[FE1] IS GIVEN THE VALUE [FE1U] 
and: [FE2] IS GIVEN THE VALUE (FE2U) 

and: [FE3] IS GIVEN THE VALUE [FE3U1 

and: [FE4] IS GIVEN THE VALUE [FE4U1 

and: [FE51 IS GIVEN THE VALUE (FE5UI 

and: (FE6) IS GIVEN THE VALUE [FE6U] 

and: Planet view factors calculations are (finished) 


Rules 
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Purpose: 

Calculate the planet view factors for an inertially oriented spacecraft. 
Premise: 

The qualifier "The spacecraft orientation is" (Q4) has the value 
"inertial". 

The qualifier "The spacecraft orbit is” (Q5) has the value "LEO (Low 
Earth Orbit)". 

Conclusion: 

Run the external program "PLANET.BAT". It returns the six planet 
view factor, [FE1], {FE2], [FE31, [FE4],. [FE5], and [FE61. 

The qualifier "Planet view factors are” is assigned the value 
"calculated". 


Purpose: 

Calculate the planet view factors for an planet oriented spacecraft. 
Premise: 

The qualifier "The spacecraft orientation is” (Q4) has the value 
"planet". 

The qualifier "The spacecraft orbit is" (Q5) has the value "LEO (Low 
Earth Orbit)". 

Conclusion: 

Assign the values to the six planet view factor, (FE1), [FE2], IFE3], 
[FE4],. [FE5], and fFE6], In the variable definition of [FET] and 
[FEN] is a call statement to an external DOS executable, 
"FETANOR.BAT". The file calculates the normal and tangential 
components of the solar view factor. 


I* RULE NUMBER: 34 
RULE: INERTIAL LEO 

IF: 

The spacecraft orientation is (inertial) 
and: The spacecraft orbit is (LEO (Low Earth Orbit)) 

THEN: 

X> RUN(PLANET.BAT (INC1, (ECCl.fWOMEGAl ,[SMA1,(V0],(RE] 
/M/Z) 

and: Planet view factors are (calculated) 


/* RULE NUMBER: 35 
RULE: PLANET LEO 

IF: 

The spacecraft orientation is (planet) 
and: The spacecraft orbit is (LEO (Low Earth Orbit)) 

THEN: 

[FE1] IS GIVEN THE VALUE [FET] 
and: [FE2] IS GIVEN THE VALUE [FET] 

and: [FE3] IS GIVEN THE VALUE 0.0 

and: (FE4] IS GIVEN THE VALUE [FET] 

and: [FE5] IS GIVEN THE VALUE [FEN] 

and: [FE6] IS GIVEN THE VALUE [FET] 

and: Planet view factors are (calculated) 


The qualifier "Planet view factors are" is assigned the value 
"calculated". 



Purpose: 

Assign the planet view factors for a Geosynchronous spacecraft. 
Premise: 

The qualifier "The spacecraft orbit is” (Q5) has the value ”GEO 
(Geosynchronous Earth Orbit)”. 

Conclusion: 

Assign the values to the six planet view factor, [FE1], [FE2], [FE3], 
[FE4],. [FE5], and [FE61. In GEO, the planet view factors are 
effectively zero. 

Assign the value of 35744 km to the variable [SMA], the semi-mqjor 
axis. 

The qualifier "Planet view factors are" is assigned the value 
"calculated". 


Purpose: 

Prevent the user from being prompted for insulated side information, 
when the user has already defined there are no insulated sides. 

Premise: 

The variable [NUMINS] has the value 0. There are no insulated sides 
during the analysis. 

Conclusion: 

Assign the value 0.0 to the variable [INS1], The screen containing the 
additional information on the insulated sides is attached to the 
variable [INS1]. If the variable has a value, it does not prompt the 
user for the information. 


/* RULE NUMBER: 36 
RULE: GEO 

IF: 

The spacecraft orbit is {GEO (Geosynchronous Earth Orbit)) 
THEN: 

[FE1] IS GIVEN THE VALUE 0.0 
and: [FE2] IS GIVEN THE VALUE 0.0 

and: (FE3] IS GIVEN THE VALUE 0.0 

and: [FE4] IS GIVEN THE VALUE 0.0 

and: [FE5] IS GIVEN THE VALUE 0.0 

and: [FE6] IS GIVEN THE VALUE 0.0 

and: [SMA] IS GIVEN THE VALUE 35744 

and: Planet view factors are (calculated) 


/* RULE NUMBER: 37 
RULE: NUMINS 

IF: 

[NUMINS] = 0 
THEN: 

fINSI] IS GIVEN THE VALUE 0.0 
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Purpose: 

Rules 38-41 have the same purpose: to define the surface properties; 
solar absorptivity and emissivity. 

Premise: 

The premise clauses for the four mles are similar, will all possible 
combinations on two qualifiers. 

The qualifier "Value to calculate" (Ql) has the value "surface 
properties". The value is selected from the Main menu. 

The qualifier "The performance time is” (Q9) may have the value 
"BOLT or” EOL". The user is shown a screen with two buttons. The 
button chosen returns the value to the qualifier. 

The qualifier "The known surface information is” (Q10) may have the 
value "surface use" or "surface application". Another two button screen 
is displayed for the user to choose from. 

Conclusion: 

The variable [EPSILON] is assigned the value of the variable as 
shown. These variables have frame calls to the frame 
"SURFACE.TXT"stated in the variable definitions . Frames are 
database calls to the small ASCII database called a frame. 

The variable [ALPHA], [MATERIAL] or [APPLICATION] are defined 
as shown with frame calls to the same file. 


/* RULE NUMBER: 38 
RULE: FRAME EPS BOL 

IF: 

Value to calculate is (surface properties) 
and: the performance time is (BOL (beginning of life)) 
and: the known surface information is (surface use) 

THEN: 

[EPSILON] IS GIVEN THE VALUE [A_EPSFB] 
and: [ALPHA] IS GIVEN THE VALUE [A.ABSFB] 

and: [MATERIAL] IS GIVEN THE VALUE FRAMECsurface.txt", 

$TYPICAL_APPLICATION$==[APPLICATION], $SURFACE$) 
and: Surface properties are (finished) 


/* RULE NUMBER: 39 
RULE: FRAME EPS EOL 

IF: 

Value to calculate is (surface properties) 
and: the performance time is (EOL (end of life)) 
and: the known surface information is (surface use) 

THEN: 

[EPSILON] IS GIVEN THE VALUE [A_EPSFE] 
and: [ALPHA] IS GIVEN THE VALUE [A.ABSFE] 

and: [MATERIAL] IS GIVEN THE VALUE FRAMECsurface.txt", 

$TYPICAL_APPLICATION$==[ APPLICATION], $SURFACE$) 
and: Surface properties are (finished) 



/* RULE NUMBER: 40 


IF: 

Value to calculate is (surface properties) 
and: the performance time is (EOL (end of life)) 
and: the known surface information is (surface material) 


THEN: 

fEPSILONl IS GIVEN THE VALUE [M.EPSFEl 
and: (ALPHA] IS GIVEN THE VALUE [M_ABSFE| 

and: [APPLICATION] IS GIVEN THE VALUE FRAMECsurface.txt", 

$SURFACE$==[MATERIAL], $TYPICAL_APPLICATION$) 
and: Surface properties are (finished) 


/* RULE NUMBER: 41 


IF: 

Value to calculate is (surface properties) 
and: the performance time is (BOL (beginning of life)) 
and: the known surface information is (surface material) 

THEN: 

(EPSILON] IS GIVEN THE VALUE [M.EPSFB] 
and: [ALPHA] IS GIVEN THE VALUE [M_ABSFB] 

and: (APPLICATION] IS GIVEN THE VALUE FRAMECsurface.txt", 

$SURFACE$==[MATERIAL], $TYPICAL_APPLICATION$) 
and: Surface properties are (finished) 


Rules 


A29 



Rules 


A30 


Purpose: 

To allow the user to modify the surface properties calculated by Help 
Desk. 

Premise: 

The qualifier "assistance is" (Q) has the value "needed". 

The qualifier "surface properties are” (Q) has the value "calculated". 
This is derived using Rules 43-46. 

The qualifier "Change of surface properties is" (Q) has the value " 
requested". The value is input by the user. The user is shown a 
screen containing the solution to the surface properties analysis (Rules 
43-46). The screen also contains two buttons one to continue, and one 
to modify the values. The modify button returns the value "requested” 
to the qualifier. 

Conclusion: 

A screen is displayed allowing the user to modify the emissivity and 
absorptivity. 


f* RULE NUMBER: 42 
RULE: CHANGE SP 

IF: 

assistance is (needed) 
and: Surface properties are (calculated) 
and: Change of surface properties is (requested) 

THEN: 

[EPSILON] IS GIVEN THE VALUE [CEPS] 
and: [ALPHA] IS GIVEN THE VALUE [CALP] 



Purpose: 

Rules 43-46 have the same purpose: to define the surface properties, 
solar absorptivity and emissivity. During a surface analysis, the user 
may ask the expert system for assistance in determining the surface 
properties, emissivity and absorptivity. 

Premise: 

The premise clauses for the four rules are similar, will all possible 
combinations on two qualifiers. 

The qualifier "Assistance is" (Q16) has the value "needed”, which is 
derived using the screen Q16. At the bottom of the screen is a button 
called "assistance". When this is pressed, the value "needed" is 
returned. 

The qualifier "The performance time is" (Q9) may have the value 
"BOL" or" EOL". The user is shown a screen with two buttons. The 
button chosen returns the value to the qualifier. 

The qualifier "The known surface information is” (Q10) may have the 
value "surface use" or "surface application". Another two button screen 
is displayed for the user to choose from. 


Conclusion: 

The variable [EPSILON] is assigned the value of the variable as 
shown. These variables have frame calls to the frame 
"SURFACE.TXT" in the variable definitions. They are database calls 
to the small ASCII database called a frame. 

The variable [ALPHA], [MATERIAL] or [APPLICATION] are defined 
as shown with frame calls to the same file. 


/* RULE NUMBER: 43 
RULE: ASSIST APP BOL 

IF: 

assistance is (needed) 

and: the performance time is {BOL (beginning of life)] 
and: the known surface information is (surface use) 

THEN: 

[EPSILON] IS GIVEN THE VALUE [A.EPSFB] 
and: [ALPHA] IS GIVEN THE VALUE [A.ABSFB] 

and: [MATERIAL] IS GIVEN THE VALUE FRAMECsurface.txt", 

$TYPICAL_APPLICATION$==[APPLICATION], $SURFACE$) 
and: Surface properties are {calculated! 


/* RULE NUMBER: 44 
RULE: ASSIST APP EOL 

IF: 

assistance is (needed) 

and: the performance time is {EOL (end of life)) 
and: the known surface information is {surface use] 

THEN: 

[EPSILON] IS GIVEN THE VALUE [A.EPSFE] 
and: [ALPHA] IS GIVEN THE VALUE [A.ABSFE] 

and: [MATERIAL] IS GIVEN THE VALUE FRAMECsurface.txt", 

$TYPICAL_APPLICATION$==[APPLICATION], $SURFACE$) 
and: Surface properties are (calculated) 
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/* RULE NUMBER: 45 
RULE: ASSIST MAT BOL 


IF: 

assistance is (needed) 

and: the performance time is (EOL (end of life)) 

and: the known surface information is (surface material) 

THEN: 

[EPSILON] IS GIVEN THE VALUE [M.EPSFE] 
and: [ALPHA) IS GIVEN THE VALUE [M.ABSFE] 

and: [APPLICATION] IS GIVEN THE VALUE FRAMECsurface.txt", 

$SURFACE$==[MATERIAL], $TYPICAL_APPLICATION$) 
and: Surface properties are (calculated) 


/* RULE NUMBER: 46 
RULE: ASSIST MAT EOL 


IF: 

assistance is (needed) 

and: the performance time is (BOL (beginning of life)) 

and: the known surface information is (surface material) 

THEN: 

[EPSILON] IS GIVEN THE VALUE [M_EPSFB] 
and: [ALPHA] IS GIVEN THE VALUE [M_ABSFB] 

and: [APPLICATION] IS GIVEN THE VALUE FRAMECsurface.txt", 

$SURFACE$==[MATERIAL], $TYPICAL_APPLICATION$) 
and: Surface properties are (calculated) 



Purpose: 

Define the values for solar flux, planet flux, planet albedo and planet 
radius. During an analysis, the user may select to run the analysis 
around a different planet than the Earth. The user can input the 
known values for the fluxes, the albedo and the radius of the planet, 
or the user can select the planet from a list. The system will use the 
defaults stored in the frame "PLANET.TXT". 

Premise: 

The qualifier "A different planet is” (Q17) has the value ’'requested”, 
which is derived using the screen Q17. At the bottom of the screen is 
a button called "different planet”. The user has the option to input the 
values of any planet in the edit boxes supplied. If the values are 
unknown, the user may select the "different planet" button. The value 
"requested" is returned. The expert system will then prompt the user 
with a list of available planets. The list is attached to the variable 
[PLANET]. 

Conclusion: 

The values for the variables are taken from the frame "PLANET.TXT”. 
It is located in the same directory as the expert system files. The 
values returned from the frame are assigned to the variables and are 
used in the rest of the analysis. 
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I* RULE NUMBER: 47 
RULE: PLANET FRAME 


a different planet is (requested) 

THEN: 

[S] IS GIVEN THE VALUE 

FRAMECPLANET.TXT" ,$PLANET$==[PLANET],#SOLAR_FLUX#) 
and: [E] IS GIVEN THE VALUE 

FRAMECPLANET.TXT", $PLANET$==[PLANET],#PLANET_FLUX#) 
and: [ALB] IS GIVEN THE VALUE 

FRAMECPLANET.TXT", $PLANET$==[PLANET],#ALBEDO#) 
and: [RE] IS GIVEN THE VALUE 

FRAMECPLANET.TXT", $PLANET$==[PLANET],#RADIUS#) 
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Qualifiers 



f* Qualifier 1 
Q> Value to calculate is 

V> radiator area (given temperature) 

V> surface temperature 
V> enclosure temperatures and heat loads 
V> orbital average temperatures and heat 
loads 

V> surface properties 
V> none of the above 

Name: CALCULATION OPT 
Maximum acceptable = 1 

/* Qualifier 2 

Q> Radiation effects to include 

V> sun 
V> planet 
V> albedo 

V> enclosure radiation 

V> none (just radiate internal power) 

Name: RADIATION EFFEC 

/* Qualifier 3 

Q> Orbit information will be 

V> user specified local sun and Planet vector 
V> user specified global sun and Planet vector 
with orbit information (orbital elements) 
V> defaults 

V> User specified Beta angle 

Name: ORBIT INFO 
Maximum acceptable = 1 

/* Qualifier 4 

Q> The spacecraft orientation is 

V> planet 
V> inertial 

Name: PLANET/INERTIAL 
Maximum acceptable = 1 


/* Qualifier 5 

Q> The spacecraft orbit is 

V> LEO (Low Earth Orbit) 

V> GEO (Geosynchronous Earth Orbit) 

Name: LEO/GEO 
Maximum acceptable = 1 

f* Qualifier 6 

Q> Solar view factors are 

V> calculated 
V> unknown 

Name: SOLAR VIEW FACT 
Default value = 2 
Maximum acceptable = 1 

/* Qualifier 7 

Q> The temperature values are 

V> calculated 
V> unknown 

Name: TEMPERATURE 
Default value = 2 
Maximum acceptable = 1 

/* Qualifier 8 

Q> Planet view factors are 

V> calculated 
V> unknown 

Name: PLANET VIEW FAC 
Maximum acceptable = 1 

/* Qualifier 9 

Q> the performance time is 

V> BOL (beginning of life) 

V> EOL (end of life) 

Name: BOL/EOL 
Maximum acceptable = 1 


Qualifiers 
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/* Qualifier 10 

Q> the known surface information is 

V> surface use 
V> surface material 

Name: SURFACE INFO 
Maximum acceptable = 1 

/* Qualifier 11 
Q> Orbital analysis is 

V> finished 
V> undone 

Name: ORBITAL ANAL 
Maximum acceptable = 1 

/* Qualifier 12 

Q> Surface temperature to be calculated is 

V> radiator temperature given area 
V> environmental temperature 

Name: SURFACE TEMP 
Maximum acceptable = 1 

/* Qualifier 13 
Q> Surface properties are 

V> calculated 
V> unknown 

Name: SURF PROPS CAL 
/* Qualifier 14 

Q> Planet view factors are to be 

V> changed 
V> left alone 

Name: CHANGE EVF 


/* Qualifier 15 

Q> Planet view factors calculations are 

V> finished 
V> undecided 


/* Qualifier 16 
Q> assistance is 

V> needed 
V> not needed 

Name: ASSISTANCE 

/* Qualifier 17 

Q> a different planet is 

V> requested 
V> not requested 

Name: DIFFERENT PLANE 


/* Qualifier 18 

Q> Change of surface properties is 

V> requested 
V> not needed 

Name: CHANGE SP 


/* Qualifier 19 

Q> The report information is to be 

V> appended to the current file 
V> written to a new file 
V> no file output 

Name: APPEND TO FILE 
Default value = 3 
Maximum acceptable = 1 

/* Qualifier 20 
Q> Thermal analysis is 

V> completed 
V> undone 

Name: THERMAL DONE 
f* Qualifier 21 

Q> Enclosure radiation effects are to 

V> calculated 
V> user input 


Name: PLANET VF FIN 


Name: CAL RAD EFF 



Appendix C 
Screens 
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EXSYS Professions! 

Choose the analysis to perfo 

rm: 


surface area {given temperature) • •!: 




surface temperature 





[ 

enclosure temperature and heat loads! • 


1 

pforbital temperature and heat tgdjdsl ' 



[ 

: surface properties 


r 

IIBSB 



Figure 23 ~Q1 


-Q1 

-SET.WIN 88, 0, 517, 366 
-LINE.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-TEXTSIZE LARGE 

-CURSET 34, 39 

Choose the analysis to perform: 

-BUTTON "surface area (given temperature)”, "ql 1”, 60, 60, 395, 40 -I 

-BUTTON "surface temperature”, "ql 2", 60, 110, 395, 40 -I 

-BUTTON "enclosure temperatures and heat loads", "ql 3”, 60, 160, 395, 40 -I 

-BUTTON "orbital temperatures and heat loads", "ql 4”, 60, 210, 395, 40 -I 

-BUTTON "surface properties”, "ql 5", 60, 260, 395, 40 -I 

-BUTTON "EXIT", "Ql 6”, 60, 310, 395, 40 -I 

-END 


Screens 
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■ • ,.K£XSYS,gro ' 

Radiation effects to include 

Choose one. or more of: the fottowing 

B;:|| □ Sun 111! 

' • □ Planet | |||i fl 

□ Albedo | 

□ Enclosure Radiation ;v(_V 

llli □ internal power only 

%yM-. continue j.i i 


Figure 24 -Q2 


~Q2 

-SET.WIN 88, 0, 375, 335 
~LINE_COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 70, 25 
Radiation effects to include: 

-CURSET 32, 59 

Choose one or more of the following 
-BUTTON "Sun”, ”q2 1”, 65, 77, 238, 35 -C 
-BUTTON "Planet", "q2 2”, 65, 120, 238, 35 -C 
-BUTTON "Albedo", ”q2 3”, 65, 160, 238, 35 -C 
-BUTTON "Enclosure Radiation", "q2 4", 65, 199, 238, 35 -C 
-BUTTON "internal power only", "q2 5”, 65, 240, 238, 35 -C 
-BUTTON "CONTINUE", 127, 289, 100, 33 -I 
-END 
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Screens 



Illllllf EXSYS ! 'Prof essionaiglMil 

Orbit information will be: 


$:£§• ond Rlanl^yecior ■ 

global Sun vector with orbit information (orbital elements) 


;4i§le gif avatlabiglj 


Figure 25 -Q3 


~Q3 

~SET_WIN 88, 0, 518, 252 
-BACKGROUND CYAN 
-CLS 

-TEXTSIZE LARGE 
-CURSET 86, 33 
Orbit information will be: 

-BUTTON 'local sun and Planet vector”, "Q3 1”, 20, 47, 476, 35 -I 

-BUTTON "global sun vector with orbit information (orbital elements)”, ”Q3 2”, 20, 90, 476, 35 -I 
-BUTTON "beta angle (not available)”, 20, 130, 476, 35 
-BUTTON "defaults”, ”Q3 3”, 20, 170, 476, 35 -I 
-END 


Screens 
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EXSVS Professional l;- :/ 
Spocecraft orientation is: 

' . Planet 


.Inertial 


Figure 26 -Q4 


-Q4 

-SET_WIN 88, 0, 249, 140 
-LINE.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 12, 30 
Spacecraft orientation is: 

-BUTTON "Planet", ”Q4 1", 69, 49, 100, 35 -I 
-BUTTON "Inertial", "Q4 2", 68, 93, 100, 35 -I 
-END 
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life': :..a :£X : 5YS Professional lliiilll 


Spacecraft orbit is: 


[ 

LEG (Low Earth; 0rbli}. 


[ 

GEO (Geosynchronous};:;:/ 


Figure 27 ~Q5 


-Q5 

~SET_WIN 88, 0, 306, 143 
-LINE.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 67, 28 
Spacecraft orbit is: 

-BUTTON "LEO (Low Earth Orbit)", ”Q5 1”, 27, 45, 243, 35 -I 
-BUTTON "GEO (Geosynchronous)”, "Q5 2", 26, 92, 243, 35 -I 
-END 


Screens 




EXSYS Rrofessionol 


'Choose the performance time for the surface: 
j BQL (Beginning GlfLilelf: 

[IrEOl (End Of Life] | 


Figure 28 ~Q9 


-Q9 

-SET.WIN 88, 0, 356, 143 
~LINE_COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 19, 31 

The perfomance time for the surface: 

-BUTTON "BOL (Beginning Of Life)", ”q9 1”, 62, 44, 229, 35 -I 
-BUTTON "EOL (End Of Life)", ”q9 2", 62, 86, 229, 35 -I 
-END 
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v. EXSYS Professtonof 

The known surface information 

is: 


. Surface material 






• Surface opptkation 






Figure 29 ~Q10 


-Q10 

~SET_WIN 88, 0, 313, 158 
-LINE.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 11, 32 

The known surface information is: 

-BUTTON "Use of the surface", "qlO 1", 34, 95, 240, 35 -I 
-BUTTON "Surface material", "qlO 2", 34, 52, 240, 35 -I 
-END 


Screens 



•:f£^SYS ::: f?raf essiono! |j 
Choose the temperature analysis 


Surface temperature {given area} 


•.environmental temperature: 


Figure 30 -Q12 


-Q12 

-SET.WIN 88, 0, 353, 145 
-LINE.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 27, 29 

Choose the temperature analysis 

-BUTTON "Surface temperature (given area)”, ”ql2 1”, 35, 52, 274, 35 -I 
-BUTTON "Environmental temperature", ”ql2 2", 35, 97, 275, 35 -I 
-END 
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. •••> ilLlXSYS Professionol^^'-;':^-- ; • 

Do you wont to modify the Planet View factors: 




Figure 31 -Q14 


-Q14 

-SET.WIN 88, 0, 312, 138 
-LINE.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 22, 29 

modify the planet view factors? 

-BUTTON "yes", "Q14 1”, 116, 47, 75, 35 -I 
-BUTTON "no”, ”Q14 2”, 116, 96, 75, 35 -I 
-END 


Screens 
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EX SYilfiiolesstancl 


Surface Constants 
Solar view factor: 

Planet view factor: 
Surface emissivity: 
Surface absorptivity 


jjj 

Figure 32 ~Q16 


; assistance 



-Q16 

~SET_WIN 88, 0, 380, 280 
-LINE.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 110, 35 
Surface Constants 
-CURSET 30, 75 
Solar view factor: 

-CURSET 30, 115 
Planet view factor: 

-CURSET 30, 155 
Emissivity: 

-CURSET 30, 195 
Absorptivity: 

-EDIT ”[MU]",.225, 50, 103, 35 -M:l 
-EDIT IFE]", 225, 90, 103, 35 -M:l 
-EDIT "[EPSILON]", 225, 130, 103, 35 -M:l 
-EDIT "[ALPHA]”, 225, 170, 103, 35 -M:l 
-BUTTON "OK”, ”Q16 2", 261, 228, 75, 35 -I 
-BUTTON "assistance", "ql6 1", 38, 228, 158, 35 -I 
-EDIT.FILL ALPSCR.DAT 
-END 
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EXSYS .Professional /• 

Planet Constants . 

Albedo coefficient: 
Solar Flux: W/m-2 
Planet Flux: W/m-2 
Planet radius: km 

.3 

1356 

236 

6378 



different J[ ;:QK | 



Figure 33 -Q17 


-Q17 

-SET.WIN 88, 0, 360, 271 
-LINE.COLOR BLACK 
-FILL.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 102, 30 
Planet Constants 
-CURSET 10, 70 
Albedo coefficient: 

-CURSET 10, 110 
Solar Flux: W/m A 2 
-CURSET 10, 150 
Planet Flux: W/m A 2 
-CURSET 10, 190 
Planet radius: km 

-EDIT "[ALB]”, 225, 45, 103, 35 -M:l 

-EDIT "[S]", 225, 85, 103, 35 -M:l 

-EDIT "[E]”, 225, 125, 103, 35 -M:l 

-EDIT "[RE]”, 225, 165, 103, 35 -M:l 

-BUTTON "OK", "Q17 2”, 241, 217, 75, 35 -I 

-BUTTON "different planet", ”ql7 1”, 23, 216, 158, 35 -I 

-EDIT.FILL ALBSCR.DAT 

-END 


Screens 
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Assistance Results 

The surface is used for interior structure 
The surface is made of black paint 

the performance time is BOL ' ' 
Emissiviiy . 9 . . 

Solor absorptivity .9 -•= .• • .• . 

Ilill l Modify: 


C animal 


Figure 34 -Q18 


-Q18 

-SET.WIN 0, 0, 584, 266 
-LINE.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 182, 29 
Assitance Results 
-CURSET 54, 62 

The surface is used for [[APPLICATION]] 
-CURSET 54, 96 

The surface is made of [[MATERIAL]] 

-CURSET 54, 132 
[f*q9]] 

-CURSET 54, 168 
Emissivity [[EPSILON]] 

-CURSET 54, 202 

Solar absorptivity [[ALPHA]] 

-BUTTON "Modify", "Q18 1", 113, 213, 106, 35 -I 
-BUTTON "Continue", "Q18 2", 304, 213, 106, 35 -I 
-END 
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r ' '• £X SYS Pro f es s eo n o i x : : i ||; ; 

Enter enclosure radiation heat load (W): 

0 



1 

^pafcuiateft 


:;;;OK 



Figure 35 -Q21 


-Q21 

-SET.WIN 88, 0, 539, 120 
-LINE.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 16, 35 

Enter enclosure radiation heat load (W): 

-EDIT "[QENCL]", 396, 11, 103, 35 -M:l 
-BUTTON "OK”, "Q21 2 ", 330, 66, 75, 35 -I 
-BUTTON "calculate", "Q21 1”, 118, 66, 158, 35 -I 
-EDIT.FILL QENCLSCR.DAT 
-END 


Screens 
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EXSYS Professional '• j 

Enclosure Instrument Paran 

Instrument inside the enclosu 

Surface area of enclosure instrument (m-2) 

! 

Emissivity of enclosure instrument 

PM 

neters 

re 

.01 

1 



Figure 36 -[AINST] 


-[AINST] 

-SET.WIN 88, 0, 546, 176 
-LINE.COLOR BLACK 
-FILL.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 137, 30 
Enclosure Instrument Parameters 
-TEXTSIZE SMALL 
-CURSET 169, 51 
Instrument inside the enclosure 
-TEXTSIZE NORMAL 
-CURSET 29, 77 

Surface area of enclosure instrument (m A 2) 
-CURSET 29, 114 
Emissivity of enclosure instrument 
-EDIT "[AINST]”, 426, 50, 100, 35 -M:l 
-EDIT "[EPSINST]", 425, 92, 100, 35 -M:l 
-BUTTON "OK", ”", 232, 128, 75, 35 -I 
-EDIT.FILL AINSTSCR.DAT 
-END 
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[ C<SYS Professional ... til 


Choose the surface use: 



interior structure ,v . 


|. amenno reflector 




:: w non h ; :pone j- • radiator /.s; . 

LV south pond rodiaior 






• • solor ppnei • 



[: ontenno structure ,i';| 


!|||F thermal, : u .^i 

apogee motor-thermal shiel#::; 

pro p.ellont •jnsutotiorv : - : "j; 




Figure 37 -[APPLICATION] 


-[APPLICATION] 

-SET.WIN 88, 0, 319, 359 
-LINE.COLOR BLACK 
-FILL.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 30, 22 

Choose the surface application: 

-BUTTON "interior structure”, "[APPLICATION] 'interior structure'", 25, 30, 272, 25 -I 
-BUTTON "antenna reflector", "[APPLICATION] 'antenna reflector’", 25, 60, 272, 25 -I 
-BUTTON "north panel radiator", "[APPLICATION] 'north panel radiator’", 25, 90, 272, 25 -I 
-BUTTON "south panel radiator", "[APPLICATION] 'south panel radiator’", 25, 120, 272, 25 -I 
-BUTTON "solar panel", "[APPLICATION] solar panel’", 25, 150, 272, 25 -I 
-BUTTON "antenna structure", "[APPLICATION] 'antenna structure’", 25, 180, 272, 25 -I 
-BUTTON "thermal insulation", "[APPLICATION] 'thermal insulation”', 25, 210, 272, 25 -I 
-BUTTON "apogee motor thermal shield”, "[APPLICATION] 'apogee motor thermal shield’", 25, 240, 
272, 25 -I 

-BUTTON "propellant insulation”, "[APPLICATION] 'propellant insulation"’, 25, 270, 272, 25 -I 
-BUTTON "surface insulation", "[APPLICATION] 'surface insulation”', 25, 300, 272, 25 -I 
-END 


Screens 
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£XSYS; ; |irofessional f : j 

Enter the area of radiating surface (m-2) 

i 



■ill 



Figure 38 -[AREA] 


-[AREA] 

~SET_WIN 88, 0, 515, 100 
~LINE_COLOR BLACK 
-FILL_COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 17, 37 

Enter the area of radiating surface (m A 2) 
-EDIT ’’[AREA]", 400, 12, 103, 35 -M:l 
-BUTTON "OK", 203, 54, 75, 35 -I 

-EDIT.FILL AREASCR.DAT 
-END 
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•■: '. : £XSYS Professional 

Surface Properiies 


Emissivity: .9 


Absorptivity: .9 

PHI 




Figure 39 -[CEPS] 


-[CEPS] 

-SET.WIN 88, 0, 306, 197 
~LINE_COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 79, 31 
Surface Properties 
-CURSET 30, 75 
Emissivity: 

-CURSET 30, 115 
Absorptivity: 

-EDIT "[CEPS]", 161, 50, 103, 35 -L 
[[EPSILON]] 

-LISTEND 

-EDIT ”[CALP]”, 161, 90, 103, 35 -L 
[[ALPHA]] 

-LISTEND 

-BUTTON "OK”, 121, 147, 75, 35 -I 
-END 


Screens 




•• EXSYS Professional ==• .. 

Specify Spacecraft Pc 
true anomaly at epoch (deg) 
semi-major axis (km) 
eccentricity j 

| OK : ip 

)sition 

o ... 

6785.0 

o 



Figure 40 Spacecraft Position screen 


~[ECC] 

-SET.WIN 88, 0, 390, 224 
-BACKGROUND CYAN 
-CLS 

-CURSET 91, 32 
Specify Spacecraft Position 
-CURSET 14, 75 
true anomaly at epoch (deg) 
-CURSET 14, 115 
semi-maoor axis (km) 

-CURSET 14, 155 
eccentricity 

-EDIT ”[v0]", 280, 50, 100, 35 -M:l 
-EDIT "[SMA]", 280, 90, 100, 35 -M:l 
-EDIT "[ECC]”, 280, 130, 100, 35 -M:l 
-BUTTON "OK', 160, 177, 75, 36 -I 

-EDIT_FILL ECCSCRJDAT 
-END 


C20 


Screens 




•7727 •• £XSYS Profegsronol 

Specify planet view factor information 

side 1 ( — Z) 

side 2 (-Y) 

side 3 ( — X) 
side 4 (+Y) 

side 5 ( + X) 
side 6 ( + Z) 

F-(^l 

Figure 41 ~[FE1U] 


.287173 


.287173 


0 


.287173 


.883693 


.287173 


-IFE1U] 

~SET_WIN 88, 0, 334, 352 
-LINE.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 10, 30 

Specify Planet view factor information 
-CURSET 35, 65 
side 1 (-Z) 

-CURSET 35, 105 
side 2 (-Y) 

-CURSET 35, 145 
side 3 (-X) 

-CURSET 35, 185 
side 4 (+Y) 

-CURSET 35, 225 
side 5 (+X) 

-CURSET 35, 265 
side 6 (+Z) 

-EDIT ”[FE1U1", 180, 40, 100, 35 -L 
[[FE1]] 

-LISTEND 

-EDIT ”[FE2U] M , 180, 80, 100, 35 -L 
[[FE2]] 

-LISTEND 

-EDIT "[FE3U]", 180, 120, 100, 35 -L 
[[FE3]] 


-LISTEND 

-EDIT "[FE4U]", 180, 160, 100, 35 -L 
[[FE4]] 

-LISTEND 

-EDIT "[FE5U]", 180, 200, 100, 35 -L 
[[FE5]] 

-LISTEND 

-EDIT ”[FE6U]”, 180, 240, 100, 35 -L 
[[FE6]] 

-LISTEND 

-BUTTON "OK", 140, 280, 75, 35 -I 

-END 


Screens 
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Professional 


Specify in global coordinates 
Solor vector Planet vector 


xs 

-1 

XE 

0 

YS 

0 ] 

YE 

0 

ZS 

0 

ZE 

b 1 


1 1 OK 


Figure 42 ~[GSU 


~[GS1] 

-SET.WIN 88, 0, 324, 248 
-LINE.COLOR BLACK 
~FILL_COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 44, 23 

Specify in global coordinates 

-CURSET 18, 53 

Solar Vector 

-CURSET 20, 90 

XS 

-CURSET 20, 130 
YS 

-CURSET 20, 170 

zs 

-EDIT ”[GS1]", 75, 65, 49, 35 -M:l 

-EDIT "[GS2]", 75, 105, 49, 35 -M:l 

-EDIT "[GS3]”, 75, 145, 49, 35 -M:l 

-CURSET 178, 50 

Planet Vector 

-CURSET 175, 90 

XE 

-CURSET 175, 130 
YE 

-CURSET 175, 170 
ZE 

-EDIT "(GEIF, 230, 65, 49, 35 -M:l 
-EDIT "[GE2]”, 230, 105, 49, 35 -Mil 
-EDIT "fGE3r, 230, 145, 49, 35 -Mil 
-BUTTON "OK”, 116, 201, 75, 35 -I 

-EDIT.FILL GSSCR.DAT 
-END 
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EX SYS Professional!! 

Specify orbit information 

inclination (deg) 
nodal crossing (deg) 
argument of periapsis (deg) 
true anomaly at epoch (deg) 
semi-major axis (km) 
eccentricity 




Figure 43 -[INC] 


-[INC] 

-SET.WIN 88, 0, 390, 351 
-BACKGROUND CYAN 
-CLS 

-CURSET 91, 32 
Specify orbit information 
-CURSET 14, 75 
inclination (deg) 

-CURSET 14, 115 
nodal crossing (deg) 

-CURSET 14, 155 
argument of periapsis (deg) 

-CURSET 14, 195 

true anomaly at epoch (deg) 

-CURSET 14, 235 
semi-major axis (km) 

-CURSET 14, 275 
eccentricity 

-EDIT ’’[INC]”, 280, 50, 100, 35 -M:l 
-EDIT ’’[OMEGA]’’, 280, 90, 100, 35 -M:l 
-EDIT "[WOMEGA]’’, 280, 130, 100, 35 -M:l 
-EDIT ’’[v0]’’, 280, 170, 100, 35 -M:l 
-EDIT ’’[SMA]”, 280, 210, 100, 35 -M:l 
-EDIT ’’[ECCr, 280, 250, 100, 35 -M:l 
-BUTTON ”OK”, 140, 302, 75, 35 -I 

-EDIT.FILL INCSCR.DAT 
-END 


Screens 
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. EXSYS Rjcofessibnal 

Insulated Sides Information 

Is; insulated side number 
2nd insulated side number 
3rd insulated side number 
4th insulated side number 
5th insulated side number 








mm 

Figure 44 ~[INS1] 


— [INS1] 

-SET.WIN 88, 0, 352, 310 
-LINE.COLOR BLACK 
-FILL.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 64, 25 

Insulated Sides Information 

-CURSET 10, 65 

1st insulated side number 

-CURSET 10, 105 

2nd insulated side number 

-CURSET 10, 145 

3rd insulated side number 

-CURSET 10, 185 

4th insulated side number 

-CURSET 10, 225 

5th insulated side number 

-EDIT "[INSIF, 261, 39, 60, 33 -M:l 

-EDIT "[INS2]", 261, 80, 60, 33 -M:l 

-EDIT "[INS3]", 261, 120, 60, 33 -M:l 

-EDIT "[INS4f, 261, 160, 60, 33 -M:l 

-EDIT "[INS5]”, 261, 200, 60, 33 -M:l 

-BUTTON "OK", 125, 255, 75, 33 -I 

-END 
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. ; EXSYS Professional 

Choose the surface material: 


| : btack' 


|. while paint 




[Itlopl'cai , solar ,ref lector : rm 


| grophtle/ epoxy 



Oluminjzeef kapton 



[ • ■ ; . ; OtomiflU m ' [ 

olumtoum lope • 

1 

| : solar celts ~ j 


1 

.gold- ;• : . y 


lisllijRaju j t;^oy 


Figure 45 -[MATERIAL] 


-[MATERIAL] 

-SETJVIN 88, 0, 343, 440 
-LINE.COLOR BLACK 
-FILL.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 42, 23 
Choose the surface material: 

-BUTTON "black paint", "[MATERIAL] black paint”', 25, 30, 285, 25 -I 
-BUTTON "white paint", "[MATERIAL] 'white paint’", 25, 60, 285, 25 -I 
-BUTTON "optical solar reflector", "[MATERIAL] 'optical solar reflector"’, 25, 90, 285, 25 -I 
-BUTTON "graphite/epoxy”, "[MATERIAL] 'graphite/epoxy'", 25, 120, 285, 25 -I 
-BUTTON "aluminized kapton", "[MATERIAL] 'aluminized kapton"’, 25, 150, 285, 25 -I 
-BUTTON "tiodized titanium", "[MATERIAL] 'tiodized titanium”', 25, 180, 285, 25 -I 
-BUTTON "aluminum", "[MATERIAL] 'aluminum'”, 25, 210, 285, 25 -I 
-BUTTON "aluminum tape", "[MATERIAL] aluminum tape'”, 25, 240, 285, 25 -I 
-BUTTON "deposited aluminum”, "[MATERIAL] 'deposited aluminum"’, 25, 270, 285, 25 -I 
-BUTTON "anodized aluminum”, "[MATERIAL] 'anodized aluminum’", 25, 300, 285, 25 -I 
-BUTTON "solar cells", "[MATERIAL] 'solar cells'", 25, 330, 285, 25 -I 
-BUTTON "gold”, "[MATERIAL] 'gold'", 25, 360, 285, 25 -I 

-BUTTON "multilayered insulation”, ’’[MATERIAL] 'multilayered insulation'”, 25, 390, 285, 25 -I 
-END 
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EXSYS Professional 

Enter number of orbit point? 

to be c 

§1111 

alculot'ed: 

12 



Figure 46 -[NUMPTS] 


-[NUMPTS] 

-SET_WIN 88, 0, 500, 120 
-LINE.COLOR BLACK 
-FILL.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 13, 42 

Enter number of orbit points to be calculated: 
-EDIT ’’[NUMPTS]", 415, 20, 75, 35 -M:l 
-BUTTON ’’OK’, 201, 64, 75, 35 -I 

-EDIT.FILL PTSSCR.DAT 
-END 
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EX SYS Professional X 

Choose the Planet: 

1 

‘ MERCURY 



,/ VENUS. 

£ARTH ; :;:. Il 




f ill MRS" X v Ef : j 


Si.'iJUBTER' :i;:| 

1 
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[ 
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Figure 47 -[PLANET] 


-[PLANET] 

-SET.WIN 88, 0, 254, 345 
-LINE.COLOR BLACK 
-FILL.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 38, 23 
Choose the planet: 

-BUTTON "MERCURY’, ’’[PLANET] ’MERCURY ”, 25, 30, 192, 25 -I 
-BUTTON "VENUS’’, ’’[PLANET] VENUS’", 25, 60, 192, 25 -I 
-BUTTON "EARTH”, ’’[PLANET] ’EARTH’", 25, 90, 192, 25 -I 
-BUTTON "MOON ”, ’’[PLANET] 'MOON'”, 25, 120, 192, 25 -I 
-BUTTON "MARS ”, ’’[PLANET] ’MARS’”, 25, 150, 192, 25 -I 
-BUTTON "JUPITER”, "[PLANET] ’JUPITER”’, 25, 180, 192, 25 -I 
-BUTTON "SATURN”, ’’[PLANET] ’SATURN"’, 25, 210, 192, 25 -I 
-BUTTON "URANUS", ’’[PLANET] ’URANUS’”, 25, 240, 192, 25 -I 
-BUTTON "NEPTUNE”, ’’[PLANET] 'NEPTUNE'", 25, 270, 192, 25 -I 
-BUTTON ’’PLUTO", ’’[PLANET] ’PLUTO’", 25, 300, 192, 25 -I 
-END 
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g EXSYS Professional }|||| 

Enter the internal heat load to dissipate (W) 

1.0 
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Figure 48 -[QPOWER] 


-[QPOWER] 

-SET.WIN 88, 0, 575, 115 
~LINE_COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 29, 38 

Enter the internal heat load to dissipate (W) 
-EDIT "[QPOWER]", 460, 15, 103, 35 -M:l 
-BUTTON ”OK', 247, 60, 75, 35 -I 

-EDIT.FILL QPSCR.DAT 
-END 
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• EXSYS Professsonol ■ . . ' ■ 

Enter temperature of radiating surface (K) 

273 



POKi 





Figure 49 Surface Temperature screen 


-[TEMP] 

-SET.WIN 88, 0, 575, 105 
-LINE.COLOR BLACK 
-FILL.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 29, 39 

Enter the temperature of radiating surface (K) 
-EDIT "[TEMP]", 450, 16, 103, 35 -M:l 
-BUTTON "OK", "", 242, 60, 75, 35 -I 
-EDIT.FILL TEMPSCR.DAT 
-END 
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-[X] 

~SET_WIN 0, 0, 608, 685 
-LINE.COLOR BLACK 
-FILL "COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-TEXTSIZE LARGE 
-CURSET 160, 35 
Enclosure Properties 
-TEXTSIZE NORMAL 
-CURSET 6, 60 
Enclosure dimensions: 
-CURSET 21, 1 00 
X: 

-CURSET 21, 140 


Y: 

-CURSET 21, 180 
Z: 

-CURSET 142, 100 
(m) 

-CURSET 142, 140 
(m) 

-CURSET 142, 180 
(m) 

-EDIT ”[X]”, 57, 75, 75, 35 -S:S -M:l 

-EDIT "[Y]'\ 57, 115, 75, 35 -S:S -M:l 

-EDIT ”[Zr, 57, 155, 75, 35 -S:S -M:l 

-CURSET 288, 60 

Solar absoprtivity 

-CURSET 206, 100 

side 1 (-Z) 
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-CURSET 206, 140 
side 2 (-Y) 

-CURSET 206, 180 
side 3 (-X) 

-CURSET 206, 220 
side 4 (+Y) 

-CURSET 206, 260 
side 5 (+X) 

-CURSET 206, 300 
side 6 (+Z) 

-EDIT "[ALPSir, 325, 75, 100, 35 -M:l 
-EDIT "[ALPS2]", 325, 115, 100, 35 -M:l 
-EDIT "[ALPS3]”, 325, 155, 100, 35 -M:l 
-EDIT "[ALPS4]", 325, 195, 100, 35 -M:l 
-EDIT "[ALPS5]", 325, 235, 100, 35 -M:l 
-EDIT "[ALPS6]’’, 325, 275, 100, 35 -M:l 
-CURSET 477, 60 
Emissivity 

-EDIT "[ALPIRl]", 473, 70, 100, 36 -M:l 
-EDIT "[ALPIR2]”, 473, 115, 100, 35 -M:l 
-EDIT "[ALPIR3]", 473, 155, 100, 35 -M:l 
-EDIT "[ALPIR4]”, 473, 195, 100, 35 -M:l 
-EDIT "[ALPIR5]", 473, 235, 100, 35 -M:l 
-EDIT "[ALPIR6]”, 473, 275, 100, 35 -M:l 
-CURSET 2, 395 
Number of decoupled sides: 

-TEXTSIZE SMALL 
-CURSET 5, 427 

a decoupled side is one with no net heat across 
it 

-EDIT "[NUMINS]”, 273, 359, 47, 35 -M:l 
-TEXTSIZE NORMAL 
-IMAGE "SSC.PCX”, 154, 478, 318, 196 
-CURSET 210, 470 
Spacecraft Coordinates 
-LINE.COLOR WHITE 
-FILL.COLOR BLACK 
-TEXTSIZE LARGE 
-FOREGROUND WHITE 
-BACKGROUND BLACK 
-CURSET 211, 347 
Internal Surfaces are Black 
-BUTTON "OK”, 402, 386, 75, 35 -I 

-TEXTSIZE NORMAL 
-EDIT.FILL XSCR.DAT 
-END 
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EXSYS Professional 


Specify in local coordinates 
Solar vector Planet vector 


XS 

-1 

XE 

1° 

YS 

0 

YE 

0 

ZS 

0 

ZE 

0 
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Figure 51 ~[XS] 


-[XS] 

-SET.WIN 88, 0, 325, 244 
-LINE.COLOR BLACK 
-FILL.COLOR BLACK 
-FOREGROUND BLACK 
-BACKGROUND CYAN 
-CLS 

-CURSET 44, 23 
Specify in local coordinates 
-CURSET 18, 53 
Solar Vector 
-CURSET 20, 90 
XS 

-CURSET 20, 130 
YS 

-CURSET 20, 170 
ZS 

-EDIT ”[XSF, 75, 65, 49, 35 -M:l 

-EDIT "[YS]", 75, 105, 49, 35 -M:l 

-EDIT "[ZS]”, 75, 145, 49, 35 -M:l 

-CURSET 178, 50 

Planet Vector 

-CURSET 175, 90 

XE 

-CURSET 175, 130 
YE 

-CURSET 175, 170 
ZE 

-EDIT "[XE]", 230, 65, 49, 35 -M:l 
-EDIT "[YE]”, 230, 105, 49, 35 -M:l 
-EDIT "[ZE]”, 230, 145, 49, 35 -M:l 
-BUTTON ”OK", 122, 196, 75, 35 -I 

-EDIT.FILL XSSCR.DAT 
-END 


C32 


Screens 




Surface Analysis Results 


Surface Area:..!:.....:. .1286 m-2 
Surface Temperture:... 273 K !; 


Rle Output 

select one value 

□ write a new file 
a append to file 

□ No file output 


Radiation effects to include sun AND planet AND albedo 


isdrfdc^^ i 

surface absorptiyity:!:|||65 : 

Solor Flux: 1 356 W/m-2 

jPtddef Flux;:S:;236 W/m-2 
Solar view factor:.. v:: ; ; 0 
Rfanet view factor:..... .248 
lAl t>e<d6! ic oeWicie : !1 


Heat load radiated from surface:... 2.02 W. 
Heat load due to the sun:.......'..!.... 0.00 W ; 

Heat. toad, due to 

Heat ldad due to tfte" otbe'do:..:!:.! !'|0!'6'4' W : 
Heat Idpd due to tHe!enciosur.e:3^0.G0 W 
Heat load to dissipate:;.;;;:!:'..!:. 4.00 W 


Continue 


Figure 52 - [REPORT 1] 


-[REPORT1] 

-SET.WIN 0, 0, 849, 451 
-FILL_COLOR BLACK 
-BACKGROUND LTGRAY 
-CLS 

-TEXTSIZE LARGE 
-CURSET 204, 42 
Surface Analysis Results 
-TEXTSIZE NORMAL 
-CURSET 9, 89 

Surface Area: [[area]8.4] m A 2 

Surface Temperature:... [[temp]7.2] K 
-CURSET 59, 180 
[[*Q 2]] 

-CURSET 350, 220 

Heat load radiated from surface: [[qrad]5.2] W 

Heat load due to the Sun: [[qsun]5.2] W 

Heat load due to the Planet: [[qplanet]5.2] W 

Heat load due to the Planet's albedo:.... [[qalbedo]5.2] W 

Heat load due to the enclosure: [[qencl]5.2] W 

Heat load to the dissipate: [[qpower]5.2] W 

-CURSET 7, 220 

surface emissivity: [[epsilon]] 

surface absorptivity:... [[alpha]] 

Solar Flux: [[S]] W/m A 2 

Planet Flux: [[E]] W/m A 2 

Solar view factor: [[mu]] 

Planet view factor: [[fell 

Albedo coefficient: [[alb]] 

-FILL.COLOR WHITE 
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-BACKGROUND WHITE 
-RECT 610, 5, 200, 155 
-CURSET 650, 26 
File Output 
-TEXTSIZE SMALL 
-CURSET 645, 48 
Select one value 

-BUTTON "write a new file", "ql9 2", 635, 60, 148, 30 -C 
-BUTTON "append to file", "ql9 1", 635, 90, 150, 30 -C 
-BUTTON "no file output", ”ql9 3", 635, 120, 150, 30 -C 
-BUTTON "Continue", "[reportl] = 1", 365, 393, 100, 35 -I 
-END 
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INVALID Surface Results 


Surfoce Area: -0.1301 m-2 


Surface Temperture:... 100 K 


File Output 

select one value 

□ write a new file 

□ append to file 

□ No file output 


Radiation effects to -include sun AND planet AND albedo 


Surface emissivity:..; 05 

Surf oce obsorptivity:..;!/ .05 
Sotor Flux;.,.,. 1356 W/m-2'- 
|; ; R janet =. 

So ta r vie wf f dc tor: ... . f • G .1* 

Pi onet view factor:..... .248 -. ■ 
Albedo coefficient;. 3 111 


Heat toad radiated from surface:... -0.04 W 

Heot load due to the sun; 0.00 W 

Heat toad d ue'to it be p ionet : . v . X . . - 0 .38. § 
Heat lood due to the albedo:. ....... .-0.65 W 

Heat food ;d|e;foi:||e:' enc ibsure; 1.1 . 0Mm% 
.Heat lobd 'tp dissipate;....! t.QQ W-. 


agbntrnuie 


Figure 53 -[REPORT2] 


-[REPORT2] 

-SET.WIN 0, 0, 850, 448 
-FILL.COLOR BLACK 
-BACKGROUND LTGRAY 
-CLS 

-FILL.COLOR RED 
-TEXTSIZE LARGE 
-BACKGROUND RED 
-CURSET 189, 44 
INVALID Surface Results 
-FILL.COLOR BLACK 
-TEXTSIZE NORMAL 
-CURSET 9, 89 

Surface Area: [[area]] m A 2 

-BACKGROUND LTGRAY 
Surface Temperature:... [[temp]] K 
-CURSET 59, 187 
[[*Q 2]] 

-CURSET 374, 220 

Heat load radiated from surface: [[qrad]5.2] W 

Heat load due to the Sun: [[qsun]5.2] W 

Heat load due to the Planet: [[qplanet]5.2] W 

Heat load due to the Planet's albedo:.... [[qalbedo]5.2] W 

Heat load due to the enclosure: [[qencl]5.2] W 

Heat load to the dissipate: [[qpower]5.2] W 

-CURSET 7, 220 

surface emissivity:!.... [[epsilon]] 

surface absorptivity:... [[alpha]] 

Solar Flux: [[S]] W/m A 2 

Planet Flux: [[E]] W/m A 2 
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Solar view factor: [[mu]] 

Planet view factor: [[fe]] 

Albedo coefficient: [[alb]] 

-FILL.COLOR WHITE 
-BACKGROUND WHITE 
-RECT 609, 14, 200, 155 
-CURSET 658, 35 
File Output 
-TEXTSIZE SMALL 
-CURSET 652, 57 
Select one value 

-BUTTON "write a new file", "ql9 2”, 635, 60, 148, 30 -C 
-BUTTON "append to file”, "ql9 1", 635, 90, 150, 30 -C 
-BUTTON "no file output”, "ql9 3", 635, 120, 150, 30 -C 
-BUTTON "Continue", "[report2] = 1”, 311, 391, 100, 35 -I 
-END 
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Temperature Results =! 

File Output 

select one value 


Surfcce Temperture:... 230.27 K 



□ write a new file 

□ append to file 

□ No file output 


Radiation effects to include sun AMD Planet AMD albedo . 


Surface emissivity:...... .05 

: Surf oce absorptivity : 05 

Solar Flux:..... 1 356 W/rrv-2 IH4 
Planet Flux:T.. 236. W/m^2 t 

■ Solar view factor:..... 0 v s 

Heat flux due to the sun:.,.. 0.00 W/m-2 

Heat flux due to. the planet:. 2.92 W/m~2 

Heat flux due to the albedo:..... 5.04 W/m-2 

Albedo coefficient:..:.. .3 i 

y GcnimuB 









Figure 54 -[REP0RT3] 


-[REP0RT3] 

-SET.WIN 0, 0, 868, 454 
-FILL.COLOR BLACK 
-BACKGROUND LTGRAY 
-CLS 

-TEXTSIZE LARGE 
-CURSET 203, 51 
Temperature Results 
-TEXTSIZE NORMAL 
-CURSET 8, 90 

Surface Temperature: [[temp]] K 
-CURSET 58, 129 
[[*q 2]] 

-CURSET 3, 210 

Surface emissivity:.... [[epsilon]] 

Surface absorptivity:., [[alpha]] 

Solar Flux: [[S]] W/m A 2 

Planet Flux: [[E]] W/m A 2 

Solar view factor: [[mu]] 

Planet view factor: [[fe]] 

Albedo Coefficient:.... [[alb]] 

-CURSET 368, 210 

Heat flux due to the Sun: [[xsun]5.2] W/m A 2 

Heat flux due to the Planet: [[xplanet]5.2] W/m A 2 

Heat flux due to the Planet's albedo:... [[xalbedo]5.2] W/m A 2 

-FILL.COLOR WHITE 

-BACKGROUND WHITE 

-RECT 609, 14, 200, 155 

-CURSET 658, 35 

File Output 

-TEXTSIZE SMALL 


Screens 
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-CURSET 654, 55 
Select one value 

-BUTTON "write a new file”, "ql9 2", 635, 60, 150, 30 -C 
-BUTTON "append to file", ”ql9 1", 635, 90, 150, 30 -C 
-BUTTON "No file output", ”ql9 3", 635, 120, 150, 30 -C 
-BUTTON "Continue", ”[report3] = 1”, 330, 376, 100, 35 -I 
-END 
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File Output 
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□ write a new file 
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□ append to file 
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Figure 55 -[REPORT4] 


-[REPORT4] 

~SET_WIN 0, 0, 1006, 550 
-FILL.COLOR BLACK 
-FOREGROUND WHITE 
-BACKGROUND LTGRAY 
-CLS 

-TEXTSIZE LARGE 
-FOREGROUND BLACK 
-CURSET 270, 39 
Enclosure Analysis Results 
-TEXTSIZE NORMAL 
-CURSET 19, 80 
Surface Temperatures 
-CURSET 39, 100 
side 1: [[Tl]6.2] K 
side 2: [[T2J6.2] K 
side 3: [[T3]6.2] K 
side 4: [[T4J6.2] K 
side 5: [[T5J6.2J K 
side 6: [[T6J6.2] K 
-CURSET 279, 80 
Net Surface Heat loads 
-CURSET 299, 100 
side 1: [[Q116.2] W 
side 2: [[Q2]6.2] W 
side 3: [[Q316.2] W 
side 4: [[Q4]6.2] W 
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side 5: [[Q5J6.2] W 
side 6: [[Q6]6.2] W 
-CURSET 85, 325 
External Surface Properties 
-CURSET 55, 350 
emissivity 
-CURSET 25, 385 
side 1: [[ALPIR1]] 
side 2: [[ALPIR2]] 
side 3: [[ALPIR3]] 
side 4: [[ALPIR4]] 
sideS: [[ALPIR5]] 
side 6: [[ALPIR6]] 

-CURSET 244, 350 
absorptivity 
-CURSET 216, 385 
side 1: [[ALPS1]] 
side 2: [[ALPS2]] 
side 3: [[ALPS3]] 
side 4: [[ALPS4]] 
side 5: [[ALPS5]] 
side 6: [[ALPS6]] 

-CURSET 610, 80 

Orbit Information 

-CURSET 559, 100 

Inclination: [[INC]] deg 

Longitude of Ascending node: [[OMEGA]] deg 

Argument of Periapsis: [[WOMEGA]] deg 

True Anomaly at Epoch: [[vO]] deg 

Semi-mejor axis: [[SMA]] km 

Eccentricity: [[ECC]] 

[[*q 5]] 

[[*q 4]] 

-CURSET 483, 326 
View Factors 
-CURSET 432, 352 
Solar 

-CURSET 387, 385 
side 1: [[MU1]] 
side 2: [[MU2]] 
side 3: [[MU3]] 
side 4: [[MU4]] 
side 5: [[MU5]] 
side 6: [[MU6]] 

-CURSET 593, 353 
Planet 

-CURSET 551, 384 
side 1: [[FE1J6.2] 
side 2: [[FE2J6.2] 
side 3: [[FE3]6.2] 
side 4: [[FE4]6.2] 
side 5: [[FE5J6.2] 
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side 6: [[FE6J6.2] 

~FILL_COLOR WHITE 
-BACKGROUND WHITE 
-RECT 791, 379, 200, 155 
-CURSET 842, 402 
File Output 
-TEXTSIZE SMALL 
-CURSET 838, 421 
Select one value 

-BUTTON "write a new file", "ql9 2", 816, 433, 150, 30 -C 
-BUTTON "append to file", "ql9 1", 816, 463, 150, 30 -C 
-BUTTON "No file output”, "ql9 3", 816, 493, 150, 30 -C 
-BUTTON "Continue", "[REPORT4] = 1", 828, 308, 106, 35 -I 
-END 
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Orbit Analysis Results 


Orbit average temperature: 193.75 K 

Higheil temperature delta during orbit; t5;0T K;;;true gnomaiy Cl dea 
Highest temperature during orbit: 196.81 K side: 5 true anomaly 0 deg 

Lowest temperature during or bit:::19i:80: 1 side; . 3 true dndmbly 0. deg 


Orbit overage heat load: 16.14 

Wight est: he at load delta . dufing. orbit:: 6175: W; 


side- 
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mm 
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: sicfe --OSinHiiiSJde ' 

side 6; : ; : ;;05 : ;;;;i : side .6:1 ;...05 : ; 


Continue" 


true anomaly 0 deg . I:!::::ll:l; : ;ll!iiii0!i 

HilH 1 1 Oilli 0 tNl t^^QE mot t q i|if ili : i|l|i|! : y ! : i Hill 
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i i i S oiT*!?* Si ! ! I ! i ^ • M ! " III ^ I ! ^ 

Eccentricity: 0 ... . 1; 

The spacecraft ofientatioo is planet " 


File Output 
select one value 

□ w ri te a new file 

□ append to file 

□ No file output 


Figure 56 -[REPORT5] 


-[REPORT5] 

~SET_WIN 0, 0, 1006, 705 
-FOREGROUND WHITE 
-BACKGROUND LTGRAY 
-CLS 

-TEXTSIZE LARGE 
-FOREGROUND BLACK 
-CURSET 283, 44 
Orbit Analysis Results 
-TEXTSIZE NORMAL 
-CURSET 42, 84 

Orbit average temperature: [[tavet]6.2] K 

Highest temperature delta during orbit: [[tdelta]6.2] K : true anomoly [[vOdelt]] deg 
Highest temperature during orbit: [[thigh]6.2] K side: [[shigh]] true anomoly: [[vOhigh]] deg 
Lowest temperature during orbit: [[tlow]6.2] K side: [[slow]] true anomoly: [[vOlow]] deg 
-CURSET 42, 200 

Orbit average heat load: [[qavet]6.2] W 

Highest heat load delta during orbit: [[qdelta]6.2] W : true anomoly [[vOqdel]] deg 

-CURSET 10, 270 

Orbit Average Temps 

-CURSET 30, 300 

side 1: [[Tl]6.2] K 

side 2: [[T2]6.2] K 
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side 3: [[T3]6.2] K 
side 4: [[T4]6.2] K 
side 5: [[T5J6.2] K 
side 6: [[T6]6.2] K 
-CURSET 254, 270 
Orbit Average Heat loads 
-CURSET 262, 301 
side 1: [[Ql]6.2] W 
side 2: [[Q216.2] W 
side 3: [[Q316.2] W 
side 4: [[Q4]6.2] W 
side 5: [[Q5]6.2] W 
side 6: [[Q6]6.2] W 
-CURSET 102, 472 
External Surface Properties 
-CURSET 51, 497 
emissivity 
-CURSET 33, 530 
side 1: [[ALPIR1]] 
side 2: [[ALPIR2]] 
side 3: [[ALPIR3]] 
side 4: [[ALPIR4]] 
side 5: [[ALPIR5]] 
side 6: [[ALPIR6]] 

-CURSET 296, 497 
absorptivity 
-CURSET 265, 530 
side 1: [[ALPS1]] 
side 2: [[ALPS2]] 
side 3: [[ALPS3]] 
side 4: [[ALPS4]] 
side 5: [[ALPS5]] 
side 6: [[ALPS6]] 

-CURSET 558, 270 
Orbit Information 
-CURSET 497, 300 

Number of orbit points analyzed: [[numpts]] 

Inclination: [[INC]] deg 

Longitude of Ascending node: [[OMEGA]] deg 
Argument of Periapsis: [[WOMEGA]] deg 
True Anomaly at Epoch: [[v0]] deg 
Semi-major axis: [[SMA]] km 
Eccentricity: [[ECC]] 

[[*Q4]] 

-BACKGROUND WHITE 
-RECT 721, 521, 200, 155 
-CURSET 773, 547 
File Output 
-TEXTSIZE SMALL 
-CURSET 770, 566 
Select one value 

-BUTTON "write a new file”, ”ql9 2”, 750, 570, 150, 30 -C 
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-BUTTON ’’append to file”, ”ql 9 1”, 750, 600, 150, 30 -C 
-BUTTON "No file output", "ql9 3”, 750, 630, 150, 30 -C 
-BUTTON "Continue”, ”[report5] = 1”, 522, 557, 106, 35 -I 
-END 
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Su rf.qce: :; iAirva lysiSjBes u I ts 

The surface is used for inferior structure 
The. surface is msde # btock^point f:; 

the perforr^bn4:^!^j5^e;:iis= BOtT 

Emissivitv .9 . ' ; ;i 
Spier obsorptivtfy .9 • 


■ 


iiCohttnue 


file Output 
select one value 

□ write a new file 

□ append to file 
O No file output 


Figure 57 -[REPORT6] 


-[REPORT6] 

-SET.WIN 0, 0, 651, 292 
-BACKGROUND LTGRAY 
-CLS 

-CURSET 160, 29 
Surface Analysis Results 
-CURSET 50, 65 

The surface is used for [[APPLICATION]] 

The surface is made of [[MATERIAL]] 

-CURSET 53, 119 
[[*Q9]] 

-CURSET 50, 190 
Emissivity [[EPSILON]] 

Solar absorptivity [[ALPHA]] 

-BACKGROUND WHITE 
-RECT 413, 138, 200, 144 
-CURSET 466, 159 
File Output 
-TEXTSIZE SMALL 
-CURSET 458, 175 
Select one value 

-BUTTON "write a new file", "q21 2”, 445, 180, 150, 30 -C 
-BUTTON "append to file", "q21 1", 445, 211, 150, 30 -C 
-BUTTON "No file output", "q213", 444, 243, 150, 30 -C 
-BUTTON "Continue”, "[REPORT6] = 1”, 252, 244, 106, 35 -I 
-END 
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Screen Default Data Files 


Q16 

ALPSCR.DAT 

Q17 

ALBSCR.DAT 

[AINST] 

AINSTSCR.DAT 

[AREA] 

AREASCR.DAT 

[ECC] 

ECCSCR.DAT 

[GS1] 

GSSCR.DAT 

[INC] 

INCSCR.DAT 

[NUMPTS] 

PTSSCR.DAT 

[QENCL] 

QENCLSCR.DAT 

[QPOWER] 

QPSCR.DAT 

[TEMP] 

TEMPSCR.DAT 

[X] 

XSCR.DAT 

[XS] 

XSSCR.DAT 


Default Data Files 
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AINSTSCR.DAT 

1 /* Surface area of the enclosure 

instrument 

1 /*emissivity of instrument surface 

inside the enclosure 


ALBSCR.DAT 
.3 /* albedo coefficient 

1356 /* solar flux W/m A 2 

236 /* planet flux W/m A 2 

6378 /* radius of the planet km 


ALPSCR.DAT 

0 /* solar view factor 

.287 /* planet view factor 

.05 /* emissivity of surface 

.05 /* solar absorptivity of surface 


AREASCR.DAT 
1 /* area of surface m A 2 


ECCSCR.DAT 

0 /* true anomaly at epoch deg 

6785.0 /* semi-msuor axis km 
0 /* eccentricity 


GSSCR.DAT 

-1 /* x component of the global sun 

vector 

0 /* y component of the global sun 

vector 

0 /* z component of the global sun 

vector 


INCSCR.DAT 
28.5 /* inclination deg 

0 /* line of nodes deg 

0 /* line of periapsis deg 

0 f* true anomaly at epoch deg 

6785.0 /* semi-msgor axis km 

0 /* eccentricity 


PTSSCR.DAT 

12 /* number of points in the orbit 

QENCLSCR.DAT 

0 /* enclosure radiation to dissipate 

QPSCR.DAT 

1.0 /* internal power to dissipate W 


TEMPSCR.DAT 

273 f* temperature of the surface 


XSCR.DAT 

1.0 /* dimension along the x axis 

1.0 /* dimension along the y axis 

1.0 f* dimension along the z axis 

.05 /* solar absorptivity for each surface 

.05 

.05 

.05 

.05 

.05 

.05 /* emissivity for each surface 

.05 

.05 

.05 

.05 

.05 

0 /* number of decoupled sides 


XSSCR.DAT 

-1 /* x component of the local sun vector 

0 /* y component of the local sun vector 

0 f* z component of the local sun vector 

0 /* x component of the local planet 

vector 

0 /* y component of the local planet 

vector 

0 /* z component of the local planet 

vector 
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Report Generator 



Q19 3 /G:END 

Q1 1 /GrREPORTl 
Q1 2 /G:REPORT2 
Q1 3 /G:REPORT3 
Q1 4 /G:REPORT4 
Q1 5 /G:REPORT5 

:REPORTl 

Q19 1 /GAPAREA 

Q19 2 /G:NWAREA 

APAREA 

FILE AREA.RPT /A 
GOTO RPT1 

: NWAREA 
FILE AREA.RPT 
GOTO RPT1 

•JIPT1 

"Surface Analysis Results" 
[AREA] /F6.2 
[TEMP] /F6.2 
Q2 

[QRAD] /F6.2 
[QSUN] /F6.2 
[QPLANET] /F6.2 
[QALBEDO] /F6.2 
[QENCL] /F6.2 
[QPOWER] /F6.2 
[EPSILON] 

[ALPHA] 

[S] 

[E] 

[MU] 

[FE] 

[ALB] 

GOTO END 

JIEPORT2 

Q19 1 /GAPTEMP 

Q19 2 /GtNWTEMP 

APTEMP 

FILE TEMP.RPT/A 
GOTO RPT2 

: NWTEMP 
FILE TEMP.RPT 
GOTO RPT2 

Report Generator 


JIPT2 

Q12 1 /G:SURTEMP 
Q12 2 /G:ENVTEMP 
:SURTEMP 

"Surface Analysis Results" 
[AREA] /F6.2 
[TEMP] /F6.2 
Q2 

[QRAD] /F6.2 
[QSUN] /F6.2 
[QPLANET] /F6.2 
[QALBEDO] /F6.2 
[QENCL] /F6.2 
[QPOWER] /F6.2 
[EPSILON] 

[ALPHA] 

[S] 

[E] 

[MU] 

[FE] 

[ALB] 

GOTO END 
:ENVTEMP 

"Temperature Analysis Results" 
[TEMP] /F6.2 

Q 2 

[EPSILON] 

[ALPHA] 

[S] /F6.2 
[E] /F6.2 
[MU] 

[FE] 

[ALB] 

[XSUN] /F6.2 
[XPLANET] /F6.2 
[XALBEDO] /F6.2 

GOTO END 

:REPORT3 
Q19 1 /GAPENCL 
Q19 2 /G:NWENCL 

APENCL 

FILE ENCL.RPT /A 
GOTO RPT3 

: NWENCL 
FILE ENCL.RPT 
GOTO RPT3 
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:RPT3 

"Enclosure Analysis Results" 

"Surface Temperatures" 

[Tl] /fB.2 
[T2] /re. 2 
[T3] /fB.2 
[T4] /fB.2 
[T5] /fB.2 
[T61 /fB.2 

"Net Surface Heat loads" 

[Ql] /fB.2 
[Q2] /ffi.2 
[Q3] /fB.2 
[Q4J /fB.2 
[Q5] /fB.2 
[Q6] /fB.2 

"External Surface Properties" 

"emissivity" 

[ALPIR1] 

[ALPIR2] 

[ALPIR3] 

[ALPIR4] 

[ALPIR5] 

[ALPIR6] 

"absorptivity” 

[ALPS1] 

[ALPS2] 

[ALPS3] 

[ALPS4] 

[ALPS5] 

[ALPS6] 

"Orbit Information" 

[INC] 

[OMEGA] 

[WOMEGA] 

[vO] 

[SMA] 

[ECC] 

q5 

q4 


"View Factors" 

"Solar" 

[MU1] 

[MU2] 

[MU3] 

[MU4] 

[MU5] 

[MU6] 

"Planet" 

[FE1] 

[FE2] 

[FE3] 

[FE4] 

[FE5] 

[FE6] 

GOTO END 

REPORT4 

Q19 1 /GrAPORBIT 

Q19 2 /G:NWORBIT 

APORBIT 

FILE ORBIT.RPT /A 
GOTO RPT4 

: NWORBIT 
FILE ORBIT.RPT 
GOTO RPT4 

:RPT4 

"Orbit Analysis Results" 

[TAVET] /F6.2 
[TDELTA] /F6.2 
[VODELT] 

[THIGH] /F6.2 
[SHIGH] 

[VOHIGH] 

[TLOW] /F6.2 
[SLOW] 

[VOLOW] 

[QAVET] /F6.2 
[QDELTA] /F6.2 
[VOQDEL] 

"Surface Temperatures" 
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[Tl] /f6.2 

[MU2] 

[T2] /fB.2 

[MU3] 

[T3] /fB.2 

[MU4] 

[T4] /fB.2 

[MU5] 

[T5] /fB.2 

[MU6] 

[T6] /f6.2 



"Planet" 

"Net Surface Heat loads" 



[FE1] 

[Ql] /fB.2 

[FE2] 

[Q2] /fB.2 

[FE3] 

[Q3] /fB.2 

[FE4] 

[Q4] /fB.2 

[FE5] 

[Q5] /fB.2 

[FE6] 

[Q6] /fB.2 



"Orbit Average Temps' 

"External Surface Properties" 

"Orbit Information" 


"emissivity” 

[NUMPTS] 


[INC] 

[ALPIR1] 

[OMEGA] 

[ALPIR2] 

[WOMEGA] 

[ALPIR3] 

[vO] 

[ALPIR4] 

[SMA] 

[ALPIR5] 

[ECC] 

[ALPIR6] 

Q4 

"absorptivity” 

GOTO END 

[ALPS1] 

JIEPORT5 

[ALPS2] 

Q19 1 /GAPSURF 

[ALPS3] 

Q19 2 /G:NWSURF 

[ALPS4] 


[ALPS5] 

APSURF 

[ALPS6] 

FILE SURF.RPT /A 


GOTO RPT5 

"Orbit Information" 



: NWSURF 

[INC] 

FILE SURF.RPT 

[OMEGA] 

GOTO RPT5 

[WOMEGA] 


[vO] 

:RPT5 

[SMA] 

"Surface Analysis Results' 

[ECC] 


Q5 

[APPLICATION] 

Q4 

[MATERIAL] 

"View Factors" 

Q9 

"Solar” 

[EPSILON] 


[ALPHA] 

[MU1] 

:END 


Report Generator 
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AREA.RPT 


Surface Analysis Results 
area of radiator (m A 2) = 0.06 
temperature of the radiator (K) = 273.00 
Radiation effects to include internal power only 
heat radiated from the surface (W) = 1.00 

heat load on surface due to sun (W) = 0.00 
heat load on surface due to planet (W) = 0.00 
heat load on the surface due to albedo (W) = 0.00 
Enclosure heat to be radiated (W) = 0.00 
heat load to be dissipated (W) = 1.00 
emissivity of the radiating surface = .05 
absorptivity of the radiating surface = .05 
solar flux W/m A 2 = 1356 
planet flux W/m A 2 = 236 

solar angle (= cos(theta) where theta is the angle (in rad) between the sun 
vector and the plane normal) of the radiating surface. -Between 0 and 1=0 
view factor to the planet of the radiating surface = .248 
albedo coefficient = .3 
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TEMP.RPT 


Temperature Analysis Results 

temperature of the radiator (K) = 238.83 

Radiation effects to include sun AND planet AND albedo 

emissivity of the radiating surface = .05 

absorptivity of the radiating surface = .05 

solar flux W/m A 2 = 1356.00 

planet flux W/m A 2 = 236.00 

solar angle (= cos(theta) where theta is the angle (in rad) between the sun 
vector and the plane normal) of the radiating surface. -Between 0 and 1 = 0 
view factor to the planet of the radiating surface = .287 
albedo coefficient = .3 

heat flux on surface due to sun (W/m A 2) = 0.00 

heat flux on surface due to radiation of planet (W/m A 2) = 3.39 

heat flux on surface due to albedo (W/m A 2) = 5.84 


Report Generator 
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ENCL.RPT 


Enclosure Analysis Results 
Surface Temperatures 
Temperature of side 1 (K) = 295.67 
Temperature of side 2 (K) = 295.67 
Temperature of side 3 (K) = 303.36 
Temperature of side 4 (K) = 295.67 
Temperature of side 5 (K) = 298.25 
Temperature of side 6 (K) = 295.67 
Net Surface Heat loads 
Net heat on side 1 = *12.44 
Net heat on side 2 = -12.44 
Net heat on side 3 = 43.79 
Net heat on side 4 = -12.44 
Net heat on side 5 = 5.97 
Net heat on side 6 = -12.44 
External Surface Properties 
emissivity 

Infrared absorptivity of side 1 (assuming gray bodies = emissivity) = .05 

Infrared absorptivity of side 2 (assuming gray bodies = emissivity) = .05 

Infrared absorptivity of side 3 (assuming gray bodies = emissivity) = .05 

Infrared absorptivity of side 4 (assuming gray bodies = emissivity) = .05 

Infrared absorptivity of side 5 (assuming gray bodies = emissivity) = .05 

Infrared absorptivity of side 6 (assuming gray bodies = emissivity) = .05 

absorptivity 

Solar absorptivity of side 1 = .05 

Solar absorptivity of side 2 = .05 

Solar absorptivity of side 3 = .05 

Solar absorptivity of side 4 = .05 

Solar absorptivity of side 5 = .05 

Solar absorptivity of side 6 = .05 

Orbit Information 

orbit inclination (deg) = 0 

longitude of ascending node (deg) = 0 

satellite position (deg) = 0 

true anomaly at epoch = 0 

Semi-major axis (km) = 6785 

eccentricity of orbit = 0 

The spacecraft orbit is LEO (Low Earth Orbit) 

The spacecraft orientation is inertial 

View Factors 

Solar 

Solar view factor for side 1 = 0 
solar view factor for side 2 = 0 
solar view factor for side 3=1 
solar view factor for side 4 = 0 
Solar view factor for side 5 = 0 
solar view factor for side 6 = 0 
Planet 

Planet view factor for side 1 = .287116 
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planet view factor for side 2 
Planet view factor for side 3 
Planet view factor for side 4 
Planet view factor for side 5 
Planet view factor for side 6 


.287116 

0 

.287116 

.883628 

.287116 
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ORBIT.RPT 


Orbit Analysis Results 

the average temperature in the enclosure over the entire orbit (K) = 

Highest temperature delta (K) = 5.01 

position with the highest temperture delta = -4.712389 

the highest temperature value (K) = 199.96 

the side with the highest temperature value = 4 

The position with the highest temperture value = -5.235988 

Lowest temperature in orbit analysis (K) = 191.80 

side with lowest temperature = 3 

Position of the lowest temperature = 0 

the average heat load in the enclosure over the entire orbit (W) = 

Highest heat load delta (W) = 7.04 

position with the highest heat load delta = -5.759587 

Surface Temperatures 

Temperature of side 1 (K) = 195.59 

Temperature of side 2 (K) = 196.11 

Temperature of side 3 (K) = 196.11 

Temperature of side 4 (K) = 196.11 

Temperature of side 5 (K) = 196.11 

Temperature of side 6 (K) = 195.59 

Net Surface Heat loads 

Net heat on side 1 = 16.30 

Net heat on side 2 = 16.34 

Net heat on side 3 = 16.34 

Net heat on side 4 = 16.34 

Net heat on side 5 = 16.34 

Net heat on side 6 = 16.30 

External Surface Properties 

emissivity 

Infrared absorptivity of side 1 (assuming gray bodies = emissivity) 
Infrared absorptivity of side 2 (assuming gray bodies = emissivity) 
Infrared absorptivity of side 3 (assuming gray bodies = emissivity) 
Infrared absorptivity of side 4 (assuming gray bodies = emissivity) 
Infrared absorptivity of side 5 (assuming gray bodies = emissivity) 
Infrared absorptivity of side 6 (assuming gray bodies = emissivity) 
absorptivity 

Solar absorptivity of side 1 = .05 
Solar absorptivity of side 2 = .05 
Solar absorptivity of side 3 = .05 
Solar absorptivity of side 4 = .05 
Solar absorptivity of side 5 = .05 
Solar absorptivity of side 6 = .05 
Orbit Information 
orbit inclination (deg) = 28.5 
longitude of ascending node (deg) = 0 
satellite position (deg) = 0 
true anomaly at epoch = 0 
Semi-major axis (km) = 6785 
eccentricity of orbit = 0 


195.94 


16.33 


= .05 
= .05 
= .05 
= .05 
- .05 
= .05 
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The spacecraft orientation is inertial 

View Factors 

Solar 

Planet 

Planet view factor for side 1 = .25 
planet view factor for side 2 = .25 
Planet view factor for side 3 = 0 
Planet view factor for side 4 = .25 
Planet view factor for side 5 = .84 
Planet view factor for side 6 = .25 
Orbit Average Temps 
Orbit Information 

Number or orbital points to be analyzed = 12 

orbit inclination (deg) = 28.5 

longitude of ascending node (deg) = 0 

satellite position (deg) = 0 

true anomaly at epoch = 0 

Semi-major axis (km) = 6785 

eccentricity of orbit = 0 

The spacecraft orientation is inertial 
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SURF.RPT 


Surface Analysis Results 

Surface APPLICATION = antenna reflector 

Surface material = white paint 

the performance time is BOL (beginning of life) 

emissivity of the radiating surface = .05 

absorptivity of the radiating surface = .05 
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The following information is organize to create quick references for the designer. 

The first section contains a list of the source code routines and a biref description on each. The second section 
contains a list of fortran files and the routines in each. The third section contians a list of executables and the files 
that must be linked to create them. 

SUBROUTINE AREA 

Calculates the areas of each side. The side axes are defined. 

SUBROUTINE CALCEB 
Calculates the blackbody radiation for each surface 

SUBROUTINE CALCQ 
Calculates the heat load on each surface 

SUBROUTINE CALCT 
Calculates the temperatures given EB 

SUBROUTINE CALCVWF 
Calculates view factors for the internal enclosure 

SUBROUTINE CALFEB 

Calculates the FEB matrix FEB = view factors * areas. It is the coefficient matrix to the blackbody radiation terms, 
EB (see documentation) 

SUBROUTINE CALFEI 

Calculates the earth view factor for inertial orientation 
SUBROUTINE CALFEP 

Calculates the earth view factor for planet orientation 
SUBROUTINE CALFLUX 

Calculates the environmental fluxes. Fills the X matrix and the XRAD matrix 
SUBROUTINE CALSOLR 

Calculates the dot product of the solar vector and the unit normals of each surface 
SUBROUTINE CALVFNOR 

Calculates the view factor from a plane element to a sphere. Used as the earth view factor normal from "thermal 
radiation heat transfer" by Siegel and Howell: APPENDIX C 

SUBROUTINE CALVFTAN 

Calculates the view factor from a plane element to a sphere. Used as the earth view factor tangential from "Thermal 
Radiation Heat Transfer" by Siegel and Howell: APPENDIX C 

FUNCTION DOTPROD 

Calculates the dot product of two vectors x and y 

SUBROUTINE EAREAD 
Reads the input for the enclosure analysis 

SUBROUTINE EAR WRIT 

FORTRAN Program List F3 
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Writes planet view factors for each side, and the normal and tangential components to "RETURN.DAT" 

SUBROUTINE ENCINIT 
Initialize the constants and sets the matrices to 0.0 

SUBROUTINE ENCLREAD 
Reads the input for the enclosure analysis 

SUBROUTINE ENCWRJT 
Writes temperature and heat load to "RETURN.DAT" 

SUBROUTINE REREAD 
Reads the input for the planet view factor analysis 

PROGRAM FETANOR 

Calculates the earth view factors normal and tangential plates. 

SUBROUTINE FEWRIT 

Writes normal and tangential view factors to "RETURN.DAT" 

SUBROUTINE IJK2PQR 

**** THIS SUBROUTINE IS CURRENTLY NOT USED, BUT LEFT FOR FUTURE REFERENCE **** 
Fills the coefficient matrix this conversion was taken from "Fundamentals of Astrodynamics" by Bate, Mueller and 
White 

SUBROUTINE MATADD 

Adds two matrices all three matrices are dimension nxm 
if IFAC =-l then subtraction if IFAC = 1 then addition 

SUBROUTINE MATMUL 

Multiplies two matrices a is dimension N,M B is dimension M,L C is dimension NX 
SUBROUTINE ORBCOMP 

Calculates the orbital average values for temperature and heat loads. 

SUBROUTINE ORB FIN 

Finishes the orbital calculations at the end of the total orbit points. 

SUBROUTINE ORBINIT 

Initialize the constants and sets the matrices to 0.0 For the orbital analysis 
PROGRAM ORBITAL 

This program solves for the temperature and heat loads on an enclosure at distinct orbit points as specified by the 
user. For inertially or planet oriented orbits. 

SUBROUTINE ORBOUT 

Writes temperatures and heat load to "RAD-ORB .DAT" 

SUBROUTINE ORBREAD 
Reads the input for the orbital analysis 

SUBROUTINE ORB WRIT 
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Writes all the final statistical analysis to "RETURN.DAT" 

SUBROUTINE ORBPROPS 

Calculates the orbital radius taken from "Fundamentals of Astrodynamics” by Bate, Mueller and White 
SUBROUTINE POSGLOB 

Calculates the spacecraft position in global coordinates 
SUBROUTINE PQR2IJK 

Fills the coefficient matrix for the rotation from perifocal system (P,Q,R) to the global system (I.J.K) see 
documentation for derivation 

PROGRAM QENCL 

Solves for the temperature and heat loads on an enclosure. 

SUBROUTINE SHADOW 

Determines if the spacecraft is in the shadow of the earth if in the shadows: returns the solar flux as 0.0 And the 
solar view factors as 0.0 If not returns the original value for solar flux 

PROGRAM SOLARV 

Solves for the dot product of the solar vector and the local unit normals on each side of the enclosure. 

SUBROUTINE SOLREAD 
Reads the input for the solar view factor analysis 

SUBROUTINE SOLWRIT 

Writes blackbody radiation, heat load and view factors to "RETURN.DAT" 

SUBROUTINE SUNREAD 
Reads the input for the sun vector analysis 

PROGRAM SUNVECT 

Calculates the sun vector in local coordinates, given inclination, longitude of ascending node and spacecraft position. 
SUBROUTINE SUNWRIT 

Writes the three components of the sun vector in local coordinates and the 6 solar view factors 
SUBROUTINE CALSUN 

Rotates the global sun vector into local coordinates solving the equation ax = b. A is the coefficient matrix from 
perifocal to global, b is the coefficient matrix from perifocal to spacecraft x is the local sun vector g is the global 
sun vector 

SUBROUTINE SURFNOR 

Calculates the dot product of x into the local surface normals as defined in documentation x is a vector in local 
coordinates y is the dot product with the surface normals 

FUNCTION VIEWFAC 

Calculates the view factor for two finite plates at 90 deg taken from "Thermal Radiation Heat Transfer" by Robert 
Siegel and John R. Howell. Appendix C #14 "two finite rectangles of same length, having on common edge, and 
at an angle of 90 deg to each other.” 
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PROGRAM EARTHV 

Calculates the earth view factors at any point in an elliptical solar inertial orbit. Given the position in orbit assuming 
at time 0 the x axis points from nadir, the y axis points in the orbit direction and the z axis completes the right hand 
rule pointing north. 

FUNCTION XMAG 
Calculates the magnitude of the vector x 
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FORTRAN FILE 


ROUTINE NAME 


PROGRAM ORBITAL 
PROGRAM OENCL 

SUBROUTINE ENCLREAO(X.Y.Z.FE.ALPIR.ALPS.EPS.XMU, 
+ ALBEDO, NUMINS.INSUL.S.E) 

SUBROUTINE ENCWRIT(O.T) 

SUBROUTINE AREA(X.Y.Z.A) 

SUBROUTINE CALFEB(F,A,FEB) 



SUBROUTINE CALFLUX(S,XMU,A,ALPS,FE,E,AIR,ALB,EPS,X,XRAD.XRADI) 
SUBROUTINE CALCEB(X,Q,XRADI,EB) 

SUBROUTINE CALCQ(XFLUX,INSUL,NUMINS,FEB,XRADI,Q) 

SUBROUTINE CALCT(EB.SIGMAJ) 

SUBROUTINE CALCVWF(X,Y,Z,A,F) 

SUBROUTINE ENCINIT(SiGMA,XRAD,XRADI,FEB) 

FUNCTION VIEWFAC(HI.'.VI.XL) 

EARTHV.FOR 

PROGRAM EARTHV 

SUBROUTINE EAREAD(W,ECC,SMA,XNU,RE) 
SUBROUTINE EARWRIT(FE.FET.FEN) 

CALVFNT.FOR 

SUBROUTINE CALVFNOR(ECC,SMA,XNU,RE,FEN) 
SUBROUTINE CALVFTAM(ECC,SMA,XNU,RE,FET) 

CALEART.FOR 

SUBROUTINE CALFEP(FEN.FET.FE) 
SUBROUTINE CALFEl(W.FEN.FET.FE) 

SOLARV.FOR 

PROGRAM SOLARV 
SUBROUTINE SOLREAD(X) 
SUBROUTINE SOLWRIT(XMU) 

CALSOLR.FOR 

SUBROUTINE CALSOLR(X.XMU) 
SUBROUTINE SURFNOR(X.Y) 

FORTRAN Program List 
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ORBITAL.FOR 

OENCL.FOR 

OENCLSUB.FOR 



SUNVECT.FOR 

PROGRAM SUNVECT 

SUBROUTINE SUNREAO(GS,XINC, OMEGA, WOMEGA.XNU) 
SUBROUTINE SUNWRIT(X.XMU) 

CALSUN.FOR 

SUBROUTINE CALSUN(G,XINC, OMEGA, WOMEGA.XNU.XS) 
SUBROUTINE POR2IJK(XII,OI.WI,XNU,A,B) 

MATMUL.FOR 

SUBROUTINE MATMUL(N,M,L,A,B,C) 

MATAOD.FOR 

SUBROUTINE MATADD(N,M,A,B,C,IFAC) 

POSGLOB.FOR 

SUBROUTINE P0$GL08(ECC,SMA,XI,0,W,XNU,RVECTG) 

ORBSUBS.FOR 

SUBROUTINE 0RBC0MP(T,0,V0.TAVE,0AVE,THIGH,TL0W,SHIGH,SL0W, 

+ VOHIGH,VCLOW,OHIGH,OLOW,SOHIGH,SOLOW,VOOHIGH,VOOLOW,TDELTA, 
+ ODELTA.VCDELT.VOODEL) 


SUBROUTINE ORBFIN(NUMPTS,TAVET,OAVET,TAVE,OAVE) 

SUBROUTINE ORBINIT(SIGMA,XRAD,XRADI,FEB,TAVE,OAVE.TOELTA, 

+ ODELTA, THIGH, OHIGH.TLOW.OLOW) 

SUBROUTINE ORBREAD(ALBEDO,S,E,GS,Xl t O, WO, ECC.SMA, INERT, X,Y,Z, 

+ RE,ALPIR,ALPS,EPS,NUMINS,INSUL,NUMPTS) 

SUBROUTINE 0RB0UT(0,T,XI,0,W,XN) 

SUBROUTINE ORBWRIT(TAVET,OAVET,TAVE,OAVE,THIGH,TLOW,SHIGH,SLOW, 
+ VOHIGH,VCLOW,OHIGH,OLOW,SOHIGH,SOLOW,VOOHIGH,VOOLOW,TDELTA, 
+ ODELTA.VCDELT.VOODEL) 

ORBPROPS.FOR 

SUBROUTINE ORBPROPS(ECC,SMA,XNU,R) 

OOTPROD.FOR 

FUNCTION DOTPROD(X.Y) 

FETANOR.FOR 

PROGRAM FETANOR 
SUBROUTINE FEREAD(ECC,SMA,XNU) 
SUBROUTINE FEWRIT(FET.FEN) 
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IJK2PQR.F0R 

MAG.FOR 

SHADOW.FOR 
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SUBROUTINE IJK2PQR(XII,0I,WI,XNU,A.B) 

FUNCTION XMAG(X) 

SUBROUTINBHADOW(ECC,XNU,XINC, OMEGA, WOMEGA,GS,SORIG,S,XMU.RE) 
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EXECUTABLE FILE 


FILES TO BE LINKED 


EARTHV.EXE 


FETANOR.EXE 


ORBITAL.EXE 


OENCL.EXE 


SOLARV.EXE 


SUNVECT.EXE 


LINK EARTHV+CALEART+CALVFNT+ORBPROPS 


LINK FETANOR+CALVFNT+ORBPROPS 


L I N K 

ORBITAL+CALEART+CALVFNT+ORBSUBS+SHADOW+CALSUN+QENCLSUBS+AXJ+ 
MATMUL+CALSOLR+DOTPROD+MATADD+ORBPROPS+MAG+POSGLOB 


LINK OENCL+OENCLSUBS+AXJ+MATMUL+MATADD 


LINK SOLARV+CALSOLR+MAG 


LINK SUNVECT+CALSOLR+MAG+CALSUN+AX.B+MATMUL+MATADD 
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QENCL.FOR G6 

QENCLSUB.FOR G9 

EARTHV.FOR G15 

CALVFNT.FOR G17 

CALEART.FOR G18 

SOLARV.FOR G19 

CALSOLR.FOR G21 

SUNVECT.FOR G22 

CALSUN.FOR G24 

MATMUL.FOR G25 

MATADD.FOR . . G25 

POSGLOB.FOR G26 

ORBSUBS.FOR G27 

ORBPROPS.FOR G33 

DOTPROD.FOR G33 

FETANOR.FOR G34 

IJK2PQR.FOR G36 

MAG.FOR G37 

SHADOW.FOR G38 
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ORBIT AL.FOR C 


PROGRAM ORBITAL 
IMPLICIT REAL*8 (A-H. 0-Z) 

C THIS PROGRAM SOLVES FOR THE TEMPERATURE AND 
C HEAT LOAOS ON AN ENCLOSURE AT DISTINCT ORBIT 
C POINTS AS SPECIFIED BY THE USER. FOR INERIIALLY 
C OR PLANET ORIENTED ORBITS. IT IS USED IN 
C CONJUCTION WITH THE EXPERT SYSTEM RADIATE.RUL 
C RACHEL ORMSBY July 20. 1993 
C 

C 0(6) THE HEAT LOAD ON EACH SIDE OF A RECTANGULAR BOX 
C EB(6) THE BLACKBODY RADIATION ON EACH SIDE OF THE BOX 
C F(6,6) THE VIEW FACTORS INSIDE THE BOX 
C FE(6) THE VIEW FACTORS TO EARTH FOR EACH SIDE 
C ALPIR 1HE INFRARED ABSORPTIVITY FOR THE EXTERNAL SURFACES 
C = EMISSIVITY 

C ALPS THE SOLAR ABSORPTIVITY FOR THE EXTERNAL SURFACES 
C EPS THE EMISSIVITY FOR THE EXTERNAL SURFACES 
C THE INTERNAL SURFACE IS ASSUMED BLACK 
C XMU THE SOLAR VIEW FACTOR TO HIE EXTERNAL SURFACES 
C A THE AREA OF EACH SIDE 
C FEB COEFFICIENTS TO EB; |FEB] |EB) = |Q) 

C X MATRIX CONTAINING THE ENVIRONMENTAL FLUX 
C = XSUN + XEARTH + XALBEDO 

C NUMINS NUMBER OF INSULATED SIDES 
C INSUL CONTAINS THE SIDES THAT ARE INSULATED (NO HEAT 
C FLOW THROUGH (E.q. RADIATOR SIDE) 

C NUMPTS NUMBER OF ORBIT POINTS 
C 

C XINC ORBIT INCLINATION (DEG) 

C OMEGA ORBIT LONGITUDE OR ASCENDING NODE (DEG) 

C WOMEGA LINE OF PERIAPSIS (DEG) 

C XNU TRUE ANAMOLY AT EPOCH (DEG) 

CECC ECCEN1RICI1Y 
CSMA SEMI- MAJOR AXIS 
C 

C INERT INERTIAL FLAG 
C I - INERTIAL 

C 0 - PLANET 

C 

C XS IS THE LOCAL SUNVECTOR 
C GS IS THE GLOBAL SUNVECTOR 


DIMENSION Q(6),EB(6).T(6) 

DIMENSION F(6,6). FE(6) 

DIMENSION ALPIR(6), ALPS(6). EPS(6) 

DIMENSION XMU(6). A(6) 

DIMENSION EEB(6.6), XRAD(6,6). XRADI(6,6), X(6) 

DIMENSION INSUL(5) 

DIMENSION XS(3), CS(3) 

DIMENSION 1AVE(6).QAVE(6) 

C OPEN THE OUTPUT FILE "RAO_ORB.DAT" CONTAINS THE 
C INFORMATION FOR EACH OF THE NUMPTS ORBIT POINTS 
OPEN(IO,FILE= 'RAO.ORB.DAT', STATUS= 'NEW') 

C READ INPUT FILE 

CALL ORBREADfAL BEDO.S.E.GS.XINC, OMEGA. WOMEGA.ECC.SMA.INERT.X.Y.Z, 
+ RE . ALPIR, ALPS.EPS.NUMINS.INSUL .NUMPTS) 

C INITIALIZE ANALYSIS VALUES 

CALL orbinit(sigma,xrad.xraoi.feb.tave.oave.tdelta.odelia, 

+ THIGH, OHIGH.TLOW.OLOW) 

C SET UP THE ORBITAL INTERVAL 
TWOPI = 2*3.141592654 
XINIERV = 1W0PI/NUMPIS 

C CALCULATE SURFACE AREAS 
CALL AREA(X.Y.Z.A) 

C CALCULATE ENCLOSURE VIEW FACTORS 
CALL CALCVWF(X,Y,Z,A,F) 

C MULTIPLIES A MODIFIED F MATRIX WITH THE AREA MATRIX TO 
C PRODUCE THE COEFFICIENTS TO THE EB MATRIX 
CALL CALFEB(F.A.FEB) 

C INITIALLY OEFINE THE TRUE ANOMALY = - LINE OF PERIAPSIS 
C THIS MAKES XNU GO FROM 0 TO 360 0 EG. THE SPACECAET MOVES 
C TO THE LEFT. 

C START AT LINE OF NODES 
XNU = -WOMEGA 
C 
C 

DO 101= I. NUMPTS 

C CALCULATE THE NORMAL AND TANGENTIAL EARTH VIEW FACTORS 
CALL CALVFNOR(ECC.SMA,XNU,RE,FEN) 

CALL CALVFTAN(ECC.SMA.XNU.RE.FET) 



C CALCULATE THE EFFECTIVE SPACECRAFT POSITION 
W = (WOMEGA+XNU) 

C CALCULATE THE EARTH VIEW FACTORS FOR INERTIAL OR 
C PLANET ORIENTEO SPACECRAFT 
IF (INERT.EO.I) THEN 
CALL CALFEI(W.FEN.FET.FE) 

ELSE 

CALL CALFEP(FEN.FEl.FE) 

ENDIF 

C DETERMINE IF THE SPACECRAFT IS IN THE SHADOW. THE SOLAR FLUX 
C RETURNED IS THE EFFECTIVE FLUX AND VIEW FACTORS 

CALL SHAOOW(ECC.SMA,XNU.XINC, OMEGA, WOMEGA.GS.S.SEFF.XMU.RE) 
C IF NOT IN THE SHADOW 
C CALCULATE THE SUN VECTOR 

C CALCULATE THE SOLAR VIEW FACTORS 

IF (SEFF.GT.O) THEN 

C IF THE SPACECRAFT IS INERTIALLY ORIENTED. THE SOLAR VIEW FACTORS 
C ARE AT THE LINE OF NODES. THEY MUSI BE RECALCULATED IN CASE 
C THE SPACECRAFT WAS PREVIOUSLY IN THE SHADOW. 

IF (INERT .EQ. I ) THEN 
CALL CALSUN(GS.XINC.OMEGA,O.O.XS) 

ELSE 

CALL CALSUN(CS.XINC.OMEGA.WOMEGA.XNU.XS) 

ENOIF 

CALL CALSOLR(XS.XMU) 

ENDIF 

C CALCULATE THE SOLAR FLUX 

CALL CALFLUX(SEFF,XMU,A,ALPS.FE,E.ALPIR .ALBEDO .EPS.X.XRAD, 

+ XRADI) 

C CALCULATE THE 0 ON EACH SIDE TAKING INTO ACCOUNT DECOUPLED SIDES 
CALL CALCO(X.INSUL.NUMINS.FEB.XRADI.O) 

C CALCULATE THE BLACKBODY RADIATION 
CALL CALCEB(X.Q.XRADI.EB) 

C CALCULATE THE TEMPERATURE 
CALL CALCT(EB,SIGMA,1) 

C PRINT THE DATA TO THE FILE "RAD.ORB.DAT" 

CALL ORBOUT(Q.T.XINC,OMEGA.WOMEGA.XNU) 

C COMPUTE THE STATISTICAL INFORMATION 

CALL ORBCOMP(T,Q.XNU.TAVE,QAVE.THIGH,TLOW.SHIGH.$LOW,VOHIGH. 


+ VOLOW.OHIGH.OLOW.SOHIGH.SOLOW.VOOHIGH.VOOLOW, 

+ TDEL T A.QDE L T A.VODEL T .VOQDE L ) 

C INCREMENT THE ANGLE 

XNU = XNU - XINTERV 
10 CONTINUE 

C CALCULATE THE FINAL SPACECRAFT STATISTICS 
CALL ORBFIN(NUMPTS.TAVET.OAVET.TAVE.OAVE) 

C WRITE DATA TO "RETURN.DAT" 

CALL ORBWRIT(TAVET .OAVET ,T AVE ,QAVE .THIGH ,T LOW.SHIGH .SLOW, 

+ VOHIGH.VOLOW.OHIGH.OLOW.SOHIGH.SOLOW.VOOHIGH.VOOLOW, 

+ TDELTA.ODELTA.VODELT.VOODEL) 

CLOSE(IO) 

END 
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QENCL.FOR 


PROGRAM OENCL 
IMPLICIT REAL *8 (A-H. 0-2) 

C THIS PROGRAM SOLVES FOR THE TEMPERATURE AND 
C HEAT LOAOS ON AN ENCLOSURE. IT IS USEO IN 
C COmC TION WITH THE EXPERT SYSTEM RADIATE.RUl 
C RACHEL ORMSBY June 29. 1993 
C 

C SOLVES FOR 0 INSTEAD OF EB 
C 

C 0(6) THE HEAT LOAD ON EACH SIDE OF A RECTANGULAR BOX 
C EB(6) THE BLACKBODY RADIATION ON EACH SIDE OF THE BOX 
C F(6.6) THE VIEW FACTORS INSIDE THE BOX 
C FE(6) THE VIEW FACTORS TO EARTH FOR EACH SIDE 
C ALPIR THE INFRARED ABSORPTIVITY FOR THE EXTERNAL SURFACES 
C = EMISSIVI1Y 

C ALPS THE SOLAR ABSORPTIVITY FOR THE EXTERNAL SURFACES 
C EPS THE EMISSIVITY FOR THE EXTERNAL SURFACES 
C THE INTERNAL SURFACE IS ASSUMEO BLACK 
C XMU THE SOLAR VIEW FACTOR TO THE EXTERNAL SURFACES 
C A THE AREA OF EACH SIDE 
C FEB COEFFICIENTS TO EB; |FEB) [EB] = JO] 

C X MATRIX CONTAINING THE ENVIRONMENTAL FLUX 
C = XSUN + XEARTH + XAL8ED0 

C NUMINS NUMBER OF INSULATED SIDES 
C INSUL CONTAINS THE SIDES THAT ARE INSULATED (NO HEAT 
C FLOW THROUGH (E.q. RADIATOR SIDE) 

DIMENSION Q(6),E8(6) 

DIMENSION F(6,6), FE(6) 

DIMENSION ALPIR(6), ALPS(6). EPS(6) 

DIMENSION XMU(6). A(6) 

DIMENSION FEB(6,6), XRAD(6.6), XRADI(6.6). XE(6) 

DIMENSION INSUl(6) 

C READ INPUT FROM EXPERT SYSTEM 

CALL ENCLREAD(XX,YY,ZZ.FE,ALPIR,AIPS,EPS.XMU.AL BE DO .NUMINS, 
+ INSUL.S.E) 

C INITIALIZES ANALYSIS VALUES 

CALL ENCINIT ($IGMA,XRAD,XRADI,FE8) 


CALL AREA(XX,TT.ZZ.A) 

C CALCULATES ENCLOSURE VIEW FAC10RS 
CALL CALCVWF(XX.YY.ZZ.A.r) 

C MULTIPLIES A MODIFIED F MATRIX WITH THE AREA MATRIX TO 
C PRODUCE THE COEFFICIENTS TO THE EB MATRIX 
CALL CALFEB(F.A.FEB) 

C CALCULATES THE ENVIRONMENTAL FLUXES 

CALL CALFLUX(S.XMU, A, ALPS, FE.E.ALPIR, ALBEDO, EPS, XE.XRAO.XRADl) 

C CALCULATES THE HEAT FLUX 

CALL CALCQ(XE,INSUL.NUMINS,FE8.XRADI.Q) 

C CALCULATES THE BLACK80DY RADIATION 
CALL CALCEB(X[,Q,XRA0I,EB) 

C CALCULATES THE TEMPERATURES 
CALL CALCT(EB.SIGMA.T) 

C WRITES THE OUTPUT TO "RETURN.DAT" 

CALL ENCWRIT(O.T) 

END 


C CALCULATES SURFACE AREAS 



SUBROUTINE ENCLREAD(X.Y.Z.FE.ALPIR,ALPS,EPS.XMU.AlBEDO. 
+ NUMINS.INSUL.S.E) 

IMPLICIT REAL *8 (A-H. O-Z) 

C 

C READS THE INPUT 
C 

C X THE X DIMENSION OF THE RECTANGULAR BOX (m) 

C Y 1HE Y DIMENSION OF THE RECTANGULAR BOX (m) 

C 2 THE Z DIMENSION OF THE RFCimUlAR BOX (m) 

C ALBEDO THE AL8EDO COEFFICIENT (USUALLY .3) 

CEE THE EARTH VIEW FACTORS 
C ALPIR THE INFRARED ABSORPTIVITY = EMISSMTY 
C ALPS THE SOLAR ABSORPTIVITY 

C XMU THE SOLAR VIEW FACTOR 

C NUMINS THE NUMBER OF INSULATED SIDES 
C INSUL THE SIDES THAT ARE INSULATED 
C S SOLAR FLUX 

C E PLANET FLUX 

C 

DIMENSION FE(6),ALPIR(6).ALPS(6).EPS(6).XMU(6),INSUL(6) 


DO 50 1=1,6 
READ(3.») XMU(I) 

50 CONTINUE 

READ(3.») NUMINS 
IF (NUMINS.GT.O) THEN 
DO 60 l=T, NUMINS 
READ(3,») INSUL(l) 
60 CONTINUE 
ENDIF 

CL0SE(3) 

RETURN 

END 


0PEN(3,FILE='RADIATE.DAT',STATUS='0LD') 

READ(3,») ALBEDO 
READ(3.*) S 
READ(3.<) E 


READ(3,») X 
READ(3.<) Y 
READ(3,») Z 


DO 101=1,6 
READ(3,0 FE(I) 

10 CONTINUE 

DO 20 1=1,6 
READ(3,») ALPIR(l) 
EPS(l) = ALPIR(I) 
20 CONTINUE 

00 301=1,6 
READ(3,«) ALPS(l) 
30 CONTINUE 
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SUBROUTINE ENCWRIT(Q.T) 

IMPLICIT REALMS (A-H, 0-Z) 

C 

C WRITES TEMPERATURE AND HEAT LOAD TO "RETURN.DAT" 
C 

DIMENSION 0(6), T(6) 

OPEN(10.EIIE='RETURN.OAT.STATUS= 'UNKNOWN') 

WRITE(I0,*) ‘ TTj’,T(l) 

WRITE] 10.') ' T2)',T(2) 

WRI1E(10,») ' 13]',T(3) 

WRIlE{tO.*) ' 14]',T(4) 

WRITE(IO,») ' T5]',T{5) 

WRITE(IO,<) ’|T6]',T(6) 

WRUEfTO,*) 

00 101 = 1,6 

IF (ABS(0(l)) .IT. I.OE-6) THEN 
0 ( 1 ) = 0.0 
ENDlf 

10 CONTINUE 

WRlTEftO,*) , (01]\0(l) 

WRITEflO,*) ' G2]'.Q(2) 

WRITE(IO,*) j03]',0(3) 

WRITE{10.«) ' 04]', 0(4) 

WRITE(IO,») ' Q5]'.Q(5) 

WRITEflO.*) ' 06]'.0(6) 

WRITEflO.*) 

RETURN 

END 
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QENCLSUB.FOR 

SUBROUTINE AREA(X.Y.Z,A) 

IMPLICIT REAL.8 (A-H. O-Z) 

C THIS SUBROUTINE CALCULATES THE AREAS Or EACH SIDE 
C THE SIOE AXES ARE OEEINEO. 


SUBROUTINE CALE£B(f,A,ffB) 

IMPLICIT REAL *8(A-H. O-Z) 

C CALCULATES THE FEB MATRIX 
C FEB = VIEW FACTORS • AREAS 
C IT IS THE COEFFICENT MATRIX TO THE BLACKBODY 
C RADIATION IERMS. EB (SEE DOCUMENIAIION) 
DIMENSION F(6,6), FEB(6.6), A(6) 


DIMENSION A(6) 

A(l) = X*Y 
A(2) = Z*X 
A(3) = Z»Y 
A(4) = Z*X 
A(5) = Z*Y 
A(6) = X*Y 

RETURN 

END 


C F THE ENCLOSURE VIEW FACTOR MATRIX. 

C A THE AREA MATRIX 

C THE DIAGONALS OF THE FEB MAT IX ARE THE SUM OF ALL THE 
C VIEW FACTORS FOR ONE SIDE (=1) TIMES THE AREA FOR THAT SIDE. 
C THE REMAINING POSITIONS ARE -F TIMES THE AREA 
DO 101=1,6 
DO 20 J= 1.6 
IF (I.EQ.J) THEN 
FEB(I.I) = A(l) 

ELSE 

FEB(I.J) = -F(l.J)* A(l) 

ENDIF 

20 CONTINUE 
10 CONTINUE 


RETURN 

END 
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SUBROUTINE CALFLUX(S.XMUAALPS.FE.E.AIR.AlB.EPS.X, 
4 XRAD, XRADI) 

IMPLICIT REAL'8(A-H.O-Z) 

C 

C CALCULATES THE ENVIRONMENTAL FLUXES 
C FILLS THE X MATRIX ANO THE XRAD MATRIX 
C 

C S SOLAR FLUX (W/m*2) 

C XMU(6) SOLAR VIEW FACTORS 
C A AREA MATRIX 
C ALPS(6) SOLAR ABSORPTIVITY 
C FE(6) EARTH VIEW FACTORS 
C AIR(6) INFRARED ABSORPTIVITY 
C ALB ALBEDO COEFFICIENT (USUALLY .3) 

C EPS(6) EMMISIVITY 

C X(6) ENVIRONMENTAL RADIATION ON EACH SIDE 
C XRA0(6.6) INTERNAL HEA1 TO BE RADIATED 
C XRAOl(6,6) I /XRAD 

DIMENSION A(6).FE(6),AIR(6).XMU(6).EPS(6). ALPS(6) 
DIMENSION X(6). XRAD(6,6). XRADI(6.6) 

DO 10 1=1,6 

OSUN = S'XMU(I).A(I)<AIPS(I) 

OEARIH = A(l)*FE(l)*E»AIR(l) 

0AL8ED0 = A(l)«FE(l)*f I -XMU(l))*S«ALB*ALPS(l) 
ORAO = (A(I).EPS(I)) 

X(l) = OSUN 4 OEARTH 4 QALBEDO 
XRAD(I.I) = -ORAD 

C XRAD AND XRADI WERE INITIALIZED TO 
C A ZERO MATRIX IN THE ROUTINE "ENCINIT" 

C SINCE XRAD IS A DIAGONAL MATRIX, THE INVERSE IS THE 
C RECIPRICOL OF EACH OF THE DIAGONAL TERMS. 

XRADI(IJ) = I /XRAD(l.l) 

10 CONTINUE 


SUBROUTINE CAICEB(X.Q.XRADI.EB) 

IMPLICIT REAL *8(A- H.O- Z) 

C 

C CALCULATES THE BALCKBODY RADIATION FOR EACH SURFACE 
C X IS THE ENVIRONMENTAL FLUXES 
C 0 IS THE CALCULATED HEAT LOAD 
C XRADI IS THE INVERSE OF THE RAOIATEO HEAT 
C OJt] = 10] - |X] 

C EB] = (XRADI] » [Q_X] 

C 

DIMENSION X(6),Q(6), EB(6,6). XRADI(6,6) 

DIMENSION 0_X(6) 

C IFAC = -1 SUBRTRACT MATRICES 
IFAC = - I 

CALL MATADD(6,1 .0.X.O-X.IFAC) 

C EB MATRIX = I /XRAD * (Q-X) OR 0 = X 4 XRAD'EB 
CALL MATMUL(6.6,I,XRADI.0-X,EB) 

RETURN 

END 


RETURN 

ENO 
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SUBROUTINE CALCO(X.INSUL.NUMINS.FEB.XRADI.O) 
IMPLICIT REAL«8(A-H,0-Z) 

CALCULATES THE HEAT LOAD ON EACH SURFACE 

DIMENSION FEBXRADI(6,6).B(6).IWK(6) 

DIMENSION X(6).0(6), FEB(6.6). XRADI(6,6) 

DIMENSION Xl(6,6), A(6.6) 

DIMENSION INSUL(6) 

C SET UP THE IDENTITY MATRIX |XI] 

DO 101=1.6 
DO 20 J= 1 ,6 
IF (l.CQ.J) THEN 
XI(I.J) = 1.0 
ELSE 

XI(I,J) = 0.0 
ENDIF 

20 CONTINUE 
10 CONTINUE 

C | FEB) = MATRIX OF VIEW FACTORS (|0) = |FEB] |EB]) 

C |XRADI) = INVERSE OF |XRAD] (|0) = |X] + )XRAD)[EB]) 
C |X] IS THE ENVIRONMENTAL LOADS 

C 

C DEFINE FOR A'Y = B 
C A) = (FE8) |XRAOI] - |XI] 

C B = (FEB) |XRADI] |X] 

C X] = 

C 

C IF NUMBER OF INSULATED SIDES > 0 THEN ZERO 
C THEN ZERO 0 BY SETTING X(SIDE) = 0 AND 
C XRAD(SIDE) = 0.0 OR XRADI(SIDE) = 9.9E+9 

IF (NUMINS.G1.0) THEN 
DO 30 K=I,NUMINS 
I = INSUL(K) 

XRADI(l.l) = 9.9E+9 
X(l) = 0.0 
30 CONTINUE 
ENDIF 

C [FEBXRADI] = [FEB) (XRADI) 

CALL MATMUL(6,6.6.FE8.XRA0I.FEBXRA0I) 
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C CALCULATE A AND B 
IF AC = -1 

CALL MATADD(6.6,FEBXRADI.XI.A.IFAC) 
CALL MATMUL(6,6. 1 .FEBXRADI.X.B) 

C SOLVE A.Y = B 
INFO = 0 
JOB =0 

CALL OGEFA(A.6.6.IWK.INFO) 

CALL DGESL(A.6.6,IWK,B,J0B) 

C FILL THE 0 VECTOR WITH THE SOLUTION 
00 801=1.6 
0(1) = 8(1) 

80 CONTINUE 

RETURN 

END 


Gil 



non 
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SUBROUTINE CALCT(EB.SIGMA.T) 

IMPLICIT REAL*8 (A-H. 0-Z) 

CALCULATES THE TEMPERATURES GIVEN EB 

DIMENSION EB(6).T(6) 


SUBROUTINE CALCVWF(X.Y.Z,A,F) 

IMPLICIT REAL'8 (A-H. O-Z) 

C 

C CALCULATES VIEW FACTORS TOR THE INTERNAL ENCLOSURE 
C 

DIMENSION F(6,6), A(6) 


C EB = SIGMA.T-4 
DO 101=1,6 

1(1) = SORT(SORT(EB(l)/SIGMA)) 
10 CONTINUE 

RETURN 

END 


E(l,l) = 0.0 
f(2,2) = 0.0 
F(3,3) = 0.0 
F(4,4) = 0.0 
F(5,5) = 0.0 
F(6.6) = 0.0 


C CALCULATE VIEW FACTOR ON THREE SIDES, USE RECIPROCITY 
C AND SUM OF VIEW FACTORS =t FOR THE REST 
F(1 ,3) = VIEWFAC(X.Z.Y) 

F(2,l) = VIEWFAC(Z.Y.X) 

F(3,2) = VIEWFAC(Y.X.Z) 

F(I,5) = F(I,3) 

F(2,5) = F(2.I) 

F(2,6) = F(2.l) 
r(3,4) = F(3,2) 

F(T.I) = F(2,l) 

F(4,3) = F(3,2) 

F(4,6) = F(2,I) 

F(5.4) = F(3,2) 

F(6,3) = F(l J) 

F(6,5) = F( 1 ,3) 


F( 1 .2) = F(2,1 )* (A(2)/A(l)) 
F(2,3) = F(3,2)< (A(3)/A(2)) 

f(3,i) = f(i,3). mm ) 

F(4,5) = F(5,4). mm) 


F(’.4) = F(1,2) 

F(3,6) = F(3,l) 
F(5,l) = F(3,l) 
F(5.2) = F( 1 .2) 
F(5,6) = F(3.l) 
F(6.2) = F( 1 ,2) 
F(6,4) = F(l,2) 



F(l.6)= I - (F{ t .2)+ F( t .3)+ F{ t .4)4 T( 1 .5)) 

F(2,4) = l-(F(2.1)+F(2.3)+F(2.6)+F(2,5)) 

F(3.5) = !-(F(3.l)+F(3.2)+F(3.4)4F(3.6)) 

F(4.2) = 1-(F(4,1)+F(4.3)4F(4.6)+F(4,5)) 

F(5,3) = 1-(F(5,1)+F(5,2)+F(!j.4)+F(5,6)) 

F(6,l) = 1-(F(6.2)+F(6.3)+F(6.4)+F(6.5)) 

C 

C THE FOLLOWING COMMENTED CODE ALLOWS THE DESIGNER 
C TO WRITE OUT THE ENCLOSURE VIEW FACTORS. FOR INSPECTION 
C OPEN ( 1 0.FILE= OUT DAV.ST ATUS= ’UNKNOWN') 

C 00 101=1.6 

C DO 20 J= 1 .6 

C WRITEf 1 0,*) l,J.F(l,J) 

C 20 CONTINUE 
C TO CONTINUE 

RETURN 

END 


SUBROUTINE ENCINIT(SIGMA.XRAD.XRADI.FEB) 

IMPLICIT REAL*8(A- H.O-Z) 

C 

C INITIALIZES THE CONSTANTS AND SETS THE MATRICES TO 0.0 
C 

DIMENSION FEB(6.6), XRAD(6.6). XRADI(6.6) 

SIGMA = 5.67E-8 

DO 101=1,6 
DO 20 J= 1.6 
XRAD(I.J) = 0.0 
XRADI(I.J) = 0.0 • . 

FEB(I.J) = 0.0 
20 CONTINUE 
10 CONTINUE 

RETURN 

END 
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FUNCTION VIEWFAC(HI.WI.XL) 

IMPLICIT REAL <8 (A-H, 0-2) 

C 

C CALCULATES THE VIEWFACTOR FOR TWO FINITE PLATES AT 90 OEG 
C TAKEN FROM "THERMAE RADIATION HEA1 IRANSfER" BY ROBERT SIEGEL AND 
C JOHN R. HOWELL. APPENDIX C #14 "TWO FINITE RECTANGLES OF SAME LENGTH, 
C HAVING ON COMMON EDGE. AND A1 AN ANGIE OF 90 DEG 10 EACH OTHER." 

C 

XPI = 3.141592654 

H = HI/XL 
W = Wl/Xl 

A = l/(XPUW) 

B = W*ATAN(1/W) 

C = H>ATAN(I/H) 

HH = H'H 
WW ^ W.W 
0 = HH + m 
E = SQRI(D) 

F = E*AIAN(I/E) 

G = l+WW 
0 = l+HH 
P r 1+HH+WW 
PP = HH+WW 
0 = ,25*0L0G(G*0/P) 

R = .25>DLOG(WW*P/(G»PP» 

S = ,25*0LOG(HH<P/(O'PP)) 

T = B+C-F+Q+R4S 
U = Atl 


VIEWFAC = U 

RETURN 

END 



EARTHV.FOR 


PROGRAM EARIHV 
IMPLICIT REAL*8 (A-H, O-Z) 

C THIS PROGRAM CALCULATES THE EARTH VIEW FACTORS 
C AT ANY POINT IN AN ELLIPTICAL SOLAR INERTIAL ORBIT. 

C GIVEN THE POSITION IN ORBIT 
C ASSUMING AT TIME 0 THE X AXIS POINTS FROM NADIR. 

C THE Y AXIS POINTS IN THE ORBIT DIRECTION AND 
C THE Z AXIS COMPLETES THE RIGHT HAND RULE POINTING NORTH. 

C 

C THIS PROGRAM IS USEO WITH THE EXPERT SYSTEM RADIATE. RUL 
C RACHEL ORMSBY July 15. 1993 
C 

DIMENSION FE(6) 

C READS INPUT 

CALL EAREAD(W,ECC .SMA.XNU .RE) 

C CALCULATES THE NORMAL AND TANGENTIAL COMPONENTS Of THE 
C EARTH VIEW FACTORS 

CALL CALVFNOR(ECC.SMA.XNU.RE.FEN) 

CALL CALVFTAN(ECC.SMA.XNU.RE.FET) 

C ULIINLS I TIL SPACECRAF I POSIIION AS MIL LINL 01 PLKIAPSIS. W. 

C MINUS THE TRUE ANOMOLY. NOTE THIS ORBITS PROCEDES TO THE LEFT. 
Wl = (W-XNU) 

C ROTATES THE VIEW FACTORS TO THE CORRECT SIDES 
CALL CALFEI(WI.FEN.FET.FE) 

CALL EARWRIT(FE.FET.FEN) 


SUBROUTINE EAREAD(W.ECC.SMA.XNU.RE) 
IMPLICIT REAL'8 (A-H. O-Z) 

C 

C READS THE INPUT 
C 

C W LINE OF PERIAPSIS 
CECC ECCENTRICITY 
CSMA SEMI- MAJOR AXIS 
C XNU TRUE ANOMOLY AT EPOCH 
C RE RADIUS OF THE PLANET 

DTR = I.7453293E-? 

0PEN(3,FILE= 'RADIATE.DAT', STA1US= 'OLD') 

READ(3,0 INC 
READ(3,<) ECC 
READ(3,») W 
READ(3,') SMA 
READ(3,«) XNU 
XNU = -XNU 
READ(3.0 RE 

INC = INC* DTR 
W W*I)IK 
XNU = XNU'DTR 

RETURN 

END 


END 
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SUBROUTINE EARWRfTfEE.fCT.rCN) 

IMPLICIT REAL *8 (A-H. 0-Z) 

c 

C WRITES PLANET VIEW FACTORS FOR EACH SIDE. AND THE 
C NORMAL AND TANGENTIAL COMPONENTS TO “RET URN .DAT” 
C 

DIMENSION FE(6) 

OPEN(lO,FILE=’RETURN.DAT',STATUS='UNKNOWN') 


WRITE! 10.0' 
WRITEflO.O ' 
WRITE! 10.0 ' 
WRITE! 1 0.*) ' 
WRITEflO.0 ' 
WRITE(IO,0 ' 
write! io,*) ' 
WRITE(10.0 ' 


FE 1 1 

FE2 

FE3 

FE4 

FE5 

FES 

FET] 

FEN’ 


’.FE(l) 

•FE(2) 

.FE(3) 

••FE(0 

\FE(5) 

.FE{6) 

'.FET 

'.FEN 


close! to) 


RETURN 

END 



CALVFNT.FOR 


SUBROUTINE CALVFNOR(ECC.SMA.XNU.RE,FEN) 

IMPLICIT REAL <8 (A-H, O-Z) 

C CALCULATES THE VIEW FACTOR FROM A PLANE ELEMENT TO 
C A SPHERE. USED AS THE EARTH VIEW FACTOR NORMAL 
C FROM "THERMAL RADIATION HEAT TRANSFER" BY SIEGEL AND 
C HOWELL: APPENDIX C 

C CALCULATE THE SPACECRAFT POSITION COEFFICIENT. R. 

CALL ORBPROPS(ECC.SMA,XNU.R) 

FEN = (RE/R)*(RE/R) 

RETURN 

END 


SUBROUTINE CALVFTAN(ECC.SMA,XNU.RE,FET) 

IMPLICIT REAL'8 (A-H, O-Z) 

C CALCULATES THE VIEW FACTOR FROM A PLANE ELEMENT TO 
C A SPHERE. USED AS THE EARTH VIEW FACTOR TANGENTIAL 
C FROM "THERMAL RADIATION HEAT TRANSFER" BY SIEGEL AND 
C HOWELL: APPENDIX C 

XPI = 3.141592654 

C CALCULATE THE SPACECRAFT POSITION COEFFICIENT. R. 

CALL ORBPROPS(ECC.SMA.XNU.R) 

H = R/RE 

H2 = SQRT(H.H-I) 

HT = ATAN(I/H2) 

FET = (l/XPl)*(HT *(H2/(H'H))) 

RETURN 

END 
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SUBROUTINE CALFEP(FEN,FET,FE) 

1MPUC11 REAL. 8 (A-H.O-Z) 

C CALCULATES THE EARTH VIEW FACTOR FOR PLANET 
C ORIENTAIION 

DIMENSION F[(6) 

FE(1) = FET 
FE(2) = FET 
FE(3) = 0.0 
FE(4) = FET 
FE(5) = FEN 
FE(6) = FEI 

RETURN 

END 
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SUBROUTINE CALFEI(W.FEN.FET.FE) 

IMPLICIT REAL«8 (A-H.O-Z) 

C CALCULATES THE EARTH VIEW FACTOR FOR INERTIAL 
C ORIENTATION 

DIMENSION FE(6) 

C W IS THE SPACECRAFT POSITION = LINE OF PERIAPSIS - 
C TRUE ANOMOLY 
SW = SIN(W) 

CW = COS(W) 

FE(1) = FET 

FE(2) = -FEN.SW + FET.ABS(CW) 

FE(3) = FEl'ABS(SW) - FEN'CW 
Ft(4) = rCN.SW + FET'ABS(CW) 

Fc(h) - m.AHs(sw) i riN>cw 
FE(6) = FEI 

C IF THE VIEW FACTOR IS LESS THAN ZERO. THE SPACECRAFT 
C SURFACE CANNOT SEE THE PLANET. 

DO 101=1,6 
IF (FE(I).LT.O.O) THEN 
FE(l) = 0.0 
ENDIF 

10 CONTINUE 

RETURN 

END 



SOLARV.FOR 


PROGRAM SOLARV 
IMPLICIT REAL*8 (A-H. O-Z) 

C THIS PROCRAM SOLVES TOR THE OOT PRODUCT OT 
C THE SOLAR VECTOR AND THE LOCAL UNIT NORMALS 
C ON EACH SIDE OF THE ENCLOSURE. 

C THIS PROGRAM IS USED WITH THE EXPERT SYSTEM RADIATE.RUL 
C RACHEL ORMSBY July 7, 1993 
C 

DIMENSION XMU(6),X(3) 


SUBROUTINE SOLREAO(X) 

IMPLICIT REAL'8 (A-H, O-Z) 

C 

C READS THE INPUT 
C 

C XS THE X COORDINATE OF HIE SOLAR VECIOR 
C YS THE Y COORDINATE OF THE SOLAR VECTOR 
C 2S THE Z COORDINATE OF THE SOLAR VECTOR 
C MAGS THE MAGNITUDE OF THE SOLAR VECTOR 

DIMENSION X(3) 

OPEN(3,FILE=’RADIATE.DAT',STATUS='OLD') 


CALL SOLREAD(X) 

CALL CALSOLR(X.XMU) 
CALL SOLWRIT(XMU) 


READ(3,*) X(l) 
READ(3,») X(2) 
READ(3.«) X(3) 


END 


RETURN 

END 
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SUBROUTINE SOLWRIT(XMU) 

IMPLICIT REAL *8 (A-H, 0-Z) 

WRITES BLACKBODY RADIATION, HEAT LOAD AND 
VIEW FACTORS TO "RETURN.OAT” 


OIMENSION XMU(6) 

OPEN(IO,FllE=’RETURN.DAT',SIATUS='UNKNOWN') 


WRI1E(I0,<) ' 
WRITEf 10.») ' 
WRI1E(I0,») ' 
WRI1E(10.») ' 
WRIIE(10,») ' 
WRITE(IO,») ' 


MU I 
MU2' 
MU3 1 
MU4 
MU5 
MU6 1 


’,XMU(l) 

\XMU(2) 

\XMU(3) 

\XMU(4) 

’,XMU(5) 

\XMU(6) 


CLOSE(IO) 

RETURN 

END 



CALSOLR.FOR 


SUBROUTINE CALSOLR(X.XMU) 

IMPLICIT REAL.8 (A-H. O-Z) 

C CALCULATES THE DOT PRODUCT OT THE SOLAR VECTOR AND THE 
C UNIT NORMALS OE EACH SURFACE 
DIMENSION XMU(6).X(3) 

C DEFINES THE SURFACE NORMAL ON EACH SIDE OT THE ENCLOSURE 
CALL SURFNOR(X.XMU) 

C IF THE VIEW FACTOR IS LESS THAN 0. THE SURFACE CANNOT SEE 
C THE SUN. 

DO 10 1=1.6 
IF (XMU(I).LT.O.O) THEN 
XMU(I) = 0.0 
ELSE 

XMU(l) = ABS(XMU(I)) 

ENDIF 

C IF THE VIEW FACTOR IS A NUMBER SMALLER THAT I.E-3. IT WILL 
C BE WRITTEN TO THE FILE IN SCIENTIFIC NOTATION. THEREFORE 
C ZERO IT OUT TO ELIMINATE THE ERROR. EXSYS* CANNOT READ 
C SCIENTIFIC NOTATION 

IT (XMU(I) .IT. I.OE-3) THEN 
XMU(l) = 0.0 

endit 

10 continue 

RETURN 

END 
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SUBROUTINE SURFNOR(X.Y) 

IMPLICIT REAL*8 (A-H, O-Z) 

C CALCULATES THE DOT PRODUCT OF X INTO THE LOCAL 
C SURFACE NORMALS AS DEFINED IN DOCUMENTATION 
C X IS A VECTOR IN LOCAL COORDINATES 
C Y IS THE DOT PRODUCT WITH THE SUREACE NORMALS 

DIMENSION Y(6),X(3) 

XS = X(l) 

YS = X(2) 

ZS = X(3) 

C CALCULATES THE MAGNITUDE OE THE SUN VECTOR. IN CASE 
C IT IS NOT I. XMAG IS A FUNCTION CALL... 

XMAGS = XMAG(X) 

C DEFINES THE SOLAR VIEW FACTORS. SINCE EACH SURFACE NORMAL 
C IS IN THE DIRECTION OF A COORDINATE AXIS. THE DOT PRODUCT IS 
C "HARD" CODED IN HERE. 

Y(l) = -ZS/XMAGS 
Y(2) = -YS/XMAGS 
Y(3) = -XS/XMAGS 
Y(4) = YS/XMAGS 
Y(5) = XS/XMAGS 
Y(6) = ZS/XMAGS 

RETURN 

END 
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SUNVECT.FOR 


PROGRAM SUNVECT 
IMPLICIT REAL'S (A-H, 0-Z) 

C 

C THIS PROGRAM CALCULATES THE SUN VECTOR IN LOCAL 
C COORDINATES, GIVEN INCLINATE. LONGITUDE OT 
C ASCENDING NODE AND SATELLITE POSITION. 

C THIS CODE IS USED WITH THE EXPERT SYSTEM 
C RADIATE.RULE. RAO 7/16/93 
C 

DIMENSION X(3).GS(3), XMU(6) 

CALL SUNREAD(GS.XINC.OMEGA.WOMEGA.XNU) 

C CALCULATES THE SUN VECTOR IN LOCAL COORDINATES 
CALL CALSUN(GS.XINC,OMCGA.WOMEGA,XNU.X) 

C CALCULATES THE SOLAR VIEW FACTORS 
CALL CALSOLR(X.XMU) 

CALL SUNWRIT(X.XMU) 

END 


SUBROUTINE SUNREAD(GS,XINC, OMEGA, WOMEGA.XNU) 
IMPLICIT REAL *8 (A-H. O-Z) 

C 

C READS THE INPUT 
C 

CGS(I) GLOBAL SUN VECTOR 
C XINC INCLINATION OT THE ORBIT (DEG) 

C OMEGA LONGITUDE OF ASCENDING NODE (DEG) 

C WOMEGA LINE OF PERIAPSIS (DEG) 

C XNU TRUE ANOMALY AT EPOCH (DEG) 

C 

DIMENSION GS(3) 

DTR = 1.7453293E-2 

OPEN(3,FILE= 'RADIATE, DA!'.STATUS= 'OLD') 

REA0(3,») GS(I) 

READ(3,») GS(2) 

REAO(3.0 GS(3) 

READ(3,<) XINC 
READ(3,<) OMEGA 
READ(3,*) WOMEGA 
READ(3,*) XNU 

C CONVERT TO RADIANS 
XINC = XINC* DTR 
OMEGA = OMEGA* DTR 
WOMEGA = WOMEGA'DIR 
XNU = XNU'DTR 

CLOSE(3) 

RETURN 

END 



o o <-> 


SUBROUTINE SUNWRIT(X.XMU) 

IMPLICIT REA1.8 (A-H. O-Z) 

WRITES THE THREE COMPONENTS OF THE SUN VECTOR 
IN LOCAL COORDINATES AND THE 6 SOLAR VIEW FACTORS 

DIMENSION X(3) 

DIMENSION XMl)(6) 

OPEN( 1 0,FIIE= 'RETURN.OAT'.STATUS= ’UNKNOWN') 

WRITE(10.*) • XS)'.X(I) 

WRITEflO.*) ' YS]'.X(2) 

WRITEflO,') ' ZS \X(3) 

WRITE(IO,') ' MUI)'.XMU(I) 

WRITE( 10.*) ' MU2]'.XMU(2) 

WRITE! I0.») ' MU3]',XMU(3) 

WRITEflO,*) ' MU4]'.XMU(4) 

WRITEflO.*) ' MU5]',XMU(5) 

WRITEflO,') ’ MU6]',XMU(6) 

CLOSE! 10) 

RETURN 

END 


FORTRAN Source Code 


G23 



FORTRAN Source Code 


G24 


CALSUN.FOR 


SUBROUTINE C ALSU N (G .X INC .OMEGA.WOME GA.XNU ,XS) 
IMPLICIT REAL *8 (A-H. O-Z) 

DIMENSION XS(3). A(3.3). B(3.3). C(3). GS(3). C(3.3) 
DIMENSION IPVT ( 3) 

C ROTATES THE GLOBAL SUNVECTOR INTO LOCAL COORDINATES 
C SOLVING THE EQUATION Ax = b 
C A IS THE COEFFICIENT MATRIX FROM PERIFOCAL TO GLOBAL 
C B IS THE COEFFICIENT MATRIX FROM PERIFOCAL TO SPACECRAFT 
C X IS THE LOCAL SUN VECTOR 
C G IS THE GLOBAL SUN VECTOR 

LDA = 3 
N = 3 
INFO = 0 
JOB = 0 

DO 101=1.3 
GS(l) = G(l) 

10 CONTINUE 

C DEFINE THE ROTATION MATRIX 

CALL PQR2UK(XINC.OMECA.WOMECA.XNU.A,B) 

CALL MATMUL(3.3.3.B.A,C) 

C SOLVE Ax = b 

CALL OGEFA(C.LDA,NJPVT.INFO) 

CALL DGESL(C,LDA.N.IPVT.CS.JOB) 

C IF THE VALUE IS LESS THAN T.OE-3 IT WILL BE WRITTEN 
C IN SCIENTIFIC NOTATION. WHICH CANNOT Ilf READ BY FXSYS* 
DO 20 1=1.3 

IF (DABS(GSfl)) IT. T.OE-3) THEN 
XS(I) = 0.0 
ELSE 

XS(I) = GS(I) 

ENDIF 

20 CONTINUE 


SUBROUTINE PQR2UK(XI.0,W,XN,A.B) 

IMPLICIT REAL *8 (A-H. O-Z) 

C THIS SUBROUTINE FILLS THE COEFFICIENT MATRIX 
C FOR THE ROTATION FROM PERIFOCAL SYSTEM (P.O.R) TO THE 
C GLOBAL SYSTEM (l.J.K) 

C SEE DOCUMENTATION FOR DERIVATION 
DIMENSION A(3,3), B(3,3) 

SXI = SIN(XI) 

SW = SIN(W) 

SO = SIN(O) 

SN = SIN(XN) 

CXI = COS(XI) 

CW = COS(W) 

CO = COS(O) 

CN = COS(XN) 

A(l,l) = CO'CW - SO<SW.SXI 
A( 1 .2) = -CO*SW - SO'CW.CXI 
A(l,3) = -(SO.SXI) 

A(2.l) = SO<CW l CO*SW.CXI 
A(2.2) = -SO*SW + COtCW.CXI 
A(2,3) = -(-CO.SXI) 

A(3.1) = -(SW.SXI) 

A(3,2) = -(CW.SXI) 

A(3,3) = CXI 

B(U) = CN 
B(l,2) = SN 
0(1.3) = 0 
B(2,l) = -SN 
8(2,2) = CN 
8(2,3) = 0 
U(3,l) = 0 
6(3,2) = 0 
B(3,3) = I 

RETURN 

END 


RETURN 

END 



MATMUL.FOR 


MATADD.FOR 


SUBROUTINE MATMUL(N,MXAB ( C) 
IMPLICIT REAL*8 (A-H.O-Z) 

C MULTIPLIES TWO MATRICES 
C A IS DIMENSIONED N,M 
C B IS DIMENSIONED MX 
C C IS DIMENSIONED NX 


SUBROUTINE MATADD(N.MAB.CJFAC) 
IMPLICIT REAL *8 (A-H.O-Z) 

C ADDS TWO MATRICES 

C ALL THREE MATRICES ARE DIMENSIONED NXM 

C IF IFAC =-1 THEN SUBRACTION 
C IF IFAC = I THEN ADDITION 


DIMENSION A(N.M), B(MX). C(NX) 


DIMENSION A(N.M), B(N.M). C(N.M) 


DO 10 I =1X 
DO 20 J = 1.N 
D = 0.0 
DO 30 K = l,M 
D = A(J,K)*B(K,I) 4 D 
30 CONTINUE 
C(J.I) = D 
20 CONTINUE 
10 CONTINUE 


DO 10 I = I.N 
DO 20 J = 1.M 

C(u) = A(l,J) 4 IFAC*B(I.J) 
20 CONTINUE 
10 CONTINUE 

RETURN 

END 


RETURN 

END 
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POSGLOB.FOR 


SUBROUTINE POSGLOB(ECC,SMA,XI,O.W,XNU, RVECTG) 

IMPLICIT REA1.8 (A-H. 0-Z) 

C CALCULATES THE ORBITAL POSITION IN GLOBAL COORDINATES 
DIMENSION RVEC1(3).RVECTG(3) 

DIMENSION A(3.3), B(3.3) 

C ECC ECCENTRICITY 
C SMA SEMI-MAJOR AXIS (KM) 

C XI ORBIT INCLINATION (RAD) 

C 0 LONCITUDE OF ASCENDING NODES (RAD) 

C W LINE OF PERIAPSIS (RAD) 

C XNU TRUE ANOMALY AT EPOCH (RAD) 

C 

C RVECTG GLOBAL POSITION VECTOR 

C ORBPROPS CALCULATES THE POSITION VECTOR MAGNITUDE R 
CALL ORBPROPS(ECC.SMA.XN.R) 

C CALCULATES THE POSITION VECTOR RVECT BY PROJECTING INTO 
C THE PERIFOCAL COORDINATE SYSTEM 
RVECT(I) = R'COS(XN) 

RVECT(2) = R'SIN(XN) 

RVECT(3) = 0.0 

C CALCULATES 1HE TRANSFORMATION MATRIX FOR LOCAL TO 
C GL08AL AND CALCULATES THE GLOBAL POSITION VECTOR 
CALL PQR2IJK(XI,0.W.XNUAB) 

CALL MATMUL(3.3.I. A, RVECT, RVECTG) 

RETURN 

END 



ORBSUBS.FOR 


10 CONTINUE 


SUBROUTINE ORBCOMP(T.Q.VO.TAVE.QAVE.1HIGH.TLOW.SHIGH.SLOW. 

4 VOHIGH.VOLOW, OHIGH, QLOW.SQHIGH.SQLOW.VOQHIGH.VOQLOW.TOELTA. 
4 OOELTA.VODELT.VOODEL) 

IMPLICIT REAL *8 (A-H.O-Z) 

C CALCULATES THE ORBITAL AVERAGE VALUES EOR TEMPERATURE 
C AND HEAT LOADS. 

C 

C 1(6) CURRENT TEMPERATURE 
C 0(6) CURRENT HEAT LOAO 
CVO TRUEANOMOLY 
C TAVE(6) AVERAGE TEMPERATURE 
C QAVE(6) AVERAGE HEAT LOAO 
C 

C THIGH HIGHEST TEMPERATURE 
C TLOW LOWEST TEMPERTURE 
CSHIGH SIDE WITH THIGH 
C SLOW SIDE WITH TLOW 
C VOHIGH VO EOR THIGH 
CVOLOW VO FOR TLOW 
C 

C OHIGH HIGHEST HEAT LOAD 
COLOW LOWEST HEAT LOAD 
CSOHIGH SIDE WITH OHIGH 
CSOLOW SIDE WITH OLOW 
C VOOHIGH VO FOR OHIGH 
C VOOLOW VO FOR OLOW 
C 

C TDELTA HIGHEST TEMPERATURE DELTA BETWEEN ANY TWO SIDES 
C ODELTA HIGHEST HEAT LOAO DELTA BETWEEN ANY TWO SIDES 
CVODELT VO FOR TDELTA 
CVOODEL VO FOR ODELTA 
C 

C TAVE AND QAVE ARE THE SUM OF ALL TEMPERATURES AND HEAT 
C LOADS DIVIDED BY THE NUMBER OF POINTS FOR EACH SIDE. 

DIMENSION T(6), 0(6) 

DIMENSION TAVE(6),OAVE(6) 

DO 101=1,6 
TAVE (I) = TAVE(I) 4 T(l) 

OAVE(I) = OAVE(l) 4 0(1) 


C INITIALIZE THE HIGHEST T AND 0 AS 0 
C LOWEST T AND 0 AS 1000.0 
TH = 0.0 
TL = 1000.0 
OH = 0.0 
OL = 1000.0 

C IF THE TEMPERATURE IS HIGHER THAN THE CURRENT HIGHEST 
C TEMP, SAVE THE VALUE AS THE CURRENT HIGH. 

DO 20 1=1,6 
IF (T(I).GT.TH) THEN 
TH = T(l) 

ENDIF 

IF (T(l).LT.Tl) THEN 
TL = 1(1) 

ENDIF 

C IF THE ABSOLUTE HEAT LOAD IS HIGHER THAN THE CURRENT HIGHEST 
C ABSOLUTE HEAT LOAD, SAVE THE VALUE AS THE CURRENT HIGH. 

IF (ABS(Q(I)).GT.ABS(QH)) THEN 
OH = 0(1) 

ENDIF 

IT (ABS(0(I)).L T .ABS(OL)) THEN 
OL = 0(1) 

ENDIF 

20 CONTINUE 

C IF DELTA IS LESS THAN THE DIFFERNCE BETWEEN THE HIGH AND 
C LOW, SAVE THE NEW DELTA AND SAVE THE POSITION 
IF (TDELTA.LT.(TH-TL)) THEN 
TDELTA = TH-TL 
V0DEL1 = VO 
ENDIF 

IF (ODELTA.LT.ABS(OH-OL)) THEN 
ODELTA = OH-OL 
VOODEL = VO 
ENDIF 

C CHECK EACH SIDE FOR A HIGHER OR LOWER TEMP THAN THE CURRENT 
C HIGHEST AND LOWEST. 

DO 30 1=1,6 
IF (l(l).GT.THIGH) THEN 
THIGH = T(I) 
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SHIGH = I 
VOHIGH = VO 
ENOIf 

If {T(I).LT.ILOW) THEN 
ROW = 1(1) 

SLOW = I 
V010W = VO 
ENOIf 

IF (A8S(Q(l)).Gl.ABS(OHIGH» 1HEN 
OHIGH = 0(1) 

SOHIGH = I 
VOOHIGH = VO 
ENOIf 

IF (ABS(O(O).LI.ABS(OLOW)) THEN 
QlOW = 0(1) 

SOI.OW = I 
VOQLOW = VO 
ENOIf 

30 C0N1INUE 

RETURN 

END 


SUBROUTINE ORBf IN(NUMPTS.TAVET.OAVET.TAVE.OAVE) 

C FINISHES THE ORBITAL CALCULATIONS 
IMPLICIT REAL *8 (A-H.O-Z) 

OIMENSION TAVE(6),OAVE(6) 

C NUMP1S NUMBER Of ORBIT POINTS CALCULATED 
C 1AVE1 TOTAL AVERAGE TEMPERATURE 
C OAVET TOTAL AVERAGE HEAT LOAD 
C TAVE(6) AVERAGE TEMPERATURE FOR EACH Of THE SIX SIDES 
C OAVE ( 6) AVERAGE HEAT LOAD FOR EACH Of THE SIX SIOES 

00 10 1 = 1,6 

TAVE(I) = 1AVE(I)/NUMP1$ 

OAVE(I) = TAVE(l)/NUMPTS 
10 CONTINUE 


TAVE1 = 0.0 
OAVE I = 0.0 
DO 20 1-1,6 

TAVEI = TAVET + TAVE(I) 
OAVET = OAVET + OAVE (I) 
20 CONTINUE 
TAVET : TAVET/6 
OAVET = OAVET/6 

RETURN 

END 



o o <-> 


SUBROUTINE ORBINIT(SIGMA,XRAD,XRADI.FEB.TAVE.QAVE.TDELTA. 
+ ODEUA.IHICH.OHIGH.ILOW.OLOW) 

IMPLICIT REAL<8(A-H.0-Z) 

INITIALIZES THE CONSTANTS AND SETS THE MATRICES TO 0.0 

DIMENSION FEB(6,6), XRAD(6.6). XRADl(6,6) 

DIMENSION TAVE(6), 0AVE(6) 

SIGMA = 5.67E-8 

THIGH = 0.0 
OHIGH = 0.0 
ROW = 1000.0 
GLOW = 1000.0 
TDELTA = 0.0 
ODELTA = 0.0 

DO 101=1,6 
DO 20 J=1,6 
XRAD(I.J) = 0.0 
XRADI(I.J) = 0.0 
FEB(I,J) = 0.0 
20 CONTINUE 

TAVE(I) = 0.0 
OAVE(I) = 0.0 
10 CONTINUE 

RETURN 

END 


SUBROUTINE 0RBREAD(ALBED0.S,E,GS,XI,0,W0.ECC.SMA.INERT .X.Y.Z. 
+ RE.ALPIR.ALPS.EPS.NUMINS.INSUL.NUMPTS) 

IMPLICIT REAL»8 (A-H, 0-Z) 

C 

C READS THE INPUT 
C 

C ALBEDO THE ALBEDO COEFFICIENT (USUALLY .3) 

C GS HIE GLOBAL SUN VLCIOR 
C INC ORBIT INCLINATION (DEG) 

C OMEGA LONGITUDE OF ASCENDING NODES (DEG) 

C WOMEGA LINE OF PERIAPSIS (DEG) 

C ECC ECCENTRICITY 
CSMA SEMI- MAJOR AXIS 

C X THE X DIMENSION OF THE RECTANGULAR BOX 
C Y THE Y DIMENSION OF THE RECTANGULAR BOX 
C Z THE Z DIMENSION OF THE RECTANGULAR BOX 
C ALPIR THE INFRARED ABSORPTIVITY = EMISSMTY 
C INERT INERTIAL FLAG I -INERTIAL O-PLANET 
C ALPS THE SOLAR ABSORPTIVITY 
C NUMINS THE NUMBER OF INSULATED SIDES 
C INS THE SIDES THAT ARE INSULATED 
C NUMPTS NUMBER OF ORBIT POINTS 
C 

DIMENSION ALPIR(6),ALPS(6),EPS(6),INSUL(5),GS(3) 

OTR = 1 .7453292E-2 

0PEN(3.FILE=’RADIATE.DAT , .STATUS='0LD’) 

REA0(3,*) ALBEDO 
READ(3,») S 
READ(3.») E 

DO 40 1=1.3 
READ(3,») GS(I) 

40 CONTINUE 

REA0(3,») XI 
READ(3,*) 0 
READ(3.*) WO 
READ(3,») ECC 
READ(3,») SMA 

XI = XHDTR 
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0 = O'DTR 
WO = WO'DTR 

READ(3.») INERT 

READ(3.») X 
REA0(3.») Y 
READ(3.«) Z 
REA0(3,') RE 

00 20 1 = 1.6 
REA0(3,<) AEPIR(I) 
EPS(I) = ALPIR(l) 
20 CONTINUE 

00 30 1=1.6 
READ(3,«) ALPS(I) 
30 CONTINUE 

READ(3,>) NUMINS 
00 60 1=1.5 
REA0(3,*) INSUL(I) 
60 CONTINUE 

READ(3.») NUMPTS 

CLOSE(3) 

RETURN 

END 
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SUBROUTINE ORBOUT(O.T.XI.O.W.XN) 

IMPLICIT REAL.8 (A-H, 0-z) 

C 

C WRITES TEMPERATURES AND HEAT LOAD 
C 10 "RAD-ORB.DAT" 

C 

DIMENSION 0(6), 1(6) 

RTO = 57.29577951 
X = XI .RID 
XO = 0 >RTD 
WO = W .RTO 
XNU = XN <RTO 

WRITEflO,.)' INC)', X 
WRITE( 10,*) ' OMEGA]'. XO 
WRITE(IO..)' WOMEGA]', WO 
WRITE(IO.<) ’ NUO]'. XNU 
WRITEOO,*) 

WRI1E(I0..)' I1J\1(1) 

WRITEflO,*) ' T2]’.l(2) 

WRITEf 10.*) ' 13]'.T(3) 

WRITEflO.*)' T 4 j'.T(4 ) 

WRITEflO,.) ' T5]',T(5) 

WRITEflO.*) ' |T6j'.T(6) 

WRITEflO,*) 

DO 101=1,6 

IE (ABS(Ofl)) IT. I.OE-6) THEN 
0 ( 1 ) = 0.0 
ENOIF 

10 CONTINUE 

WRITEflO.*) ' |0l]',0(l) 

WRITEflO,.) ' 02 '.0(2) 

WRITEflO,.) ' Q3]'.Q(3) 

WRITEflO,*) ' 04]'.0(4) 

WRITEflO,.)' 05]', 0(5) 

WRITEflO,.)' 06]', Of6) 

WRITEflO.*) 

c 

RETURN 

END 



SUBROUTINE ORBWRITfTAVET.QAVET.IAVE.QAVE.IHIGH.TLOW.SHICH.SLOW. 
+ VOHIGH.VOLOW.QHIGH.QLOW.SQHIGH.SQLOW.VOQHIGH.VOQLOW.TDELTA. 
+ QDELTA.VODELT.VOQDEL) 

IMPLICIT REAL.8 (A-H. O-Z) 

DIMENSION TAVE(6).QAVE(6) 

C WRITES ALL THE FINAL STATISTICAL ANALYSIS TO RETURN.DAT 

OPEN( 1 1 ,FILE= 'RETURN.DAT-, STATUS= ’UNKNOWN’) 

RTD = 57.29577951 

WRITE(I t.«) ' |TAVET]'.TAVET 
WRITEfl 1 ,*) ‘ |OAVET]'.OAVET 
WRITE(I i.*) 

WRITEfl I.*) ' THIGHj'.THIGH 
WRITEfl l.») ' SHIGHj'.SHIGH 
WRITEfl 1.0 ’ VOHIGHj'.VOHIGH'RTO 
WRITEfl I.*) 

WRITEfl I.*)' TLOWj'.TLOW 
WRITEfl 1.0' SLOWj'.SLOW 
WRITEfl I.*) 1 VOLOWj'.VOLOW.RTD 
WRITEfl 1.0 

WRITEfl 1.0 ' |OHIGH]',OHIGH 
WRITEfl I.*) ' jsOHIGII j’.SOHIGH 
WRIlEft I.*) ' j VOOHIGI 1 1‘ .VOOHIGH • RTD 
WRITEfl I.*) 

WRITEfl 1.0' OLOWj'.OLOW 
WRITEfl I.*)' SOLOwj'.SOLOW 
WRITEfl 1.0 ' VOOLOW]’.VOOLOWtRTO 
WRITEfl 1.0 

WRITEfl 1,0 ' |TDELTA]',TDELTA 
WRITEfl 1.0 ’ VODELT]',VODELI«RTD 
WRITEfl 1.0' OOELTA]',ODELTA 
WRITEfl 1.0’ VOO DEL ] '.VOODE L » R1 D 
WRITEfl 1.0 

WRITEfl 1.0' TI]'.TAVE(l) 

WRITEfl 1.0 ' T2]'.TAVEf2) 

WRITEfl I.*) ' T3J',TAVE(3) 


WRITEfl 1.0' T4]',TAVE(4) 

WRITEfl 1.0 ‘ T5]’.1AVE(5) 

WRITEfl I.*)* T6 j'.TAVE(6) 

WRITEfl 1.0 

DO 10 1=1.6 

IF (ABS(OAVEfl)) IT. I.OE-6) THEN 
OAVE(I) = 0.0 


ENDIF 


10 CONTINUE 


WRITEfl 1.0 

01 J'.OAVEf 1) 

WRITEfl 1,0 

Q2]'.QAVE(2) 

WRITEfl 1.0 

03|'.0AVE(3) 

WRIIEfll.0 

04 '.0AVLf4) 

WRITEfl 1.0 

' 05]',0AVE(5) 

WRITEfl 1,0 

' 06]',0AVEf6) 

WRITEfl 1.0 


c 

C WRITE TO DATA FILE 


WRITEfl 0,0 ' TAVET J’.TAVET 
WRITEfl 0,0' OAVETj'.OAVET 
WRITEflO.O 

WRITEflO.O ' 1AVE 1 ]',TAVE( I ) 
WRITEflO.O ' IAVE2]'.TAVEf2) 
WRITEflO.O ' IAVE3]'.1AVE{3) 
WRITEflO.O ' !AVE4]',!AVE(4) 
WRITEflO.O ' TAVE5]',TAVEf5) 
WRITEflO.O' TAVEej'.TAVtfe) 
WRITEflO.O 


WRITEflO.O' OAVE 1 J'.OAVEf 1 ) 
WRITEflO.O ' QAVE2]',QAVE(2) 
WRITEflO.O ' |QAVE3]'.QAVEf3) 
WRITEflO.O ' OAVE4]'.OAVEf4) 
WRITEflO.O ' OAVE5]'.OAVEf5) 
WRITEflO.O ' 0AVE6 ,QAVEf6 
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WRIIEflO,*) • 
WRIIEflO,*) ‘ 
WRIIEflO,*)’ 

wriie(io,*) • 

WRIIEflO,*) 


VOHIGH]',VOHICH*RID 
U0W]'.110W 
SLOW]' .SLOW 
V010W|',V010W*R1D 


WRIIEflO.*) ' 
WRIIEflO.*) ' 
WRIIEflO.*) • 
WRIIEflO,*) ' 
WRIIEflO,*) ' 
WRIIEflO,*) • 
WRIIEflO,*) 


OHIGHJ'.OHIGH 

SOHIGH j'.SOKIGH 

VOOHIGH]'.VOOHIGH*RlD 

OLOW)',OLOW 

SOLOwj'.SOlOW 

VOOLOW]',VOOLOW.RIO 


WRIIEflO,*) 

WRIIEflO,*) 

WRIIEflO,*) 

WRIIEflO,*) 


I0ELTA)',10ELIA 
VODElllWODEll 
ODEUAj'.OOELIA 
VOODEL ] 1 .VOOOEL • R1 D 


CLOSEft 1) 

RETURN 

END 



ORBPROPS.FOR 


DOTPROD.FOR 


SUBROUTINE ORBPROPS(ECC.SMA.XNU.R) 

IMPLICIT REAl«8 (A-H. O-Z) 

C CALCULATES THE ORBITAL RADIUS 
C TAKEN FROM "FUNDAMENTALS OF ASTRODYNAMICS" BY 
C BATE. MUELLER AND WHITE 

XN = XNU » 0.0 1 7453293 

SLR = SMA.(I-ECC'ECC) 

R = SLR/(l+ECC*COS(XN)) 

RETURN 

END 


FUNCTION DOTPROD(X.Y) 

IMPLICIT REAL»8 (A-H. O-Z) 

C CALCULATES THE DOT PRODUCT OF TWO VECTORS X AND Y 
DIMENSION Y(3).X(J) 

Z = 0 

DO 10 1=1.3 
Z = X(I).Y(I) + Z 
10 CONTINUE 

DOTPROD = Z 

RETURN 

END 
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FETANOR.FOR 


PROGRAM FETANOR 
IMPLICIT REA1.8 (A-H, 0-Z) 

C THIS PROGRAM CALULATES THE EARTH VIEW FACIORS 
C NORMAL ANO TANGENTIAL PLATES. IT IS USED IN 
C CONJUCTION WITH IHE EXPERT SYSTEM RAOIAIE.RUL 
C RACHEL ORMSBY July 29. 1993 
C 

C READS THE INPUT OATA 
C ECC ECCENTRICITY 
C SMA SEMI- MAJOR AXIS (KM) 

C XNU TRUE ANOMALY AT EPOCH (RAD) 

CALL FEREAD(ECC.SMA.XNU.RE) 

C CALCULATES THE EARTH VIEW FACTOR FOR 
C A PLANE NORMAL TO THE EARTH 

CALL CALVFNOR(ECC,SMA,XNU,RE.F[N) 

C CALCULATES THE EARTH VIEW FACTOR FOR 
C A PLANE TANGENTIAL TO THE EARTH 
CALL CALVFTAN(ECC.SMA.XNU.RE.FET) 


SUBROUTINE FEREADfECC.SMA.XNU.RE ) 
IMPLICIT REAL *8 (A-H. O-Z) 

C 

C READS THE INPUT 
C ECC ECCENTRICITY 
C SMA SEMI- MAJOR AXIS 
C 

OPEN(3,FILE= 'RADIATE OAT ',STATUS= 'OL 0 ') 

DTR = 1.7453293E-2 

READ(3.*) ECC 
READ(3,*) SMA 
READ(3,0 XNU 
XNU = -XNU'DTR 
READ(3,«) RE 

CLOSE(3) 

RETURN 

END 


CALL FEWRIT(FEl.FEN) 
END 



o <~> o 


SUBROUTINE FEWRIT(FET.FEN) 

IMPLICIT REAL'8 (A-H, O-Z) 

WRITES NORMAL AND TANGENGTIAL 
VIEW FACTORS TO "REIURN.DAI" 

OPEN( 1 0.FILE = •RETURN.OAT.STATUS= ’UNKNOWN') 

WRITE(lO,») '|FET]'.FET 
WRITEflO,*) '|FEN]'.FEN 

CLOSE(IO) 

RETURN 

END 
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SUBROUTINE UK2PGR(XI.0.W.XNAB) 

IMPLICIT REAL»8 (A-H. 0-Z) 

C 

C •»•• THIS SUBROUTINE IS CURRENTLY NOT USED. BUT 
C LEFT HERE FOR FUTURE REFERENCE .... 

C 

C THIS SUBROUTINE FILLS THE COEFFICIENT MATRIX 
C THIS CONVERSION WAS TAKEN FROM "FUNDAMENTALS 
C OF AST ROOYAN M ICS" BY BATE. MUELLER AND WHITE 
C 

CXI) INCLINATION (RAO) 

C 01 LONGITUDE OF ASCENDING NODE (RAD) 

C Wl LINE OF PERIAPSIS (RAD) 

C XNU TRUE ANOMOLY AT EPOCH (RAD) 

C A TRANSFORMATION MATRIX FROM GLOBAL TO PERIFOCAL 
C B TRANSFORMATION MATRIX FROM PERIFOCAL TO 
C SPACECRAFT 

O DIMENSION A(3.3), 8(3,3) 

C>3 SXI = SIN(XI) 

SW = SIN(W) 

SO = SIN(O) 

SN = SIN(XN) 

CXI = COS(XI) 

CW = COS(W) 

CO = COS(O) 

CN = COS(XN) 

A(l.1) = sxi+sw 

A(l,2) = CXI.CXI.CW.SO 4 CW.SXI.SXI.SO 4 
CXI.CO.SW 

A( 1 .3) = -(CXKXI.CW.CO 4 CW.SXI.SXI.CO - 
CXI.SO.SW) 

A(2,l) = CW.SXI 

A(2.2) = CXI.CO.CW - CXKXI.SO.SW - 
SXI.SXI.SO«SW 

A(2.3) = -(-CXI.CW.SO - CXI.CXI.CO«SW - 
CO.SXI.SXI.SW) 

A(3,l) = -CXI 
A(3,2) = -(-CO.SXI) 

A(3,3) = SXI.SO 
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B(l.l) = CN 
B(1,2) = SN 
B(1.3) = 0 
B(2.l) = -SN 
B(2,2) = CN 
B(2,3) = 0 
B(3,l) = 0 
B(3,2) = 0 
6(3,3) = 1 

RETURN 

END 



MAG.FOR 


FUNCTION XMAC(X) 

IMPLICIT REAL *8 (A-H. O-Z) 

C CALCULATES THE MAGNITUDE OE THE VECTOR X 

DIMENSION X(3) 

A = X(1)*X(I) 

8 = X(2)»X(2) 

C = X(3)«X(3) 

D = A+B+C 

XMAG = SORT(D) 

RETURN 

END 
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SHADOW.FOR 


SUBROUTINE SHADOW(ECC.XNU.XINC.OMEGA,WOMEGA.GS.SORIG.S.XMU RE) 
C DETERMINES IE THE SPACECRAFT IS IN THE SHADOW Or IHE EARTH 
C ir IN IHE SHADOWS: 

C RE I URNS IHE SOLAR FLUX AS 0.0 
C AND IHE SOLAR VIEW FACTORS AS 0.0 
C IF NOT 

C RETURNS THE ORIGINAL VALUE FOR SOLAR FLUX 
C 

C ECC ECCENTRICITY 
C SMA SEMIrMAJOR AXIS (KM) 

C XNU I RUE ANOMALY Al EPOCH (DEG) 

C XINC INCLINATION (DEG) 

C OMEGA ORBIT LONGITUDE OF ASCENDING NODE (DEG) 

C WOMEGA LINE OF PERIAPSIS (DEG) 

C GS(3) GLOBAL SUN VECTOR 
C SORIG ORIGINAL NON ZERO SOLAR FLUX (W/m-2) 

C S CURRENT SOLAR FLUX (W/m-2) 

C = SORIG NOT IN SHADOW 
C 0 IN SHADOW 

C XMU(6) SOLAR VIEW FACTORS 
C RE RADIUS OF THE PLANET 

C 

IMPLICIT REAL.8 (A-H, 0-Z) 

DIMENSION GS(3), 0RIGIN(3) 

DIMENSION V(3), XMU(6) 

DIMENSION P0S(3) 

C CALCULATES THE CURRENT SPACECRAFT POSITION IN LOCAL COORDINATES 
CALL POSCLOB(ECC.SMA.XINC.OMEGA.WOMEGA,XNU.POS) 

C THE ORIGIN IS CURRENTLY DEFINED AT (0,0.0) 

ORIGIN(l) = 0.0 
0RIGIN(2) = 0.0 
0RIGIN(3) = 0.0 

C MATADO SUBTRACTS THE ORIGIN FROM THE POS VECTOR GIVING 
C A CORRECTED POSITION (V) 

CALL MATADD(3. 1 .POS.ORIGIN.V - 1 ) 

C DISTANCE FROM THE POSITION 10 THE VERTICAL SHADOW PLANE 
Dt = DOTPROD(V.GS) 
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C XMAG CALCULATES THE MAGNITUDE OF A VECTOR 
XMAGV = XMAG(V) 

C DETERMINES IF THE SPACECRAFT FALLS WITHIN THE SHADOW 
C DETERMINED BY THE CYLINDRICAL SHADOW PLANE 
D2 = SQRT(XMAGV»XMAGV - DI'Dl) 

IF (D2.LT.RE .AND. DI.L1.0) THEN 
S = 0.0 
DO 101=1,6 
XMU(l) = 0.0 
10 CONTINUE 

fi sr 

S = SORIG 

ENDIF 

RETURN 

END 
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